Télécharger cmaxot.eso

Retour à la liste

Numérotation des lignes :

  1. C CMAXOT SOURCE BP208322 17/03/01 21:15:43 9325
  2. SUBROUTINE CMAXOT(wrk52,wrk53,wrk54,WR12,IB,IGAU,NCHAIN,DT,NWA,
  3. & EPSFLU)
  4. *
  5. *
  6. *==============================================================
  7. * modele maxott
  8. *==============================================================
  9. *
  10. * entrees
  11. *
  12. * WRK12
  13. * EMi(2,NWA(i+1)) i=0 a 8
  14. * = EVOLUTION DU MODULE D'YOUNG DE LA BRANCHE i
  15. * SMi(NSTRS) i=1 a 8
  16. * = VARIABLES INTERNES AU DEBUT DU PAS D'INTEGRATION
  17. *
  18. * IVALMA(NCXMAT) = INDEX CARACTERISTIQUES MECANIQUES DU MATERIAU
  19. *
  20. *
  21. * IB = NUMERO DE L'ELEMENT
  22. * IGAU = NUMERO DU POINT DE GAUSS
  23. * NCHAIN = NOMBRE DE BRANCHE DU MODELE DE MAXWELL
  24. * DT = INCREMENT DE TEMPS
  25. * NWA = DIMENSION DES EVOLUTIONS EMi
  26. * EPSFLU = DEFORMATION DE FLUAGE
  27. *
  28. * sorties
  29. *
  30. * SIGF(NSTRS)= CONTR. A LA FIN DU PAS D'INTEGRATION
  31. * DEFP(NSTRS)= INCREMENT DES DEFORM. PLASTIQUES A LA FIN DU PAS
  32. * D'INTEGRATION
  33. * EPINF(NSTRS) = DEFORMATION INELASTIQUE A LA FIN DU PAS D'INTEGRATION
  34. *
  35. *==============================================================
  36. *
  37. IMPLICIT INTEGER(I-N)
  38. IMPLICIT REAL*8(A-H,O-Z)
  39. -INC CCOPTIO
  40. -INC DECHE
  41. *
  42. *
  43. SEGMENT WR12
  44. REAL*8 EM0(2,NWA(1)),EM1(2,NWA(2)),EM2(2,NWA(3))
  45. REAL*8 EM3(2,NWA(4)),EM4(2,NWA(5)),EM5(2,NWA(6))
  46. REAL*8 EM6(2,NWA(7)),EM7(2,NWA(8)),EM8(2,NWA(9))
  47. REAL*8 SM0(NSTRS),SM1(NSTRS),SM2(NSTRS),SM3(NSTRS)
  48. REAL*8 SM4(NSTRS),SM5(NSTRS),SM6(NSTRS),SM7(NSTRS)
  49. REAL*8 SM8(NSTRS)
  50. ENDSEGMENT
  51. *
  52. *
  53. DIMENSION NWA(9)
  54. DIMENSION EPSFLU(8)
  55. DIMENSION EPSOTT(8)
  56. DIMENSION DEPFIS(8)
  57. DIMENSION DEPSF(8)
  58. DIMENSION SFLU(8)
  59. DIMENSION SIGM(8)
  60. REAL*8 ETILDE,COE
  61. INTEGER IND
  62.  
  63. NCXMAT = XMAT(/1)
  64. NSTRS = SIG0(/1)
  65. NVARI = VAR0(/1)
  66. ICARA=XCARB(/1)
  67. ANU=XMAT(2)
  68. UNANU=1.D0-ANU
  69. TPS1=TEMP0
  70. TPS2=TEMP0+DT
  71. SE1=0.D0
  72. SE2=0.D0
  73. ETILDE = 0.D0
  74.  
  75. IF (IFOUR.EQ.2.AND.MFR.EQ.1) THEN
  76. IND = 21
  77. ELSE IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-3
  78. & .OR.IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
  79. IND = 16
  80. ELSE
  81. IND = 13
  82. ENDIF
  83.  
  84. *======================================================================
  85. * Formulation TRIDIMENSIONNELLE MASSIVE
  86. * =====================================
  87. * Formulation DEFORMATIONS PLANES /AXISYMETRIQUE
  88. * ==============================================
  89. *----------------------------------------------------------------------
  90. *
  91. IF ((IFOUR.EQ.2.AND.MFR.EQ.1).OR.
  92. & (IFOMOD.EQ.-1.AND.IFOUR.NE.-2).OR.
  93. & (IFOMOD.EQ.0.OR.IFOMOD.EQ.1)) THEN
  94. CONS=1.D0/((1.D0+ANU)*(1.D0-2.D0*ANU))
  95. DO M=0,NCHAIN-1
  96. *
  97. *----------------------------------------------------------------------
  98. * module d'YOUNG apparent
  99. *----------------------------------------------------------------------
  100.  
  101. CALL CFICMO(TPS1,TPS2,M,WR12,wrk52,wrk53,NWA,EFIC,TR,E1,E2)
  102.  
  103. ETILDE = ETILDE + EFIC
  104.  
  105. SE1=SE1+E1
  106. SE2=SE2+E2
  107.  
  108. IF (E1.EQ.0.) THEN
  109. COE = 1.
  110. ELSE
  111. COE = E2 / E1
  112. ENDIF
  113. *
  114. *----------------------------------------------------------------------
  115. * deformations de fluage
  116. *----------------------------------------------------------------------
  117. *
  118. FLUA=EXP(-TR*DT)
  119. DO II=1,NSTRS
  120. IF (M.EQ.0) THEN
  121. SFLU(II)=SM0(II)*(COE*FLUA-1.0D0)
  122. ELSE IF (M.EQ.1) THEN
  123. SFLU(II)=SFLU(II)+SM1(II)*(COE*FLUA-1.0D0)
  124. ELSE IF (M.EQ.2) THEN
  125. SFLU(II)=SFLU(II)+SM2(II)*(COE*FLUA-1.0D0)
  126. ELSE IF (M.EQ.3) THEN
  127. SFLU(II)=SFLU(II)+SM3(II)*(COE*FLUA-1.0D0)
  128. ELSE IF (M.EQ.4) THEN
  129. SFLU(II)=SFLU(II)+SM4(II)*(COE*FLUA-1.0D0)
  130. ELSE IF (M.EQ.5) THEN
  131. SFLU(II)=SFLU(II)+SM5(II)*(COE*FLUA-1.0D0)
  132. ELSE IF (M.EQ.6) THEN
  133. SFLU(II)=SFLU(II)+SM6(II)*(COE*FLUA-1.0D0)
  134. ELSE IF (M.EQ.7) THEN
  135. SFLU(II)=SFLU(II)+SM7(II)*(COE*FLUA-1.0D0)
  136. ELSE IF (M.EQ.8) THEN
  137. SFLU(II)=SFLU(II)+SM8(II)*(COE*FLUA-1.0D0)
  138. ENDIF
  139. END DO
  140. END DO
  141. *
  142. EPSFLU(1)=(SFLU(1) - ANU*SFLU(2) - ANU*SFLU(3)) / ETILDE
  143. EPSFLU(2)=(SFLU(2) - ANU*SFLU(3) - ANU*SFLU(1)) / ETILDE
  144. EPSFLU(3)=(SFLU(3) - ANU*SFLU(1) - ANU*SFLU(2)) / ETILDE
  145. DO II=4,NSTRS
  146. EPSFLU(II)=(SFLU(II)*(1.0D0+ANU))/ETILDE
  147. END DO
  148. *
  149. *----------------------------------------------------------------------
  150. * stockage du module d'elasticite
  151. *----------------------------------------------------------------------
  152. *
  153. VARF(IND) = ETILDE
  154. XMAT(1) = ETILDE
  155. *
  156. *----------------------------------------------------------------------
  157. * deformation corrigee pour Ottosen
  158. *----------------------------------------------------------------------
  159. *
  160. DO I=1,NSTRS
  161. EPSOTT(I) = DEPST(I) + EPSFLU(I)
  162. END DO
  163. *
  164. *----------------------------------------------------------------------
  165. * appel a Ottosen
  166. *----------------------------------------------------------------------
  167.  
  168. *
  169. CALL OTTOSE(INPLAS,SIG0,NSTRS,EPSOTT,VAR0,XMAT,IVALMA,NCXMAT,
  170. & XCARB,ICARA,NVARI,SIGF,VARF,DEPFIS,MFR,KERRE,IB,IGAU)
  171. *
  172. *----------------------------------------------------------------------
  173. * deformation corrigee pour Maxwell
  174. *----------------------------------------------------------------------
  175. *
  176. DO II=1,NSTRS
  177. DEPSF(II) = DEPST(II) - DEPFIS(II)
  178. END DO
  179. *
  180. *----------------------------------------------------------------------
  181. * modele de Maxwell
  182. *----------------------------------------------------------------------
  183. *
  184. DO M=0,NCHAIN-1
  185. CALL CFICMO(TPS1,TPS2,M,WR12,wrk52,wrk53,NWA,EFIC,TR,E1,E2)
  186. IF (E1.EQ.0.) THEN
  187. COE = 1.
  188. ELSE
  189. COE = E2 / E1
  190. ENDIF
  191. FLUA=EXP(-TR*DT)
  192. IF (M.EQ.0) THEN
  193. SM0(1)=SM0(1)*COE*FLUA
  194. 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  195. SM0(2)=SM0(2)*COE*FLUA
  196. 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  197. SM0(3)=SM0(3)*COE*FLUA
  198. 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  199. DO II = 4,NSTRS
  200. SM0(II)=SM0(II)*COE*FLUA
  201. 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  202. END DO
  203. DO II=1,NSTRS
  204. SIGM(II)=SM0(II)
  205. END DO
  206. ELSE IF (M.EQ.1) THEN
  207. SM1(1)=SM1(1)*COE*FLUA
  208. 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  209. SM1(2)=SM1(2)*COE*FLUA
  210. 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  211. SM1(3)=SM1(3)*COE*FLUA
  212. 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  213. DO II = 4,NSTRS
  214. SM1(II)=SM1(II)*COE*FLUA
  215. 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  216. END DO
  217. DO II=1,NSTRS
  218. SIGM(II)=SIGM(II)+SM1(II)
  219. END DO
  220. ELSE IF (M.EQ.2) THEN
  221. SM2(1)=SM2(1)*COE*FLUA
  222. 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  223. SM2(2)=SM2(2)*COE*FLUA
  224. 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  225. SM2(3)=SM2(3)*COE*FLUA
  226. 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  227. DO II = 4,NSTRS
  228. SM2(II)=SM2(II)*COE*FLUA
  229. 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  230. END DO
  231. DO II=1,NSTRS
  232. SIGM(II)=SIGM(II)+SM2(II)
  233. END DO
  234. ELSE IF (M.EQ.3) THEN
  235. SM3(1)=SM3(1)*COE*FLUA
  236. 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  237. SM3(2)=SM3(2)*COE*FLUA
  238. 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  239. SM3(3)=SM3(3)*COE*FLUA
  240. 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  241. DO II = 4,NSTRS
  242. SM3(II)=SM3(II)*COE*FLUA
  243. 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  244. END DO
  245. DO II=1,NSTRS
  246. SIGM(II)=SIGM(II)+SM3(II)
  247. END DO
  248. ELSE IF (M.EQ.4) THEN
  249. SM4(1)=SM4(1)*COE*FLUA
  250. 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  251. SM4(2)=SM4(2)*COE*FLUA
  252. 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  253. SM4(3)=SM4(3)*COE*FLUA
  254. 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  255. DO II = 4,NSTRS
  256. SM4(II)=SM4(II)*COE*FLUA
  257. 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  258. END DO
  259. DO II=1,NSTRS
  260. SIGM(II)=SIGM(II)+SM4(II)
  261. END DO
  262. ELSE IF (M.EQ.5) THEN
  263. SM5(1)=SM5(1)*COE*FLUA
  264. 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  265. SM5(2)=SM5(2)*COE*FLUA
  266. 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  267. SM5(3)=SM5(3)*COE*FLUA
  268. 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  269. DO II = 4,NSTRS
  270. SM5(II)=SM5(II)*COE*FLUA
  271. 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  272. END DO
  273. DO II=1,NSTRS
  274. SIGM(II)=SIGM(II)+SM5(II)
  275. END DO
  276. ELSE IF (M.EQ.6) THEN
  277. SM6(1)=SM6(1)*COE*FLUA
  278. 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  279. SM6(2)=SM6(2)*COE*FLUA
  280. 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  281. SM6(3)=SM6(3)*COE*FLUA
  282. 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  283. DO II = 4,NSTRS
  284. SM6(II)=SM6(II)*COE*FLUA
  285. 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  286. END DO
  287. DO II=1,NSTRS
  288. SIGM(II)=SIGM(II)+SM6(II)
  289. END DO
  290. ELSE IF (M.EQ.7) THEN
  291. SM7(1)=SM7(1)*COE*FLUA
  292. 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  293. SM7(2)=SM7(2)*COE*FLUA
  294. 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  295. SM7(3)=SM7(3)*COE*FLUA
  296. 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  297. DO II = 4,NSTRS
  298. SM7(II)=SM7(II)*COE*FLUA
  299. 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  300. END DO
  301. DO II=1,NSTRS
  302. SIGM(II)=SIGM(II)+SM7(II)
  303. END DO
  304. ELSE IF (M.EQ.8) THEN
  305. SM8(1)=SM8(1)*COE*FLUA
  306. 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  307. SM8(2)=SM8(2)*COE*FLUA
  308. 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  309. SM8(3)=SM8(3)*COE*FLUA
  310. 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  311. DO II = 4,NSTRS
  312. SM8(II)=SM8(II)*COE*FLUA
  313. 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  314. END DO
  315. DO II=1,NSTRS
  316. SIGM(II)=SIGM(II)+SM8(II)
  317. END DO
  318. ENDIF
  319. END DO
  320. DEFP(1)=DEPSF(1)-((SIGM(1)-ANU*(SIGM(2)+SIGM(3)))/SE2
  321. 1 -(SIG0(1)-ANU*(SIG0(2)+SIG0(3)))/SE1)
  322. DEFP(2)=DEPSF(2)-((SIGM(2)-ANU*(SIGM(1)+SIGM(3)))/SE2
  323. 1 -(SIG0(2)-ANU*(SIG0(1)+SIG0(3)))/SE1)
  324. DEFP(3)=DEPSF(3)-((SIGM(3)-ANU*(SIGM(2)+SIGM(1)))/SE2
  325. 1 -(SIG0(3)-ANU*(SIG0(2)+SIG0(1)))/SE1)
  326. DO II = 4,NSTRS
  327. DEFP(II)=DEPSF(II)-(1+ANU)*(SIGM(II)/SE2-SIG0(II)/SE1)
  328. END DO
  329. DO I=1,NSTRS
  330. EPINF(I)=EPIN0(I)+DEFP(I)+DEPFIS(I)
  331. END DO
  332. *
  333. *======================================================================
  334. * Formulation CONTRAINTES PLANES
  335. * ==============================
  336. *----------------------------------------------------------------------
  337. *
  338. ELSE IF (IFOMOD.EQ.-1.AND.IFOUR.EQ.-2) THEN
  339. CONS=1.D0/(1.D0-(ANU*ANU))
  340. DO M=0,NCHAIN-1
  341. *
  342. *----------------------------------------------------------------------
  343. * module d'YOUNG apparent
  344. *----------------------------------------------------------------------
  345.  
  346. CALL CFICMO(TPS1,TPS2,M,WR12,wrk52,wrk53,NWA,EFIC,TR,E1,E2)
  347. ETILDE = ETILDE + EFIC
  348. SE1=SE1+E1
  349. SE2=SE2+E2
  350. IF (E1.EQ.0.) THEN
  351. COE = 1.
  352. ELSE
  353. COE = E2 / E1
  354. ENDIF
  355. *
  356. *----------------------------------------------------------------------
  357. * deformations de fluage
  358. *----------------------------------------------------------------------
  359. *
  360. FLUA=EXP(-TR*DT)
  361. DO II=1,NSTRS
  362. IF (M.EQ.0) THEN
  363. SFLU(II)=SM0(II)*(COE*FLUA-1.0D0)
  364. ELSE IF (M.EQ.1) THEN
  365. SFLU(II)=SFLU(II)+SM1(II)*(COE*FLUA-1.0D0)
  366. ELSE IF (M.EQ.2) THEN
  367. SFLU(II)=SFLU(II)+SM2(II)*(COE*FLUA-1.0D0)
  368. ELSE IF (M.EQ.3) THEN
  369. SFLU(II)=SFLU(II)+SM3(II)*(COE*FLUA-1.0D0)
  370. ELSE IF (M.EQ.4) THEN
  371. SFLU(II)=SFLU(II)+SM4(II)*(COE*FLUA-1.0D0)
  372. ELSE IF (M.EQ.5) THEN
  373. SFLU(II)=SFLU(II)+SM5(II)*(COE*FLUA-1.0D0)
  374. ELSE IF (M.EQ.6) THEN
  375. SFLU(II)=SFLU(II)+SM6(II)*(COE*FLUA-1.0D0)
  376. ELSE IF (M.EQ.7) THEN
  377. SFLU(II)=SFLU(II)+SM7(II)*(COE*FLUA-1.0D0)
  378. ELSE IF (M.EQ.8) THEN
  379. SFLU(II)=SFLU(II)+SM8(II)*(COE*FLUA-1.0D0)
  380. ENDIF
  381. END DO
  382. END DO
  383. *
  384. EPSFLU(1)=(SFLU(1) - ANU*SFLU(2) - ANU*SFLU(3)) / ETILDE
  385. EPSFLU(2)=(SFLU(2) - ANU*SFLU(3) - ANU*SFLU(1)) / ETILDE
  386. EPSFLU(3)=(SFLU(3) - ANU*SFLU(1) - ANU*SFLU(2)) / ETILDE
  387. EPSFLU(4)=(1.0D0+ANU) * SFLU(4) / ETILDE
  388. *
  389. *----------------------------------------------------------------------
  390. * stockage du module d'elasticite
  391. *----------------------------------------------------------------------
  392. *
  393. VARF(IND) = ETILDE
  394. XMAT(1) = ETILDE
  395. *
  396. *----------------------------------------------------------------------
  397. * deformation corrigee pour Ottosen
  398. *----------------------------------------------------------------------
  399. *
  400. DO I=1,NSTRS
  401. EPSOTT(I) = DEPST(I) + EPSFLU(I)
  402. END DO
  403. *
  404. *----------------------------------------------------------------------
  405. * appel a Ottosen
  406. *----------------------------------------------------------------------
  407. *
  408.  
  409. CALL OTTOSE(INPLAS,SIG0,NSTRS,EPSOTT,VAR0,XMAT,IVALMA,NCXMAT,
  410. & XCARB,ICARA,NVARI,SIGF,VARF,DEPFIS,MFR,KERRE,IB,IGAU)
  411. *
  412. *----------------------------------------------------------------------
  413. * deformation corrigee pour Maxwell
  414. *----------------------------------------------------------------------
  415. *
  416. DO II=1,NSTRS
  417. DEPSF(II) = DEPST(II) - DEPFIS(II)
  418. END DO
  419. *----------------------------------------------------------------------
  420. * modele de Maxwell
  421. *----------------------------------------------------------------------
  422. DO M=0,NCHAIN-1
  423. CALL CFICMO(TPS1,TPS2,M,WR12,wrk52,wrk53,NWA,EFIC,TR,E1,E2)
  424. IF (E1.EQ.0.) THEN
  425. COE = 1.
  426. ELSE
  427. COE = E2 / E1
  428. ENDIF
  429. FLUA=EXP(-TR*DT)
  430. IF (M.EQ.0) THEN
  431. SM0(1)=SM0(1)*COE*FLUA
  432. 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  433. SM0(2)=SM0(2)*COE*FLUA
  434. 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  435. SM0(3)=0.D0
  436. SM0(4)=SM0(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  437. DO II=1,NSTRS
  438. SIGM(II)=SM0(II)
  439. END DO
  440. ELSE IF (M.EQ.1) THEN
  441. SM1(1)=SM1(1)*COE*FLUA
  442. 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  443. SM1(2)=SM1(2)*COE*FLUA
  444. 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  445. SM1(3)=0.D0
  446. SM1(4)=SM1(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  447. DO II=1,NSTRS
  448. SIGM(II)=SIGM(II)+SM1(II)
  449. END DO
  450. ELSE IF (M.EQ.2) THEN
  451. SM2(1)=SM2(1)*COE*FLUA
  452. 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  453. SM2(2)=SM2(2)*COE*FLUA
  454. 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  455. SM2(3)=0.D0
  456. SM2(4)=SM2(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  457. DO II=1,NSTRS
  458. SIGM(II)=SIGM(II)+SM2(II)
  459. END DO
  460. ELSE IF (M.EQ.3) THEN
  461. SM3(1)=SM3(1)*COE*FLUA
  462. 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  463. SM3(2)=SM3(2)*COE*FLUA
  464. 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  465. SM3(3)=0.D0
  466. SM3(4)=SM3(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  467. DO II=1,NSTRS
  468. SIGM(II)=SIGM(II)+SM3(II)
  469. END DO
  470. ELSE IF (M.EQ.4) THEN
  471. SM4(1)=SM4(1)*COE*FLUA
  472. 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  473. SM4(2)=SM4(2)*COE*FLUA
  474. 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  475. SM4(3)=0.D0
  476. SM4(4)=SM4(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  477. DO II=1,NSTRS
  478. SIGM(II)=SIGM(II)+SM4(II)
  479. END DO
  480. ELSE IF (M.EQ.5) THEN
  481. SM5(1)=SM5(1)*COE*FLUA
  482. 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  483. SM5(2)=SM5(2)*COE*FLUA
  484. 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  485. SM5(3)=0.D0
  486. SM5(4)=SM5(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  487. DO II=1,NSTRS
  488. SIGM(II)=SIGM(II)+SM5(II)
  489. END DO
  490. ELSE IF (M.EQ.6) THEN
  491. SM6(1)=SM6(1)*COE*FLUA
  492. 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  493. SM6(2)=SM6(2)*COE*FLUA
  494. 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  495. SM6(3)=0.D0
  496. SM6(4)=SM6(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  497. DO II=1,NSTRS
  498. SIGM(II)=SIGM(II)+SM6(II)
  499. END DO
  500. ELSE IF (M.EQ.7) THEN
  501. SM7(1)=SM7(1)*COE*FLUA
  502. 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  503. SM7(2)=SM7(2)*COE*FLUA
  504. 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  505. SM7(3)=0.D0
  506. SM7(4)=SM7(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  507. DO II=1,NSTRS
  508. SIGM(II)=SIGM(II)+SM7(II)
  509. END DO
  510. ELSE IF (M.EQ.8) THEN
  511. SM8(1)=SM8(1)*COE*FLUA
  512. 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  513. SM8(2)=SM8(2)*COE*FLUA
  514. 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  515. SM8(3)=0.D0
  516. SM8(4)=SM8(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  517. DO II=1,NSTRS
  518. SIGM(II)=SIGM(II)+SM8(II)
  519. END DO
  520. END IF
  521. END DO
  522. DEFP(1)=DEPSF(1)-((SIGM(1)-ANU*(SIGM(2)+SIGM(3)))/SE2
  523. 1 -(SIG0(1)-ANU*(SIG0(2)+SIG0(3)))/SE1)
  524. DEFP(2)=DEPSF(2)-((SIGM(2)-ANU*(SIGM(1)+SIGM(3)))/SE2
  525. 1 -(SIG0(2)-ANU*(SIG0(1)+SIG0(3)))/SE1)
  526. DEFP(3)=DEPSF(3)-((SIGM(3)-ANU*(SIGM(2)+SIGM(1)))/SE2
  527. 1 -(SIG0(3)-ANU*(SIG0(2)+SIG0(1)))/SE1)
  528. DEFP(4)=DEPSF(4)-(1+ANU)*(SIGM(4)/SE2-SIG0(4)/SE1)
  529. DO I=1,NSTRS
  530. EPINF(I)=EPIN0(I)+DEFP(I)+DEPFIS(I)
  531. END DO
  532. ENDIF
  533.  
  534. *======================================================================
  535. * QUELQUES IMPRESSIONS
  536. *======================================================================
  537. *
  538. IF (IIMPI.EQ.103) THEN
  539. WRITE(*,66770) IB,IGAU
  540. WRITE(*,*) dt,tps2
  541. WRITE(*,66775) XMAT(1)
  542. WRITE(*,66775) (XMAT(I),I=1,NCXMAT)
  543. WRITE(*,66776) (XCARB(I),I=1,ICARA)
  544. WRITE(*,66773) (VAR0(I),I=1,NVARI)
  545. WRITE(*,66774) (DEPST(I),I=1,NSTRS)
  546. WRITE(*,66777) (EPSFLU(I),I=1,NSTRS)
  547. WRITE(*,66778) (EPSOTT(I),I=1,NSTRS)
  548. WRITE(*,66779) (DEPFIS(I),I=1,NSTRS)
  549. WRITE(*,66771) (DEPSF(I),I=1,NSTRS)
  550. WRITE(*,66794) (DEFP(I),I=1,NSTRS)
  551. WRITE(*,66784) (SIGM(I),I=1,NSTRS)
  552. WRITE(*,66772) (SIG0(I),I=1,NSTRS)
  553. WRITE(*,66783) (SIGF(I),I=1,NSTRS)
  554. WRITE(*,66795) (EPIN0(I),I=1,NSTRS)
  555. WRITE(*,66782) (EPINF(I),I=1,NSTRS)
  556. ENDIF
  557. 66770 FORMAT(////////2X,'ELEMENT ',I6,2X,'POINT ',I3//)
  558. 66771 FORMAT(2X,' DEPSF '/(6(1X,1PE12.5)))
  559. 66772 FORMAT(2X,' SIG0 '/(6(1X,1PE12.5)))
  560. 66773 FORMAT(2X,' VAR0 '/(6(1X,1PE12.5)))
  561. 66774 FORMAT(2X,' DEPST '/(6(1X,1PE12.5)))
  562. 66775 FORMAT(2X,' XMAT '/(6(1X,1PE12.5)))
  563. 66776 FORMAT(2X,' XCARB '/(6(1X,1PE12.5)))
  564. 66777 FORMAT(2X,' EPSFLU '/(6(1X,1PE12.5)))
  565. 66778 FORMAT(2X,' EPSOTT '/(6(1X,1PE12.5)))
  566. 66779 FORMAT(2X,' DEPFIS '/(6(1X,1PE12.5)))
  567. 66782 FORMAT(2X,' EPINF '/(6(1X,1PE12.5)))
  568. 66783 FORMAT(2X,' SIGF '/(6(1X,1PE12.5)))
  569. 66784 FORMAT(2X,' SIGM '/(6(1X,1PE12.5)))
  570. 66794 FORMAT(2X,' DEFP '/(6(1X,1PE12.5)))
  571. 66795 FORMAT(2X,' EPIN0 '/(6(1X,1PE12.5)))
  572.  
  573. *----------------------------------------------------------------------
  574. RETURN
  575. END
  576.  
  577.  
  578.  
  579.  
  580.  
  581.  
  582.  
  583.  
  584.  
  585.  
  586.  
  587.  
  588.  
  589.  
  590.  
  591.  
  592.  
  593.  
  594.  
  595.  
  596.  
  597.  
  598.  
  599.  
  600.  
  601.  
  602.  
  603.  
  604.  
  605.  
  606.  
  607.  

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