Télécharger betaxi.eso

Retour à la liste

Numérotation des lignes :

  1. C BETAXI SOURCE CHAT 05/01/12 21:38:40 5004
  2. SUBROUTINE BETAXI(SIG0,NSTRS,DSIGT,VAR0,XMAT,IVAL,NCOMAT,XCAR,
  3. 1 SIGF,VARF,DEFP,MFR,KERRE,ECOU,NECOU)
  4. C----------------------------------------------------------------------
  5. C PLASTICITE MODELE BETON
  6. C
  7. C ENTREES
  8. C SIG0(NSTRS) = CONTRAINTES INITIALES
  9. C NSTRS = NOMBRE DE CONTRAINTES
  10. C DSIGT(NSTRS) = INCREMENT DE CONTRAINTES TOTALES
  11. C VAR0(NVARI) = VARIABLES INTERNES DEBUT
  12. C VAR0( 1 ) = EPSE
  13. C VAR0( 2 ) = DAM1
  14. C VAR0( 3 ) = DAM2
  15. C VAR0( 4 ) = DAM3
  16. C VAR0( 5 ) = ANGL
  17. C VAR0( 6 ) = LAM1
  18. C VAR0( 7 ) = LAM2
  19. C VAR0( 8 ) = LAM3
  20. C VAR0( 9 ) = LAM4
  21. C VAR0( 10 ) = LAM5
  22. C XMAT(NCOMAT) = COMPOSANTES DE MATERIAU
  23. C IVAL(NCOMAT) = INDICE DES COMPOSANTES DE MATERIAU
  24. C NCOMAT = NOMBRE DE COMPOSANTES DE MATERIAU
  25. C XCAR(ICARA) = CARACTERISTIQUES
  26. C MFR = NUMERO DE LA FORMULATION DE L'ELEMENT FINI
  27. C = 1 MASSIF
  28. C = 3 COQUE MINCE ( COQ2 , COQ3 ET DKT )
  29. C = 5 COQUE EPAISSE ( COQ6 , COQ8 )
  30. C = 7 POUTRE
  31. C = 9 COQUE MINCE AVEC CISAILLEMENT TRANSVERSE
  32. C ( COQ4 ET DST )
  33. C = 11 LIQUIDE
  34. C = 13 TUYAU
  35. C = 15 LINESPRING
  36. C = 17 TUYAU FISSURE
  37. C = 19 RACCORD MASSIF
  38. C = 21 RACCORD COQUE
  39. C = 23 SURFACE LIBRE
  40. C = 25 MEMBRANE
  41. C = 27 UNIAXIALE
  42. C = 29 THERMIQUE
  43. C = 31 INCOMPRESSIBLES
  44. C = 33 POREUX
  45. C = 35 JOINT
  46. C = 37 HOMOGENEISE
  47. C = 39 TUYO
  48. C = 41 TUYAU ACOUSTIQUE PURE
  49. C = 43 RACCORD TUYAU FLUIDE
  50. C SORTIES
  51. C SIGF(NSTRS) = CONTRAINTES FINALES
  52. C VARF(NVARI) = VARIABLES INTERNES FINALES
  53. C DEFP(NSTRS) = DEFORMATIONS PLASTIQUES
  54. C KERRE = 0 TOUT OK
  55. C 1 SI DLAMBDA NEGATIF
  56. C 2 NOMBRE MAX D ITERATIONS INTERNES DEPASSE
  57. C 21 ON NE TROUVE PAS L INTERSECTION AVEC LA SURFACE DE CHARGE
  58. C 22 SIG0 A L EXTERIEUR DE LA SURFACE DE CHARGE
  59. C
  60. C-----------------------------------------------------------------------
  61. C VARIABLES PASSEES PAR LES COMMONS COPTIO , ECOU ET NECOU
  62. C
  63. C IFOUR INDICE DU TYPE DE PROBLEME
  64. C -3 DEFORMATIONS PLANES GENERALISEES
  65. C -2 CONTRAINTES PLANES
  66. C -1 DEFORMATIONS PLANES
  67. C 0 AXISYMETRIQUE
  68. C 1 SERIE DE FOURIER
  69. C 2 TRIDIMENSIONNEL
  70. C ITYP TYPE DE FORMULATION MECANIQUE
  71. C --------------- ATTENTION ---------------
  72. C IL EST ACTIF APRES L APPEL DE VISAVI
  73. C -----------------------------------------
  74. C ITYP=1 CAS DES ELEMENTS MASSIFS
  75. C ITYP=2 CAS DES COQUES
  76. C ITYP=3 CAS DES MEMBRANES
  77. C ITYP=4 CAS DES CABLES ET DES BARRES
  78. C ITYP=5 CAS QUELCONQUE
  79. C ITYP=6 CAS DES CONTRAINTES PLANES
  80. C ITYP=7 CAS DES COQUES A NU=0. OU CONTRAINTES PLANES
  81. C ITYP=8 CAS DES MEMBRANES A NU=0. OU CONTRAINTES PLANES
  82. C ITYP=9 CAS DES COQUES EPAISSES
  83. C ITYP=10 CAS DES JOINTS
  84. C ITYP=11 CAS DES POUTRES
  85. C ITYP=12 CAS DES TUYAUX
  86. C ITYP=13 CAS DES COQUES AVEC CISAILLEMENT TRANSVERSE
  87. C
  88. C-----------------------------------------------------------------------
  89. C
  90. IMPLICIT INTEGER(I-N)
  91. IMPLICIT REAL*8(A-H,O-Z)
  92.  
  93. -INC PPARAM
  94. -INC CCOPTIO
  95.  
  96. SEGMENT ECOU
  97. *** COMMON/ECOU/TEST,ALFAH,
  98. REAL*8 TEST, ALFAH,
  99. 1 ecow1,ecow2,ecow3(6),ecow4(9),ecow5(6),
  100. C 1 HPAS,TEMPS,COVNMS(6),VECPRO(9),VALPRO(6),
  101. 2 ecow6(12),ecow7(6),ecow8(6),ecow9(6),ecow10(6),ecow11(6),
  102. 2 ecow12(6),
  103. C 2 CVNMSD(12),STOT(6),SIGEL(6),DSIGP(6),SIGT(6),W1(6),W2(6),
  104. 1 ecow13(6),ecow14(6),ecow15(12),ecow16(3),
  105. C 1 DALPHA(6),EPSPLA(6),E(12),XINV(3),
  106. 2 ecow17(6),ecow18(6),ecow19,ecow20
  107. C 2 SIPLAD(6),DSIGP0(6),TET,TETI
  108. ENDSEGMENT
  109.  
  110. SEGMENT NECOU
  111. * COMMON/NECOU/NCOURB,IPLAST,IT,IMAPLA,ISOTRO,
  112. * INTEGER ncow1, ncow2,ncow3,ncow4, ncow5,
  113. INTEGER NCOURB,IPLAST,IT, IMAPLA,ISOTRO,
  114. 1 ITYP,ncow7, ncow8,
  115. C . ITYP, IFOURB,IFLUAG,
  116. 2 ncow9,ncow10,ncow11,ncow12,ncow13,
  117. C . ICINE,ITHER, IFLUPL,ICYCL, IBI,
  118. 3 ncow14,ncow15,ncow16,
  119. C . JFLUAG,KFLUAG,LFLUAG,
  120. 4 ncow17,ncow18,ncow19,ncow20,ncow21
  121. C . IRELAX,JNTRIN,MFLUAG,JSOUFL,JGRDEF
  122. ENDSEGMENT
  123.  
  124.  
  125. * COMMON/ECOU/TEST,ALFAH,ZOZO(112)
  126. * COMMON/NECOU/NCOURB,IPLAST,IT,IMAPLA,ISOTRO,ITYP,IZOZO(15)
  127. C
  128. DIMENSION SIG0(*),DSIGT(*),VAR0(*),XMAT(*),XCAR(*),IVAL(*)
  129. DIMENSION SIGF(*),VARF(*),DEFP(*)
  130. DIMENSION SIGMA(6),DSIGMA(6),STOT(6),SIGFIN(6)
  131. DIMENSION XLAMBD(5),DEFPLA(6),DSIGP(6)
  132. DIMENSION SPHER(6),AUXIL(6),DSIGZE(6),ORMAT(1)
  133. CHARACTER*8 CMATE
  134. CMATE='ISOTROPE'
  135. C
  136. C TEST DE CONSISTANCE DES DONNEES
  137. C
  138. DO 100 I=1,NCOMAT
  139. IF(XMAT(I).LT.0.D0) KERRE=437
  140. 100 CONTINUE
  141. IF(KERRE.EQ.437) RETURN
  142. C
  143. YOUN= 0.D0
  144. XNU = 0.D0
  145. RHO = 0.D0
  146. ALPH= 0.D0
  147. CTR1= 0.D0
  148. CTR2= 0.D0
  149. CTR3= 0.D0
  150. ANGL= 0.D0
  151. CLCS= 0.D0
  152. EPCS= 0.D0
  153. CLCT= 0.D0
  154. PCT = 0.D0
  155. EPCT= 0.D0
  156. CLCD= 0.D0
  157. PCD = 0.D0
  158. EPCD= 0.D0
  159. CLPH= 0.D0
  160. PENT= 0.D0
  161. C
  162. YOUN= XMAT(1)
  163. XNU = XMAT(2)
  164. RHO = XMAT(3)
  165. ALPH= XMAT(4)
  166. CTR1= XMAT(5)
  167. CTR2= XMAT(6)
  168. CTR3= XMAT(7)
  169. ANGL= XMAT(8)
  170. CLCS= XMAT(9)
  171. EPCS= XMAT(10)
  172. CLCT= XMAT(11)
  173. PCT = XMAT(12)
  174. CLCD= XMAT(13)
  175. PCD = XMAT(14)
  176. EPCD= XMAT(15)
  177. CLPH= XMAT(16)
  178. PENT= XMAT(17)
  179. C
  180. IF(IIMPI.EQ.9) WRITE(IOIMP,2007) (XMAT(I),I=1,17)
  181. IF(IIMPI.EQ.9) WRITE(IOIMP,2008) (IVAL(I),I=1,17)
  182. C
  183. C
  184. C PETIT TEST SUR NU POUR CERTAINS CAS
  185. C
  186. IF(MFR.EQ.1.AND.(IFOUR.GE.-1.OR.IFOUR.EQ.-3)) THEN
  187. C
  188. C VALEURS PAR DEFAUT
  189. C
  190. C DONNEES GENERALES
  191. C
  192. IF(IVAL( 3).EQ.0) RHO =2.4D+3
  193. IF(IVAL( 4).EQ.0) ALPH=1.2D-5
  194. C
  195. C POUR LA TRACTION
  196. C
  197. IF(IVAL( 5).EQ.0) CTR1=YOUN*1.2E-4
  198. IF(IVAL( 6).EQ.0) CTR2=CTR1
  199. IF(IVAL( 7).EQ.0) CTR3=CTR1
  200. C
  201. C POUR LA COMPRESSION
  202. C
  203. IF(IVAL( 9).EQ.0) CLCS=YOUN*1.2D-3
  204. IF(IVAL(10).EQ.0) EPCS=10.D0*CLCS/YOUN
  205. IF(IVAL(12).EQ.0) PCT=1.6D0*CLCS
  206. C
  207. C BI-COMPRESSION
  208. C
  209. REP=EPCS*YOUN/CLCS
  210. AAA=CLCT/CLCS
  211. BBB=PCT/CLCS
  212. CCC=CLCT/PCT
  213. UMN=1.D0-XNU
  214. ALF=((REP-UMN-(SQRT(0.5D0*REP*UMN)))/(2.D0*REP-UMN))*0.95D0
  215. DDD=(1.D0-ALF)/(1.D0-2.D0*ALF)
  216. CLBI=DDD*CLCS
  217. C
  218. C CISAILLEMENT FRAGILE
  219. C
  220. IF(IVAL(11).EQ.0) THEN
  221. CCC=((1.D0+BBB)/BBB)+((3.D0*ALF)/(1.D0-ALF))
  222. CLCT=CCC*PCT
  223. ELSE IF(IVAL(12).EQ.0) THEN
  224. BBB=((AAA-1.D0)*(1.D0-ALF))/(1.D0+2.D0*ALF)
  225. PCT=BBB*CLCS
  226. ENDIF
  227. C
  228. C CISAILLEMENT DUCTILE
  229. C
  230. IF(IVAL(14).EQ.0) PCD=1.125D0*PCT
  231. CCC=CLCT/PCT
  232. EEE=PCD/PCT
  233. FFF=(((1.D0+2.D0*ALF)*(EEE-1.D0))/(EEE*(1.D0-ALF)))+(CCC/EEE)
  234. FFF=FFF/1.05D0
  235. IF(IVAL(13).EQ.0) CLCD=FFF*PCD
  236. IF(IVAL(15).EQ.0) EPCD=5.D0*EPCS
  237. C
  238. C PRESSION HYDROSTATIQUE
  239. C
  240. BBB=PCT/CLCS
  241. IF(IVAL(16).EQ.0) CLPH=(((CCC+2.D0)*BBB)/3.D0)*1.05D0*CLCS
  242. IF(IVAL(17).EQ.0) PENT=YOUN/(9.D0*(1.D0-2.D0*XNU))
  243. C
  244. C INITIALISATION
  245. C
  246. JEPRIN = 0
  247. PREC=1.D-2
  248. RFPR=1.D-3
  249. RFEP=1.D-7
  250. RFSG=YOUN*RFEP
  251. C
  252. C CALCUL ET TEST DES PARAMETRES DU BETON
  253. C
  254. C CALCUL DES PENTES
  255. C
  256. XK1=YOUN/(3.D0*(1.D0-2.D0*XNU))
  257. XK2=PENT
  258. EPSA=CLCS/YOUN
  259. EPSB=(CLCT-PCT)/YOUN
  260. EPSC=(CLCD-PCD)/YOUN
  261. PENTE1=(CLCD-PCD-CLCT+PCT)/(EPCD-EPSB)
  262. PENTE1=(YOUN*PENTE1)/(YOUN-PENTE1)
  263. PENTE2=-CLCS/(EPCS-EPSA)
  264. PENTE2=(YOUN*PENTE2)/(YOUN-PENTE2)
  265. PENTE3=XK1*XK2/(XK1-XK2)
  266. C
  267. CLCS=-CLCS
  268. CLCSA=CLCS*RFEP
  269. CLCT=-CLCT
  270. CLCD=-CLCD
  271. PCT=-PCT
  272. PCD=-PCD
  273. EPSA=-EPSA
  274. EPSB=-EPSB
  275. EPSC=-EPSC
  276. EPCD=-EPCD
  277. EPCS=-EPCS
  278. C
  279. C CALCUL DES CONSTANTES
  280. C
  281. XI1A=CLCS
  282. XI1AA=CLCSA
  283. XI1B=CLCT+2.D0*PCT
  284. XI1C=CLCD+2.D0*PCD
  285. SEQA=ABS(CLCS)
  286. SEQAA=ABS(CLCSA)
  287. SEQB=ABS(CLCT-PCT)
  288. SEQC=ABS(CLCD-PCD)
  289. C
  290. VMELAS=SEQB
  291. C
  292. ALFADV=(SEQC-SEQB)/(XI1B-XI1C)
  293. DPELAS=ALFADV*XI1B+SEQB
  294. C
  295. ALFAD1=(SEQB-SEQAA)/(XI1AA-XI1B)
  296. DPELA1=ALFAD1*XI1B+SEQB
  297. C
  298. ALFAD2=(SEQB-SEQA)/(XI1A-XI1B)
  299. DPELA2=ALFAD2*XI1B+SEQB
  300. PENTE2=PENTE2*(ALFAD2-1.D0)*(ALFAD2-1.D0)
  301. C
  302. PORELA=CLPH
  303. C
  304. UMU=YOUN/(1.D0+XNU)/2.D0
  305. ALAM=2.D0*UMU*XNU/(1.-2.D0*XNU)
  306. CPLA=YOUN/(1.D0-(XNU*XNU))
  307. C
  308. TRTRAC=CTR1+CTR2+CTR3
  309. XLAMER=(ALFAD2*TRTRAC-DPELA2)/PENTE2
  310. XLAMAX=(ALFAD2*DPELA1-ALFAD1*DPELA2)/(ALFAD1*PENTE2)
  311. DP2MIN=DPELA1*ALFAD2/ALFAD1
  312. C
  313. SIGT=TRTRAC*0.333333333333333333333333333333333D0
  314. PMIN=-XI1B*0.333333333333333333333333333333333D0
  315. SIG1=(DPELA2-(2.D0*ALFAD2+1.D0)*PCD)/(ALFAD2-1.D0)
  316. SIG3=(DPELA2-(ALFAD2-1.D0)*CLCD)/(2.D0*ALFAD2+1.D0)
  317. C
  318. C TEST SUR LES VALEURS DES PARAMETRES
  319. C
  320. IF(CLCT.GT.CLCS.OR.CLCT.LT.CLCD) THEN
  321. KERRE=640
  322. WRITE(IOIMP,5010) ABS(CLCT),ABS(CLCS),ABS(CLCD)
  323. ENDIF
  324. C
  325. IF(PCT.LT.PCD) THEN
  326. KERRE=640
  327. WRITE(IOIMP,5011) ABS(PCT),ABS(PCD)
  328. ENDIF
  329. C
  330. IF(EPCS.GE.EPSA) THEN
  331. KERRE=640
  332. WRITE(IOIMP,5012) ABS(EPCS),ABS(EPSA)
  333. ENDIF
  334. C
  335. IF(EPCD.GE.EPSC) THEN
  336. KERRE=640
  337. WRITE(IOIMP,5013)ABS(EPCD),ABS(EPSC)
  338. ENDIF
  339. C
  340. IF(ALFADV.LT.0.D0) THEN
  341. KERRE=640
  342. WRITE(IOIMP,5017)
  343. ENDIF
  344. C
  345. IF(ALFAD1.LT.0.D0.OR.ALFAD2.LT.0.D0) THEN
  346. KERRE=640
  347. WRITE(IOIMP,5018)
  348. ENDIF
  349. C
  350. IF(ABS(CLCS).LT.SIGT) THEN
  351. KERRE=640
  352. WRITE(IOIMP,5014) ABS(CLCS),SIGT
  353. ENDIF
  354. C
  355. IF(PMIN.GE.PORELA) THEN
  356. KERRE=640
  357. WRITE(IOIMP,5015) PMIN
  358. ENDIF
  359. C
  360. IF(ALFAD1.GE.0.5D0) THEN
  361. KERRE=640
  362. WRITE(IOIMP,5019) ABS(CLCT),(4*ABS(PCT))
  363. ENDIF
  364. C
  365. IF(DPELA2.GE.DPELAS) THEN
  366. KERRE=640
  367. WRITE(IOIMP,5020) ABS(CLCD),ABS(SIG1),ABS(PCD),ABS(SIG3)
  368. ENDIF
  369. C
  370. IF(XK2.GE.XK1) THEN
  371. KERRE=640
  372. WRITE(IOIMP,5021) XK2,XK1
  373. ENDIF
  374. C
  375. ECR2=(2.D0*YOUN*(0.5D0-ALFAD2)*(0.5D0-ALFAD2))/(1.D0-XNU)
  376. IF(ABS(PENTE2).GT.ECR2) THEN
  377. KERRE=640
  378. WRITE(IOIMP,5016) ABS(PENTE2),ECR2
  379. ENDIF
  380. C
  381. ECR3=(1.5D0/(1.D0+XNU))+(3.D0*ALFAD2*ALFAD2/(1.D0-2.D0*XNU))
  382. ECR3=ECR3*YOUN
  383. IF(ABS(PENTE2).GT.ECR3) THEN
  384. KERRE=640
  385. WRITE(IOIMP,5022) ABS(PENTE2),ECR3
  386. ENDIF
  387. C
  388. IF(KERRE.NE.0) RETURN
  389. C
  390. C ECOULEMENT PLASTIQUE
  391. C
  392. 333 MCOD=1
  393. CALL VISAVI(SIG0,DSIGT,VAR0,SIGMA,DSIGMA,SPHER,AUXIL,
  394. . SIGF,DEFP,VARF,SIGFIN,DEFPLA,
  395. . DSIGZE,ICENT2,MCOD,IBOU,MFR,NSTRS,XCAR,CMATE,ecou,necou)
  396. IF(ITYP.EQ.0) THEN
  397. KERRE=269
  398. RETURN
  399. ENDIF
  400. C
  401. C IMPRESSION DES PARAMETRES DU MODELE
  402. C
  403. IF(IIMPI.EQ.9) THEN
  404. WRITE(IOIMP,1000)
  405. WRITE(IOIMP,1001) YOUN
  406. WRITE(IOIMP,1002) XNU
  407. WRITE(IOIMP,1003) RHO
  408. WRITE(IOIMP,1004) ALPH
  409. WRITE(IOIMP,1005)
  410. WRITE(IOIMP,1006) CTR1,CTR2,CTR3
  411. WRITE(IOIMP,1007) ANGL
  412. WRITE(IOIMP,1008)
  413. WRITE(IOIMP,1009) ABS(CLCS),ABS(EPCS)
  414. WRITE(IOIMP,1010) ABS(CLBI)
  415. WRITE(IOIMP,1011)
  416. WRITE(IOIMP,1012) ABS(CLCT),ABS(PCT)
  417. WRITE(IOIMP,1013) ABS(CLCD),ABS(PCD),ABS(EPCD)
  418. WRITE(IOIMP,1014)
  419. WRITE(IOIMP,1015) CLPH,PENT
  420. ENDIF
  421. C
  422. C INITIALISATIONS DES PARAMETRES
  423. C
  424. EPSE=VAR0(1)
  425. IF(CTR1.GT.0.D0) RT1=(1.D0-VAR0(2))*CTR1
  426. IF(CTR1.LE.0.D0) RT1=0.D0
  427. IF(CTR2.GT.0.D0) RT2=(1.D0-VAR0(3))*CTR2
  428. IF(CTR2.LE.0.D0) RT2=0.D0
  429. IF(CTR3.GT.0.D0) RT3=(1.D0-VAR0(4))*CTR2
  430. IF(CTR3.LE.0.D0) RT3=0.D0
  431. ANGLE=VAR0(5)
  432. IF(IVAL( 8).NE.0) ANGLE=ANGL
  433. XLAMBD(1)=VAR0(6)
  434. XLAMBD(2)=VAR0(7)
  435. XLAMBD(3)=VAR0(8)
  436. XLAMBD(4)=VAR0(9)
  437. XLAMBD(5)=VAR0(10)
  438. C
  439. IFIS=0
  440. IF((RT1.EQ.0.D0.OR.RT2.EQ.0.D0).AND.ANGLE.NE.0.D0) IFIS=1
  441. IF(IVAL( 8).NE.0) IFIS=1
  442. C
  443. IDED=0
  444. IF(XLAMBD(3).GT.0.D0) IDED=IDED+1
  445. IF(XLAMBD(5).GT.0.D0) IDED=IDED+2
  446. C
  447. DO 200 IB=1,IBOU
  448. STOT(IB)=SIGMA(IB)+DSIGMA(IB)
  449. 200 CONTINUE
  450. C
  451. C ECOULEMENT PLASTIQUE
  452. C
  453. IF(IIMPI.EQ.9) THEN
  454. WRITE(IOIMP,2000)
  455. WRITE(IOIMP,2001) (SIGMA(I),I=1,IBOU)
  456. WRITE(IOIMP,2002) (DSIGMA(I),I=1,IBOU)
  457. WRITE(IOIMP,2003) (STOT(I),I=1,IBOU)
  458. WRITE(IOIMP,2006) (XLAMBD(I),I=1,5)
  459. WRITE(IOIMP,2010) RT1,RT2,RT3
  460. ENDIF
  461. C
  462. CALL DAMAGE(ALFADV,ALFAD1,ALFAD2,VMELAS,DPELAS,
  463. . DPELA1,DPELA2,PORELA,RT1,RT2,RT3,
  464. . PENTE1,PENTE2,PENTE3,XLAMBD,SIGMA,DSIGMA,STOT,
  465. . SIGFIN,ANGLE,YOUN,XNU,XLAMER,XLAMAX,DP2MIN,
  466. . PREC,RFSG,RFEP,RFPR,IFIS,IDED,IDAM,KERRE)
  467.  
  468. C
  469. IF(IIMPI.EQ.9) WRITE(IOIMP,2011) IDAM
  470. C
  471. IF(KERRE.NE.0) THEN
  472. IF(JEPRIN.EQ.0) THEN
  473. KERRE=0
  474. JEPRIN=1
  475. IIMPI=9
  476. GO TO 333
  477. ELSE
  478. IIMPI=0
  479. RETURN
  480. ENDIF
  481. ENDIF
  482. C
  483. C PAS DE PLASTICITE COMPORTEMENT ELASTIQUE
  484. C
  485. IF(IDAM.EQ.1) THEN
  486. DO 300 IB=1,IBOU
  487. SIGFIN(IB)=STOT(IB)
  488. DEFPLA(IB)=0.D0
  489. 300 CONTINUE
  490. DEPSEQ=0.D0
  491. C
  492. ELSE
  493. C
  494. C NOUS AVONS ENDOMMAGE LE MATERIAU
  495. C
  496. C CALCUL DE DEFORMATION PLASTIQUE
  497. C
  498. DO 400 IB=1,IBOU
  499. DSIGP(IB)=STOT(IB)-SIGFIN(IB)
  500. 400 CONTINUE
  501. CALL EPSIG(DSIGP,DEFPLA,IFOUR,YOUN,XNU,ITYP,ORMAT,XCAR)
  502. DEPSEQ=SQRT(AVM(DEFPLA,DEFPLA))
  503. C
  504. ENDIF
  505. C
  506. IF(IIMPI.EQ.9) THEN
  507. WRITE(IOIMP,2000)
  508. WRITE(IOIMP,2001) (SIGMA(I),I=1,IBOU)
  509. WRITE(IOIMP,2002) (DSIGMA(I),I=1,IBOU)
  510. WRITE(IOIMP,2003) (STOT(I),I=1,IBOU)
  511. WRITE(IOIMP,2004) (SIGFIN(I),I=1,IBOU)
  512. WRITE(IOIMP,2005) (DEFPLA(I),I=1,IBOU)
  513. WRITE(IOIMP,2006) (XLAMBD(I),I=1,5)
  514. ENDIF
  515. C
  516. C MISE A JOUR DES VARIABLES INTERNES
  517. C
  518. VARF(1)=EPSE+DEPSEQ
  519. IF(CTR1.GT.0.D0) VARF(2)=1.D0-RT1/CTR1
  520. IF(CTR1.LE.0.D0) VARF(2)=1.D0
  521. IF(CTR2.GT.0.D0) VARF(3)=1.D0-RT2/CTR2
  522. IF(CTR2.LE.0.D0) VARF(3)=1.D0
  523. IF(CTR3.GT.0.D0) VARF(4)=1.D0-RT3/CTR3
  524. IF(CTR3.LE.0.D0) VARF(4)=1.D0
  525. VARF(5)=ANGLE
  526. VARF(6)=XLAMBD(1)
  527. VARF(7)=XLAMBD(2)
  528. VARF(8)=XLAMBD(3)
  529. VARF(9)=XLAMBD(4)
  530. VARF(10)=XLAMBD(5)
  531. C
  532. IF(IIMPI.EQ.9) WRITE(IOIMP,2009) EPSE,DEPSEQ,VARF(1)
  533. C
  534. MCOD=2
  535. CALL VISAVI(SIG0,DSIGT,VAR0,SIGMA,DSIGMA,SPHER,AUXIL,
  536. . SIGF,DEFP,VARF,SIGFIN,DEFPLA,
  537. . DSIGZE,ICENT2,MCOD,IBOU,MFR,NSTRS,XCAR,CMATE,ecou,necou)
  538. C
  539. ELSE
  540. KERRE=639
  541. RETURN
  542. ENDIF
  543. C
  544. C
  545. C LES FORMATS D IMPRESSION
  546. C
  547. 1000 FORMAT(//,31X,'MATERIAU BETON',//)
  548. 1001 FORMAT(6X,'MODULE D YOUNG ',T51,'= ',1PD12.5)
  549. 1002 FORMAT(6X,'COEF. DE POISSON',T51,'= ',1PD12.5)
  550. 1003 FORMAT(6X,'MASSE VOLUMIQUE ',T51,'= ',1PD12.5)
  551. 1004 FORMAT(6X,'COEF. DE DILATATION THERMIQUE',T51,'= ',1PD12.5)
  552. 1005 FORMAT(/,1X,'DOMMAGE PAR TRACTION :',/)
  553. 1006 FORMAT(6X,'CONTRAINTES LIMITES :',/
  554. * 11X,'DIRECTION(1) = ',1PD12.5,
  555. * 5X,'DIRECTION(2) = ',1PD12.5,5X,'DIRECTION(3) = ',1PD12.5/)
  556. 1007 FORMAT(6X,'ANGLE DE FISSURATION',T51,'= ',1PD12.5/)
  557. 1008 FORMAT(/,1X,'DOMMAGE PAR COMPRESSION :',/)
  558. 1009 FORMAT(6X,'COMPRESSION SIMPLE :',/
  559. * 11X,'CONTRAINTE LIMITE = ',1PD12.5,
  560. * 5X,'DEFORMATION A RUPTURE = ',1PD12.5/)
  561. 1010 FORMAT(6X,'COMPRESSION BIAXIALE :',/
  562. * 11X,'CONTRAINTE LIMITE = ',1PD12.5/)
  563. 1011 FORMAT(/,1X,'DOMMAGE PAR CISAILLEMENT :',/)
  564. 1012 FORMAT(6X,'SOUS CONFINEMENT DE TRANSITION :',/
  565. * 11X,'CONTRAINTE LIMITE = ',1PD12.5,
  566. * 5X,'PRESSION DE CONFINEMENT = ',1PD12.5/)
  567. 1013 FORMAT(6X,'SOUS CONFINEMENT DUCTILE :',/
  568. * 11X,'CONTRAINTE LIMITE = ',1PD12.5,
  569. * 5X,'PRESSION DE CONFINEMENT = ',1PD12.5,/
  570. * 11X,'DEFORMATION AU DEBUT DU PALIER PLASTIQUE = ',1PD12.5/)
  571. 1014 FORMAT(/,1X,'DOMMAGE PAR PRESSION HYDROSTATIQUE :',/)
  572. 1015 FORMAT(6X,'PRESSION LIMITE = ',1PD12.5,
  573. * 12X,'MODULE TANGENT = ',1PD12.5,//)
  574. C
  575. 2000 FORMAT(1X,'RESULTATS DE L ECOULEMENT PLASTIQUE')
  576. 2001 FORMAT(1X,'SIGMA =',6(1X,1PD12.5))
  577. 2002 FORMAT(1X,'DSIGMA =',6(1X,1PD12.5))
  578. 2003 FORMAT(1X,'STOT =',6(1X,1PD12.5))
  579. 2004 FORMAT(1X,'SIGFIN =',6(1X,1PD12.5))
  580. 2005 FORMAT(1X,'DEFPLA =',6(1X,1PD12.5))
  581. 2006 FORMAT(1X,'XLAMBD =',5(1X,1PD12.5))
  582. 2007 FORMAT(7(6(1X,1PD12.5),/))
  583. 2008 FORMAT(7(6(1X,I12),/))
  584. 2009 FORMAT(1X,'EPSE0 =',1X,1PD12.5,1X,'DEPSEQ=',1X,1PD12.5,
  585. . 1X,'EPSE =',1X,1PD12.5)
  586. 2010 FORMAT(1X,'RT1 =',1X,1PD12.5,1X,'RT2 =',1X,1PD12.5,
  587. . 1X,'RT3 =',1X,1PD12.5)
  588. 2011 FORMAT(1X,'INDICE D ENDOMMAGEMENT IDAM =',I4)
  589. C
  590. 5010 FORMAT(1X,'ERREUR LA LIMITE LCT =',1PD12.5,' DOIT ETRE COMPRISE',
  591. . /,1X,'ENTRE LCS =',1PD12.5,' ET LCD =',1PD12.5)
  592. 5011 FORMAT(1X,'ERREUR LE CONFINEMENT PCT =',1PD12.5,' DOIT ETRE ',/,
  593. . 1X,'INFERIEUR A PCD =',1PD12.5)
  594. 5012 FORMAT(1X,'ERREUR LA DEFORMATION A RUPTURE EN COMPRESSION SIMPLE'
  595. . ,' EST MAUVAISE',/,
  596. . 1X,'EPCS =',1PD12.5,' DOIT ETRE SUPERIEURE A ',1PD12.5)
  597. 5013 FORMAT(1X,'ERREUR LA DEFORMATION A RUPTURE EN COMPRESSION DUCTILE'
  598. . ,' EST MAUVAISE',/,
  599. . 1X,'EPCD =',1PD12.5,' DOIT ETRE SUPERIEURE A ',1PD12.5)
  600. 5014 FORMAT(1X,'ERREUR LA LIMITE EN COMPRESSION SIMPLE LCS =',1PD12.5,
  601. . /,1X,'DOIT ETRE SUPERIREURE A LA LIMITE EN TRACTION LCT =',
  602. . D12.5)
  603. 5015 FORMAT(1X,'ERREUR LA LIMITE EN PRESSION HYDROSTATIQUE SE TROUVE',
  604. . ' DANS LE DOMAINE FRAGILE',/,
  605. . 1X,'IL FAUT AUGMENTER LPH POUR QU''ELLE SOIT SUPERIEURE ',
  606. . 'A PMIN =',1PD12.5)
  607. 5016 FORMAT(1X,'ERREUR MODULE D ECROUISSAGE DU DRUCKER EST TROP GRAND',
  608. . ' =',1PD12.5,/,
  609. . 1X,'PENTE2 < 2 * E * (0.5 - ALFAD2 ) ** 2 / (1-NU) ',
  610. . 'DONC PENTE2 MAX VAUT =',1PD12.5)
  611. 5017 FORMAT(1X,'ERREUR LA VALEUR DE ( LCD - PCD ) DOIT ETRE ',
  612. . 'SUPERIEURE A ( LCT - PCT )')
  613. 5018 FORMAT(1X,'ERREUR LA VALEUR DE ( LCT - PCT ) DOIT ETRE ',
  614. . 'SUPERIEURE A LCS')
  615. 5019 FORMAT(1X,'ERREUR LA LIMITE LCT =',1PD12.5,'DOIT ETRE INFERIEURE',
  616. . ' A',1PD12.5)
  617. 5020 FORMAT(1X,'ERREUR LA LIMITE LCD =',1PD12.5,'DOIT ETRE INFERIEURE',
  618. . ' A',1PD12.5,' OU BIEN :',/,
  619. . 1X,'LE CONFINEMENT PCD =',1PD12.5,' DOIT ETRE SUPERIEURE',
  620. . ' A',1PD12.5)
  621. 5021 FORMAT(1X,'ERREUR LA PENTE DE LA PRESSION HYDROSTATIQUE ',
  622. . 'EST MAUVAISE',/,
  623. . 1X,'PENT =',1PD12.5,' DOIT ETRE INFERIEURE A ',1PD12.5)
  624. 5022 FORMAT(1X,'ERREUR MODULE D ECROUISSAGE DU DRUCKER EST TROP GRAND',
  625. . ' =',1PD12.5,/,
  626. . 1X,'PENTE2 < 1.5*E/(1+NU) + 3*ALFAD2**2*E/(1-2*NU) ',
  627. . 'DONC PENTE2 MAX VAUT ',1PD12.5)
  628. C
  629. RETURN
  630. END
  631.  
  632.  
  633.  
  634.  
  635.  
  636.  
  637.  
  638.  
  639.  
  640.  

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