Télécharger beinsa.eso

Retour à la liste

Numérotation des lignes :

  1. C BEINSA SOURCE PV 05/04/25 21:15:01 5077
  2. C BEINSA SOURCE INSL 24/10/96
  3. SUBROUTINE BEINSA(SIG0,NSTRSS,DEPST,VAR0,XMAT,IVAL,NMATT,
  4. 1 SIGF,VARF,KERRE,MELE,IFOUR,NVARI,XCAR,NCARR,MFR,EPIN0,EPINF,
  5. 2 DTEMP,XE,NBNN,CMATE,IBB,IGAU,wrk12)
  6. C---------------------------------------------------------------------
  7. C PLASTICITE MODELE BETON
  8. C
  9. C ENTREES
  10. C SIG0(NSTRS) = CONTRAINTES INITIALES
  11. C NSTRS = NOMBRE DE CONTRAINTES
  12. C DEPST(NSTRS) = INCREMENT DE DEFORMATION TOTALES
  13. C VAR0(NVARI) = VARIABLES INTERNES DEBUT
  14. C VAR0( 1 ) = EPPE :Deformation plastique equivalente
  15. C VAR0( 2 ) = ALPT1 :Coef. multi. Defo. elas. au pic dir. 1
  16. C VAR0( 3 ) = ALPT2 :Coef. multi. Defo. elas. au pic dir. 2
  17. C VAR0( 4 ) = IPLA1 :Indicateur plas dir. 1
  18. C VAR0( 5 ) = IFISU2 :Indicateur de fissure actuelle (0 ou 1)
  19. C VAR0( 6 ) = DTR1 :Resistance en traction dir. 1
  20. C VAR0( 7 ) = DTR2 :Resistance en traction dir. 2
  21. C VAR0( 8 ) = IFIS0 :Indicateur de fissuration etat precedent
  22. C VAR0( 9 ) = JFISU :Indicateur de fissure (0 ou 1)
  23. C VAR0( 10 ) = ANGL :Angle de fissuration
  24. C VAR0( 11 ) = IPLA :Indicateur d'etat en bicompres. (0 2 4)
  25. C VAR0( 12 ) = EQSTRC1 :Contrainte equiv. de compres. dir. 1
  26. C VAR0( 13 ) = EPSEQC1 :Deformation equiv. de compres. dir. 1
  27. C VAR0( 14 ) = EQSTRT1 :Contrainte equiv. de traction dir. 1
  28. C VAR0( 15 ) = EPSEQT1 :Deformation equiv. de traction dir. 1
  29. C VAR0( 16 ) = EQSTRC2 :Contrainte equiv. de compres. dir. 2
  30. C VAR0( 17 ) = EPSEQC2 :Deformation equiv. de compres. dir. 2
  31. C VAR0( 18 ) = EQSTRT2 :Contrainte equiv. de traction dir. 2
  32. C VAR0( 19 ) = EPSEQT2 :Deformation equiv. de traction dir. 2
  33. C VAR0( 20 ) = IDECH1 :Indicateur de decharge beton intergre
  34. C VAR0( 21 ) = OUV1 :Indicateur ouverture fissure 1
  35. C VAR0( 22 ) = EDC1 :Module (d'elasti.) initiale equiva. dir.1
  36. C VAR0( 23 ) = EDC2 :Module (d'elasti.) initiale equiva. dir.2
  37. C VAR0( 24 ) = ETS1 :Pente courbe de traction (Post-pic) dir.1
  38. C VAR0( 25 ) = ETS2 :Pente courbe de traction (Post-pic) dir.2
  39. C VAR0( 26 ) = EPF01 :Deformation au pic courbe de traction dir.1
  40. C VAR0( 27 ) = EPF02 :Deformation au pic courbe de traction dir.2
  41. C VAR0( 28 ) = OUV2 :Indicateur ouverture fissure 2
  42. C VAR0( 29 ) = DEFR1 :defo residuelle en compres.1 (beton fissure)
  43. C VAR0( 30 ) = DEFR2 :defo residuelle en compres.2 (beton fissure)
  44. C VAR0( 31 ) = IPAS0 :Initial. variables internes 0=oui/1=non
  45. C VAR0( 32 ) = EPSR(1) :Deform. initi. EpsX (repere local elem.)
  46. C VAR0( 33 ) = EPSR(2) :Deform. initi. EpsY ...
  47. C VAR0( 34 ) = EPSR(3) :Deform. initi. GamXY ...
  48. C VAR0( 35 ) = EPSR(4) :Deform. initi. GamXZ ...
  49. C VAR0( 36 ) = EPSR(5) :Deform. initi. GamYZ ...
  50. C VAR0( 37 ) = DEP(1,1):Terme de matrice elastoplastique tangente
  51. C VAR0( 38 ) = DEP(1,2):
  52. C VAR0( 39 ) = DEP(1,3):
  53. C VAR0( 40 ) = DEP(1,4):
  54. C VAR0( 41 ) = DEP(1,5):
  55. C VAR0( 42 ) = DEP(2,1):
  56. C VAR0( 43 ) = DEP(2,2):
  57. C VAR0( 44 ) = DEP(2,3):
  58. C VAR0( 45 ) = DEP(2,4):
  59. C VAR0( 46 ) = DEP(2,5):
  60. C VAR0( 47 ) = DEP(3,1):
  61. C VAR0( 48 ) = DEP(3,2):
  62. C VAR0( 49 ) = DEP(3,3):
  63. C VAR0( 50 ) = DEP(3,4):
  64. C VAR0( 51 ) = DEP(3,5):
  65. C VAR0( 52 ) = DEP(4,1):
  66. C VAR0( 53 ) = DEP(4,2):
  67. C VAR0( 54 ) = DEP(4,3):
  68. C VAR0( 55 ) = DEP(4,4):
  69. C VAR0( 56 ) = DEP(4,5):
  70. C VAR0( 57 ) = DEP(5,1):
  71. C VAR0( 58 ) = DEP(5,2):
  72. C VAR0( 59 ) = DEP(5,3):
  73. C VAR0( 60 ) = DEP(5,4):
  74. C VAR0( 61 ) = DEP(5,5):
  75. C XMAT(NMATT) = COMPOSANTES DE MATERIAU
  76. C IVAL(NMATT) = INDICE DES COMPOSANTES DE MATERIAU
  77. C NMATT = NOMBRE DE COMPOSANTES DE MATERIAU
  78. C SORTIES
  79. C SIGF(NSTRS) = CONTRAINTES FINALES
  80. C VARF(NVARI) = VARIABLES INTERNES FINALES
  81. C DEFP(NSTRS) = DEFORMATIONS PLASTIQUES
  82. C KERRE = 0 TOUT OK
  83. C---------------------------------------------------------------------
  84. C VARIABLES PASSEES PAR LES COMMONS COPTIO , ECOU ET NECOU
  85. C
  86. C IFOUR INDICE DU TYPE DE PROBLEME
  87. C -2 CONTRAINTES PLANES
  88. C -1 DEFORMATIONS PLANES
  89. C 0 AXISYMETRIQUE
  90. C 1 SERIE DE FOURIER
  91. C 2 TRIDIMENSIONNEL
  92. C---------------------------------------------------------------------
  93. C COMPOSANTES DE MATERIAU
  94. C=====================================================================
  95. C EX : Module d'Young
  96. C XNU : Coeficient de Poisson
  97. C RHO : Masse volumique (Facultatif)
  98. C ALPH : Coeficient de dilation thermique (Facultatif)
  99. C ALPHA: Resis. tract. simple / resis. compr. simple
  100. C RB : Resis. compr. simple
  101. C EMAX : Deformation ultime en compression (sigma = 0)
  102. C EPUT : Deformation ultime en traction (sigma = 0)
  103. C FTC : Facteur de transfert de cisaillement
  104. C
  105. C ICAL : Choix du comportement Pre-pic
  106. C ICAL 0 = Comportement elastoplastique ecrouissable
  107. C 1 = Comportement elastique
  108. C
  109. C ICU : Choix de la courbe de compression (Calibrage)
  110. C 1 = Parabole pre-pic et post-pic
  111. C 2 = Parabole pre-pic et lineaire post-pic
  112. C 3 = courbe INSA pre-pic et lineaire post-pic
  113. C 5 = courbe de Krishnan pre-pic et post-pic
  114. C 4 = courbe de Krishnan pre-pic et lineaire post-pic
  115. C
  116. C ILOI : Choix de la loi d'ecoulement
  117. C 0 = Plasticite Associee
  118. C 1 = Plasticite Non Associee
  119. C
  120. C IOUV : Nombre d'iteration interne
  121. C
  122. C IMOD : Choix du modele Beton
  123. C 0 = Modele BETON_INSA ( critere d'OTTOSEN)
  124. C 1 = Modele BETON_INSA ( critere de NADAI partout)
  125. C 2 = Modele BETON_INSA ( NADAI + Sigma Max en trac.)
  126. C
  127. C---------------------------------------------------------------------
  128. C
  129. IMPLICIT INTEGER(I-N)
  130. IMPLICIT REAL*8(A-H,O-Z)
  131. CHARACTER*8 CMATE
  132. DIMENSION SIG0(NSTRSS),DEPST(NSTRSS),VAR0(NVARI),XMAT(NMATT)
  133. DIMENSION IVAL(*),SIGF(NSTRSS),VARF(NVARI),XE(3,NBNN)
  134. DIMENSION EPIN0(NSTRSS),EPINF(NSTRSS),DEPIN(6),EPFLU(6)
  135. DIMENSION STRN(6),SIGM(6),D(6,6),D1(6,6),D2(6,6),XCAR(NCARR)
  136. DIMENSION EPSR(6),SIGR(6),VART(100),VV(36)
  137.  
  138.  
  139. C
  140. segment wrk12
  141. real*8 AA,BB,DK1,DK2,RB,ALPHA,EX,XNU,EMAX
  142. real*8 EPUT,FTC,EPO,EPO1,ENGF,RMOY,PHIF,TEMP0
  143. real*8 DTEMP1,TEMP1,POAR,SCT,TETA,DTR1,DTR2,EDC1
  144. real*8 EDC2,ETS1,ETS2,EDT1,EDT2,OUV1,OUV2,TANG1
  145. real*8 TANG2,DEFR1,DEFR2,EPSC1,EPSC2,EPST1,EPST2,EQSTR1
  146. real*8 EQSTR2,EPSEQ1,EPSEQ2,EQSTR3,EPSEQ3,EPST3,EPSC3,DEFR3
  147. real*8 RTM3,EDC3,ETS3,EDT3,OUV3,TANG3
  148. integer ICU,ILOI,IOUV,ICAL,IFLU,IPLA2,IPLA1,IFISU2
  149. integer IFISU1,JFISU,JFISU2,IPLA3,IFISU3,JFISU3,IBB1,IGAU1
  150. endsegment
  151. c 1 2 3 4 5 6 7 8 9 10 11 12
  152. c COMMON /CINSA/ AA,BB,DK1,DK2,RB,ALPHA,EX,XNU,EMAX,EPUT,FTC,EPO,
  153. c 13 14 15 16 17 18 19 20 21 22
  154. c 1 EPO1,ENGF,RMOY,PHIF,TEMP0,DTEMP1,TEMP1,POAR,SCT,TETA,
  155. c 23 24 25 26 27 28 29 30 31 32 33
  156. c 2 DTR1,DTR2,EDC1,EDC2,ETS1,ETS2,EDT1,EDT2,OUV1,OUV2,TANG1,
  157. c 34 35 36 37 38 39 40 41 42
  158. c 3 TANG2,DEFR1,DEFR2,EPSC1,EPSC2,EPST1,EPST2,EQSTR1,EQSTR2,
  159. C 43 44 45 46 47 48 49 50 51
  160. c 4 EPSEQ1,EPSEQ2,EQSTR3,EPSEQ3,EPST3,EPSC3,DEFR3,RTM3,EDC3,
  161. c 52 53 54 55
  162. c 5 ETS3,EDT3,OUV3,TANG3,
  163. C 1 2 3 4 5 6 7 8 9
  164. c 6 ICU,ILOI,IOUV,ICAL,IFLU,IPLA2,IPLA1,IFISU2,IFISU1,
  165. C 10 11 12 13 14 15 16
  166. c 7 JFISU,JFISU2,IPLA3,IFISU3,JFISU3,IBB1,IGAU1
  167. *
  168. IRTD =1
  169. IBB1 =IBB
  170. IGAU1=IGAU
  171. C
  172. C TEST DE CONSISTANCE DES DONNEES
  173. C
  174. EX = XMAT( 1)
  175. XNU = XMAT( 2)
  176. RHO = XMAT( 3)
  177. ALPH = XMAT( 4)
  178. IKA=4
  179. IF((MFR.EQ.1.OR.MFR.EQ.3.OR.MFR.EQ.31).AND.IFOUR.EQ.-2) IKA=5
  180. *
  181. IFOU=IFOUR
  182. IF(MELE.EQ.28.OR.MFR.EQ.3) IFOU=-2
  183. C
  184. ALPHA=0.1D0
  185. RB =EX*1.D-3
  186. EMAX =10.D0*RB/EX
  187. EPUT =3.D0*RB*ALPHA/EX
  188. FTC =0.1D0
  189. ICU =3
  190. ILOI =0
  191. IOUV =2
  192. IFLU =0
  193. PHIF =2.D0
  194. TEMP1=0.D0
  195. POAR =0.D0
  196. ENGF =0.D0
  197. *
  198. ICAL = 0
  199. DTEMP1=DTEMP
  200. C
  201. AA=0.D0
  202. BB=0.D0
  203. DK1=0.D0
  204. DK2=0.D0
  205. *
  206. ALPHA= XMAT(IKA+1)
  207. RB = XMAT(IKA+2)
  208. EMAX = XMAT(IKA+3)
  209. EPUT = XMAT(IKA+4)
  210. FTC = XMAT(IKA+5)
  211. ICAL = INT(REAL(XMAT(IKA+6)))
  212. ENGF = XMAT(IKA+7)
  213. *
  214. IFLU = INT(REAL(XMAT(IKA+8)))
  215. RMOY = XMAT(IKA+9)
  216. PHIF = XMAT(IKA+10)
  217. TEMP1= XMAT(IKA+11)
  218. POAR = XMAT(IKA+12)
  219. EPO = 2.4D0*RB/EX
  220. C
  221. C VALEURS PAR DEFAUT
  222. C
  223. IF(MELE.EQ.4.OR.MELE.EQ.8.OR.MELE.EQ.6.OR.MELE.EQ.10) THEN
  224. IF(IVAL(IKA+7).EQ.0.AND.IVAL(IKA+4).EQ.0) ENGF =.15D0
  225. ELSE
  226. ENGF =.0D0
  227. ENDIF
  228. *
  229. IF(IVAL(IKA+8).EQ.0) IFLU=0
  230. IF(IVAL(IKA+10).EQ.0) PHIF=2.D0
  231. IF(IVAL(IKA+11).EQ.0) TEMP1=0.D0
  232. IF(IVAL(IKA+12).EQ.0) POAR =.0D0
  233. C---------------------------------------------------------------------
  234. IECR=0
  235. IF(IECR.EQ.10) THEN
  236. * WRITE(*,*) ' ** XMAT='
  237. * WRITE(*,1991) (XMAT(I),I=1,NMATT)
  238. * WRITE(*,*) ' ** IVAL='
  239. * WRITE(*,1992) (IVAL(I),I=1,NMATT)
  240. WRITE(*,404) EX,XNU,RHO,ALPH,ALPHA,RB,EMAX,EPUT,FTC,ICAL,ENGF
  241. 404 FORMAT('EX=',E9.3,' XNU=',E9.3,' RHO=',E9.3,' ALPH=',E9.3,
  242. *' ALPHA=',E9.3,/,' RB=',E9.3,' EMAX=',E9.3,' EPUT=',E9.3,
  243. &' FTC=',E9.3,' ICAL=',I1,' ENGF=',E9.3)
  244. WRITE(*,405) NSTRSS,NMATT,MELE,NVARI,NCARR,MFR,NBNN
  245. 405 FORMAT(' NSTRSS=',I1,' NMATT=',I2,' MELE=',I2,/,' NVARI=',I2,
  246. *' NCARR=',I2,' MFR=',I2,' NBNN=',I2)
  247. ENDIF
  248. C---------------------------------------------------------------------
  249. C
  250. IF(IFLU.EQ.20.OR.IFLU.EQ.21) THEN
  251. IF(IVAL(IKA+10).EQ.0) PHIF= 60.D0
  252. ENDIF
  253. *
  254. IF(IFLU.EQ.30.OR.IFLU.EQ.31) THEN
  255. IF(IVAL(IKA+9).EQ.0) RMOY= 0.159D0
  256. IF(IVAL(IKA+10).EQ.0) PHIF=-0.88D0
  257. ENDIF
  258. *
  259. IF((IFLU.EQ.10.OR.IFLU.EQ.11.OR.IFLU.EQ.20.OR.IFLU.EQ.21).AND.
  260. & IVAL(IKA+9).EQ.0) THEN
  261. WRITE(*,*)'!! ATTENTION POUR FAIRE DU FLUAGE BETON (LOI BPEL91)'
  262. WRITE(*,*)' IL FAUT DONNER LA VALEUR DU RAYON MOYEN (RMOY) '
  263. STOP
  264. ENDIF
  265. C---------------------------------------------------------------------
  266. C
  267. C CORRESPONDANCE DES VARIABLES POUR BETDJE
  268. C
  269. CALL ZERO(D,6,6)
  270. CALL ZERO(D1,6,6)
  271. CALL ZERO(EPSR,6,1)
  272. CALL ZERO(DEPIN,6,1)
  273. CALL ZERO(EPFLU,6,1)
  274. CALL ZERO(VART,100,1)
  275. CALL ZERO(SIGR,6,1)
  276. CALL ZERO(SIGM,6,1)
  277. CALL ZERO(STRN,6,1)
  278. C
  279. EPAIST=1.D0
  280. NSTRS=NSTRSS
  281. *
  282. IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.6) THEN
  283. NSTRS=3
  284. EPAIST=XCAR(1)
  285. ENDIF
  286. *
  287. DO 1 I=1,NSTRS
  288. SIGR(I)=SIG0(I)/EPAIST
  289. STRN(I)=DEPST(I)
  290. EPFLU(I)=EPIN0(I)
  291. 1 CONTINUE
  292. *
  293. IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.4) THEN
  294. NSTRS=3
  295. SIGR(3) =SIG0(4)
  296. STRN(3) =DEPST(4)
  297. ENDIF
  298. C---------------------------------------------------------------------
  299. IF(IECR.EQ.10) THEN
  300. WRITE(*,*) ' ** SIGR='
  301. WRITE(*,1991) (SIGR(I),I=1,NSTRS)
  302. WRITE(*,*) ' ** STRN='
  303. WRITE(*,1991) (STRN(I),I=1,NSTRS)
  304. ENDIF
  305. C------------------------------------------------------
  306. C INITIALISATION DES VARIABLES INTERNES
  307. C------------------------------------------------------
  308. DO II=1,NVARI
  309. VART(II)=VAR0(II)
  310. END DO
  311. C
  312. EPSPL = VART( 1)
  313. TANG1 = VART( 2)
  314. TANG2 = VART( 3)
  315. IPLA1 = INT(REAL(VART( 4)))
  316. IFISU2 = INT(REAL(VART( 5)))
  317. DTR1 = VART( 6)
  318. DTR2 = VART( 7)
  319. IFISU1 = INT(REAL(VART( 8)))
  320. JFISU = INT(REAL(VART( 9)))
  321. ANGL = VART(10)
  322. JFISU2 = INT(REAL(VART(11)))
  323. EQSTR1 = VART(12)
  324. EPSEQ1 = VART(13)
  325. EPST1 = VART(14)
  326. EPST2 = VART(15)
  327. EQSTR2 = VART(16)
  328. EPSEQ2 = VART(17)
  329. EPSC1 = VART(18)
  330. EPSC2 = VART(19)
  331. IPLA2 = INT(REAL(VART(20)))
  332. OUV1 = VART(21)
  333. EDC1 = VART(22)
  334. EDC2 = VART(23)
  335. ETS1 = VART(24)
  336. ETS2 = VART(25)
  337. EDT1 = VART(26)
  338. EDT2 = VART(27)
  339. OUV2 = VART(28)
  340. DEFR1 = VART(29)
  341. DEFR2 = VART(30)
  342. IPAS0 = INT(REAL(VART(31)))
  343. TEMP0 = VART(32)
  344. C
  345. KD = 32
  346. IF(NSTRS.EQ.4.AND.(IFOU.EQ.-1.OR.IFOU.EQ.0)) THEN
  347. EQSTR3 = VART(KD+1)
  348. EPSEQ3 = VART(KD+2)
  349. EPST3 = VART(KD+3)
  350. EPSC3 = VART(KD+4)
  351. DEFR3 = VART(KD+5)
  352. RTM3 = VART(KD+6)
  353. EDC3 = VART(KD+7)
  354. ETS3 = VART(KD+8)
  355. EDT3 = VART(KD+9)
  356. TANG3 = VART(KD+10)
  357. OUV3 = VART(KD+11)
  358. IPLA3 = INT(REAL(VART(KD+12)))
  359. IFISU3 = INT(REAL(VART(KD+13)))
  360. JFISU3 = INT(REAL(VART(KD+14)))
  361. KD = 46
  362. ENDIF
  363. C
  364. EPSR(1)= VART(KD+1)
  365. EPSR(2)= VART(KD+2)
  366. EPSR(3)= VART(KD+3)
  367. C -------------------------------------
  368. IF(IFLU.NE.0.AND.(TEMP0+DTEMP-TEMP1).LE.0.D0) THEN
  369. WRITE(*,*)'!! ATTENTION DANS LE MODELE DE FLUAGE BETON '
  370. TEMPF = TEMP0+DTEMP-TEMP1
  371. WRITE(*,*)' (LOI BPEL91) TEMPS NEGATIF TEMPF=',TEMPF
  372. WRITE(*,*)' TEMP0=',TEMP0,' TEMP1=',TEMP1,' DTEMP=',DTEMP
  373. STOP
  374. ENDIF
  375. C -------------------------------------
  376. TETA=ANGL
  377. IF((NSTRS.EQ.4.OR.NSTRS.EQ.6).AND.IFOU.EQ.-2) THEN
  378. NSTRS=3
  379. SIGR(3)=SIG0(4)
  380. STRN(3)=DEPST(4)
  381. STRN(4)=DEPST(3)
  382. SIGR(4)=0.D0
  383. EPFLU(3)=EPIN0(4)
  384. EPFLU(4)=EPIN0(3)
  385. ENDIF
  386. C
  387. GO TO (10,10,30,40,50,60),NSTRS
  388. 10 CONTINUE
  389. KERRE=437
  390. WRITE(*,*) '!! ATTENTION DANS BETON NSTRS=',NSTRS
  391. STOP
  392. c+mdj
  393. C GO TO 1000
  394. c+mdj
  395. C
  396. 60 CONTINUE
  397. C
  398. C COMPORTEMENT TRIDIMENSIONNEL
  399. C
  400. EPSR(4)= VART(KD+4)
  401. EPSR(5)= VART(KD+5)
  402. EPSR(6)= VART(KD+6)
  403. GO TO 30
  404. C
  405. 50 CONTINUE
  406. C
  407. C CONTRAINTES OU DEFORMATIONS PLANES
  408. C AVEC TOUS LES CISAILLEMENTS
  409. C
  410. EPSR(4)= VART(KD+4)
  411. EPSR(5)= VART(KD+5)
  412. GO TO 30
  413. C
  414. 40 CONTINUE
  415. C
  416. C COMPORTEMENT AXISYMETRIQUE
  417. C
  418. EPSR(4)= VART(KD+4)
  419. C
  420. 30 CONTINUE
  421. C
  422. C---------------------------------------------------------------------
  423. IF(IFLU.EQ.0) CALL ZERO(EPFLU,6,1)
  424. IF(ICAL.EQ.1) EPO=RB/EX
  425. EPO1=EPO
  426. IF(ICAL.NE.1.AND.ICU.EQ.3) EPO1=0.85D0*EPO
  427. * *-----------------------*
  428. * * MODELE D' OTTOSEN *
  429. * *-----------------------*
  430. *
  431. CALL PAMAOT(ALPHA,AA,BB,DK1,DK2)
  432. CALL BETDJE(EPSR,SIGR,SIGM,NSTRS,D,D1,IFOU,STRN,EPSPL,
  433. A DEPIN,EPFLU,XE,NBNN,MELE,wrk12)
  434. C---------------------------------------------------------------------
  435. NSTRS1=NSTRSS
  436. IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.GT.3) NSTRS1=3
  437. C
  438. IF((NSTRS1.EQ.4.OR.NSTRS1.EQ.6).AND.IFOU.EQ.-2) THEN
  439. SIGR(4)=SIGR(3)
  440. SIGR(3)=0.D0
  441. *
  442. DEPFZ =DEPIN(4)
  443. DEPIN(4)=DEPIN(3)
  444. DEPIN(3)=DEPFZ
  445. ENDIF
  446. C
  447. DO 2 I=1,NSTRS1
  448. SIGF(I) =SIGR(I)*EPAIST
  449. * EPINF(I)=EPIN0(I)+DEPIN(I)
  450. EPINF(I)=DEPIN(I)
  451. 2 CONTINUE
  452. *
  453. IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.6) THEN
  454. CALL ZERO(D1,6,6)
  455. CALL ZERO(D2,6,6)
  456. CALL ZERO(SIGM,6,1)
  457. *
  458. CALL DOHCOM(XMAT,NMATT,CMATE,IFOUR,NSTRSS,D1,IRTD)
  459. CALL HOOKMU(EPAIST,0.D0,NSTRSS,D1,D2)
  460. CALL BST(D2,DEPST,NSTRSS,NSTRSS,SIGM)
  461. *
  462. DO I = 4,NSTRSS
  463. SIGF(I)=SIG0(I)+SIGM(I)
  464. * EPINF(I)=EPIN0(I)+DEPIN(I)
  465. EPINF(I)=DEPIN(I)
  466. END DO
  467. ENDIF
  468. *
  469. IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.4) THEN
  470. SIGF(3) =0.D0
  471. SIGF(4) =SIGR(3)*EPAIST
  472. DEPFZ =DEPIN(4)
  473. DEPIN(4)=DEPIN(3)
  474. DEPIN(3)=DEPFZ
  475. ENDIF
  476. C---------------------------------------------------------------------
  477. IF(IECR.EQ.20) THEN
  478. WRITE(*,*) ' ** IFLU =',IFLU,' RMOY=',RMOY,' PHIF=',PHIF,
  479. & 'TEMP0= ',TEMP0,' DTEMP =',DTEMP
  480. WRITE(*,*) ' ** DEPIN *** '
  481. WRITE(*,1991) (DEPIN(IC),IC=1,NSTRSS)
  482. WRITE(*,*) ' ** EPIN0 / EPINF '
  483. WRITE(*,1991) (EPIN0(IC),IC=1,NSTRSS),(EPINF(IC),IC=1,NSTRSS)
  484. ENDIF
  485. C-----------------------------------------------------------------------
  486. C
  487. VART( 1)=EPSPL
  488. VART( 2)=TANG1
  489. VART( 3)=TANG2
  490. c+mdj
  491. VART( 4)=IPLA1
  492. VART( 5)=IFISU2
  493. c+mdj
  494. VART( 6)=DTR1
  495. VART( 7)=DTR2
  496. c+mdj
  497. VART( 8)=IFISU1
  498. VART( 9)=JFISU
  499. c+mdj
  500. VART(10)=TETA
  501. c+mdj
  502. VART(11)=JFISU2
  503. c+mdj
  504. VART(12)=EQSTR1
  505. VART(13)=EPSEQ1
  506. VART(14)=EPST1
  507. VART(15)=EPST2
  508. VART(16)=EQSTR2
  509. VART(17)=EPSEQ2
  510. VART(18)=EPSC1
  511. VART(19)=EPSC2
  512. c+mdj
  513. VART(20)=IPLA2
  514. c+mdj
  515. VART(21)=OUV1
  516. VART(22)=EDC1
  517. VART(23)=EDC2
  518. VART(24)=ETS1
  519. VART(25)=ETS2
  520. VART(26)=EDT1
  521. VART(27)=EDT2
  522. VART(28)=OUV2
  523. VART(29)=DEFR1
  524. VART(30)=DEFR2
  525. VART(31)=1.D0
  526. VART(32)=TEMP0
  527. C
  528. KD = 32
  529. IF(NSTRSS.EQ.4.AND.(IFOU.EQ.-1.OR.IFOU.EQ.0)) THEN
  530. VART(KD+1) =EQSTR3
  531. VART(KD+2) =EPSEQ3
  532. VART(KD+3) =EPST3
  533. VART(KD+4) =EPSC3
  534. VART(KD+5) =DEFR3
  535. VART(KD+6) =RTM3
  536. VART(KD+7) =EDC3
  537. VART(KD+8) =ETS3
  538. VART(KD+9) =EDT3
  539. VART(KD+10)=TANG3
  540. VART(KD+11)=OUV3
  541. VART(KD+12)=IPLA3
  542. VART(KD+13)=IFISU3
  543. VART(KD+14)=JFISU3
  544. KD = 46
  545. ENDIF
  546. C
  547. VART(KD+1)=EPSR(1)
  548. VART(KD+2)=EPSR(2)
  549. VART(KD+3)=EPSR(3)
  550. K1=KD+3
  551. C
  552. GO TO (11,11,31,41,51,61),NSTRS
  553. 11 CONTINUE
  554. KERRE=437
  555. WRITE(*,*) '!! ATTENTION DANS BETON NSTRS=',NSTRS
  556. STOP
  557. c+mdj
  558. C GO TO 1000
  559. c+mdj
  560. C
  561. 61 CONTINUE
  562. C
  563. C COMPORTEMENT TRIDIMENSIONNEL
  564. C
  565. VART(KD+4)=EPSR(4)
  566. VART(KD+5)=EPSR(5)
  567. VART(KD+6)=EPSR(6)
  568. K1=KD+6
  569. GO TO 31
  570. C
  571. 51 CONTINUE
  572. C
  573. C CONTRAINTES OU DEFORMATIONS PLANES
  574. C AVEC TOUS LES CISAILLEMENTS
  575. C
  576. VART(KD+4)=EPSR(4)
  577. VART(KD+5)=EPSR(5)
  578. K1=KD+5
  579. GO TO 31
  580. C
  581. 41 CONTINUE
  582. C
  583. C COMPORTEMENT AXISYMETRIQUE
  584. C
  585. VART(KD+4)=EPSR(4)
  586. K1=KD+4
  587. GO TO 31
  588. 31 CONTINUE
  589. C
  590. C CONTRAINTES OU DEFORMATIONS PLANES
  591. C SANS CISAILLEMENTS TRANSVERSAL
  592. C
  593. L1=0
  594. DO 66 I1 = 1,6
  595. DO 67 J1 = 1,6
  596. L1=L1+1
  597. VV(L1) = D(J1,I1)
  598. 67 CONTINUE
  599. 66 CONTINUE
  600. *
  601. L1=0
  602. DO 62 I1 = 1,NSTRS
  603. DO 63 J1 = 1,NSTRS
  604. K=K1+(I1-1)*NSTRS+J1
  605. L1=L1+1
  606. VART(K)=VV(L1)
  607. 63 CONTINUE
  608. 62 CONTINUE
  609. C
  610. DO 23 II=1,NVARI
  611. VARF(II)=VART(II)
  612. 23 CONTINUE
  613. C---------------------------------------------------------------------
  614. IF(IECR.EQ.10) THEN
  615. WRITE(*,66770) IBB,IGAU
  616. 66770 format(///'SORTIE DE BEINSA : element ',i6,2x,'point ',i3//)
  617. WRITE(*,66771) MELE,IFOU,MFR,NSTRS
  618. 66771 format('0 MELE=',i4,2x,'ifour=',i4,'mfr=',i2,'nstrs=',i2/)
  619. WRITE(*,66774) (DEPST(I),I=1,NSTRSS)
  620. 66774 format(2x,' depst '/(6(1x,1pe12.5)))
  621. WRITE(*,66772) (SIGF(I),I=1,NSTRSS)
  622. 66772 format(2x,' sigf '/(6(1x,1pe12.5)))
  623. WRITE(*,66773) (VARF(I),I=1,NVARI)
  624. 66773 format(2x,' varf '/(6(1x,1pe12.5)))
  625. ENDIF
  626. C
  627. 1991 FORMAT(18(1X,E12.5))
  628. 1992 FORMAT(18(1X,I6))
  629. C---------------------------------------------------------------------
  630. C
  631. 1000 CONTINUE
  632. C
  633. RETURN
  634. END
  635.  
  636.  
  637.  
  638.  
  639.  
  640.  
  641.  
  642.  

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