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
268 : CALL ZERO(D,6,6)
268 : CALL ZERO(D1,6,6)
268 : CALL ZERO(EPSR,6,1)
268 : CALL ZERO(DEPIN,6,1)
268 : CALL ZERO(EPFLU,6,1)
268 : CALL ZERO(VART,100,1)
268 : CALL ZERO(SIGR,6,1)
268 : CALL ZERO(SIGM,6,1)
268 : CALL ZERO(STRN,6,1)
269 : C
270 : EPAIST=1.D0
271 : NSTRS=NSTRSS
272 : *
273 : IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.6) THEN
274 : NSTRS=3
275 : EPAIST=XCAR(1)
276 : ENDIF
277 : *
278 : DO 1 I=1,NSTRS
279 : SIGR(I)=SIG0(I)/EPAIST
280 : STRN(I)=DEPST(I)
281 : EPFLU(I)=EPIN0(I)
282 : 1 CONTINUE
283 : *
284 : IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.4) THEN
285 : NSTRS=3
286 : SIGR(3) =SIG0(4)
287 : STRN(3) =DEPST(4)
288 : ENDIF
289 : C---------------------------------------------------------------------
290 : IF(IECR.EQ.10) THEN
291 : WRITE(*,*) ' ** SIGR='
292 : WRITE(*,1991) (SIGR(I),I=1,NSTRS)
293 : WRITE(*,*) ' ** STRN='
294 : WRITE(*,1991) (STRN(I),I=1,NSTRS)
295 : ENDIF
296 : C------------------------------------------------------
297 : C INITIALISATION DES VARIABLES INTERNES
298 : C------------------------------------------------------
299 : DO II=1,NVARI
300 : VART(II)=VAR0(II)
301 : END DO
302 : C
303 : EPSPL = VART( 1)
304 : TANG1 = VART( 2)
305 : TANG2 = VART( 3)
306 : IPLA1 = INT(REAL(VART( 4)))
307 : IFISU2 = INT(REAL(VART( 5)))
308 : DTR1 = VART( 6)
309 : DTR2 = VART( 7)
310 : IFISU1 = INT(REAL(VART( 8)))
311 : JFISU = INT(REAL(VART( 9)))
312 : ANGL = VART(10)
313 : JFISU2 = INT(REAL(VART(11)))
314 : EQSTR1 = VART(12)
315 : EPSEQ1 = VART(13)
316 : EPST1 = VART(14)
317 : EPST2 = VART(15)
318 : EQSTR2 = VART(16)
319 : EPSEQ2 = VART(17)
320 : EPSC1 = VART(18)
321 : EPSC2 = VART(19)
322 : IPLA2 = INT(REAL(VART(20)))
323 : OUV1 = VART(21)
324 : EDC1 = VART(22)
325 : EDC2 = VART(23)
326 : ETS1 = VART(24)
327 : ETS2 = VART(25)
328 : EDT1 = VART(26)
329 : EDT2 = VART(27)
330 : OUV2 = VART(28)
331 : DEFR1 = VART(29)
332 : DEFR2 = VART(30)
333 : IPAS0 = INT(REAL(VART(31)))
334 : TEMP0 = VART(32)
335 : C
336 : KD = 32
337 : IF(NSTRS.EQ.4.AND.(IFOU.EQ.-1.OR.IFOU.EQ.0)) THEN
338 : EQSTR3 = VART(KD+1)
339 : EPSEQ3 = VART(KD+2)
340 : EPST3 = VART(KD+3)
341 : EPSC3 = VART(KD+4)
342 : DEFR3 = VART(KD+5)
343 : RTM3 = VART(KD+6)
344 : EDC3 = VART(KD+7)
345 : ETS3 = VART(KD+8)
346 : EDT3 = VART(KD+9)
347 : TANG3 = VART(KD+10)
348 : OUV3 = VART(KD+11)
349 : IPLA3 = INT(REAL(VART(KD+12)))
350 : IFISU3 = INT(REAL(VART(KD+13)))
351 : JFISU3 = INT(REAL(VART(KD+14)))
352 : KD = 46
353 : ENDIF
354 : C
355 : EPSR(1)= VART(KD+1)
356 : EPSR(2)= VART(KD+2)
357 : EPSR(3)= VART(KD+3)
358 : C -------------------------------------
359 : IF(IFLU.NE.0.AND.(TEMP0+DTEMP-TEMP1).LE.0.D0) THEN
360 : WRITE(*,*)'!! ATTENTION DANS LE MODELE DE FLUAGE BETON '
361 : TEMPF = TEMP0+DTEMP-TEMP1
362 : WRITE(*,*)' (LOI BPEL91) TEMPS NEGATIF TEMPF=',TEMPF
363 : WRITE(*,*)' TEMP0=',TEMP0,' TEMP1=',TEMP1,' DTEMP=',DTEMP
364 : STOP
365 : ENDIF
366 : C -------------------------------------
367 : TETA=ANGL
368 : IF((NSTRS.EQ.4.OR.NSTRS.EQ.6).AND.IFOU.EQ.-2) THEN
369 : NSTRS=3
370 : SIGR(3)=SIG0(4)
371 : STRN(3)=DEPST(4)
372 : STRN(4)=DEPST(3)
373 : SIGR(4)=0.D0
374 : EPFLU(3)=EPIN0(4)
375 : EPFLU(4)=EPIN0(3)
376 : ENDIF
377 : C
378 : GO TO (10,10,30,40,50,60),NSTRS
379 : 10 CONTINUE
380 : KERRE=437
381 : WRITE(*,*) '!! ATTENTION DANS BETON NSTRS=',NSTRS
382 : STOP
383 : c+mdj
384 : C GO TO 1000
385 : c+mdj
386 : C
387 : 60 CONTINUE
388 : C
389 : C COMPORTEMENT TRIDIMENSIONNEL
390 : C
391 : EPSR(4)= VART(KD+4)
392 : EPSR(5)= VART(KD+5)
393 : EPSR(6)= VART(KD+6)
394 : GO TO 30
395 : C
396 : 50 CONTINUE
397 : C
398 : C CONTRAINTES OU DEFORMATIONS PLANES
399 : C AVEC TOUS LES CISAILLEMENTS
400 : C
401 : EPSR(4)= VART(KD+4)
402 : EPSR(5)= VART(KD+5)
403 : GO TO 30
404 : C
405 : 40 CONTINUE
406 : C
407 : C COMPORTEMENT AXISYMETRIQUE
408 : C
409 : EPSR(4)= VART(KD+4)
410 : C
411 : 30 CONTINUE
412 : C
413 : C---------------------------------------------------------------------
413 : IF(IFLU.EQ.0) CALL ZERO(EPFLU,6,1)
414 : IF(ICAL.EQ.1) EPO=RB/EX
415 : EPO1=EPO
416 : IF(ICAL.NE.1.AND.ICU.EQ.3) EPO1=0.85D0*EPO
417 : * *-----------------------*
418 : * * MODELE D' OTTOSEN *
419 : * *-----------------------*
420 : *
420 : CALL PAMAOT(ALPHA,AA,BB,DK1,DK2)
420 : CALL BETDJE(EPSR,SIGR,SIGM,NSTRS,D,D1,IFOU,STRN,EPSPL,
421 : A DEPIN,EPFLU,XE,NBNN,MELE,wrk12)
422 : C---------------------------------------------------------------------
423 : NSTRS1=NSTRSS
424 : IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.GT.3) NSTRS1=3
425 : C
426 : IF((NSTRS1.EQ.4.OR.NSTRS1.EQ.6).AND.IFOU.EQ.-2) THEN
427 : SIGR(4)=SIGR(3)
428 : SIGR(3)=0.D0
429 : *
430 : DEPFZ =DEPIN(4)
431 : DEPIN(4)=DEPIN(3)
432 : DEPIN(3)=DEPFZ
433 : ENDIF
434 : C
435 : DO 2 I=1,NSTRS1
436 : SIGF(I) =SIGR(I)*EPAIST
437 : * EPINF(I)=EPIN0(I)+DEPIN(I)
438 : EPINF(I)=DEPIN(I)
439 : 2 CONTINUE
440 : *
441 : IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.6) THEN
441 : CALL ZERO(D1,6,6)
441 : CALL ZERO(D2,6,6)
441 : CALL ZERO(SIGM,6,1)
442 : *
442 : CALL DOHCOM(XMAT,NMATT,CMATE,IFOUR,NSTRSS,D1,IRTD)
442 : CALL HOOKMU(EPAIST,0.D0,NSTRSS,D1,D2)
442 : CALL BST(D2,DEPST,NSTRSS,NSTRSS,SIGM)
443 : *
444 : DO I = 4,NSTRSS
445 : SIGF(I)=SIG0(I)+SIGM(I)
446 : * EPINF(I)=EPIN0(I)+DEPIN(I)
447 : EPINF(I)=DEPIN(I)
448 : END DO
449 : ENDIF
450 : *
451 : IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.4) THEN
452 : SIGF(3) =0.D0
453 : SIGF(4) =SIGR(3)*EPAIST
454 : DEPFZ =DEPIN(4)
455 : DEPIN(4)=DEPIN(3)
456 : DEPIN(3)=DEPFZ
457 : ENDIF
458 : C---------------------------------------------------------------------
459 : IF(IECR.EQ.20) THEN
460 : WRITE(*,*) ' ** IFLU =',IFLU,' RMOY=',RMOY,' PHIF=',PHIF,
461 : & 'TEMP0= ',TEMP0,' DTEMP =',DTEMP
462 : WRITE(*,*) ' ** DEPIN *** '
463 : WRITE(*,1991) (DEPIN(IC),IC=1,NSTRSS)
464 : WRITE(*,*) ' ** EPIN0 / EPINF '
465 : WRITE(*,1991) (EPIN0(IC),IC=1,NSTRSS),(EPINF(IC),IC=1,NSTRSS)
466 : ENDIF
467 : C-----------------------------------------------------------------------
468 : C
469 : VART( 1)=EPSPL
470 : VART( 2)=TANG1
471 : VART( 3)=TANG2
472 : c+mdj
473 : VART( 4)=IPLA1
474 : VART( 5)=IFISU2
475 : c+mdj
476 : VART( 6)=DTR1
477 : VART( 7)=DTR2
478 : c+mdj
479 : VART( 8)=IFISU1
480 : VART( 9)=JFISU
481 : c+mdj
482 : VART(10)=TETA
483 : c+mdj
484 : VART(11)=JFISU2
485 : c+mdj
486 : VART(12)=EQSTR1
487 : VART(13)=EPSEQ1
488 : VART(14)=EPST1
489 : VART(15)=EPST2
490 : VART(16)=EQSTR2
491 : VART(17)=EPSEQ2
492 : VART(18)=EPSC1
493 : VART(19)=EPSC2
494 : c+mdj
495 : VART(20)=IPLA2
496 : c+mdj
497 : VART(21)=OUV1
498 : VART(22)=EDC1
499 : VART(23)=EDC2
500 : VART(24)=ETS1
501 : VART(25)=ETS2
502 : VART(26)=EDT1
503 : VART(27)=EDT2
504 : VART(28)=OUV2
505 : VART(29)=DEFR1
506 : VART(30)=DEFR2
507 : VART(31)=1.D0
508 : VART(32)=TEMP0
509 : C
510 : KD = 32
511 : IF(NSTRSS.EQ.4.AND.(IFOU.EQ.-1.OR.IFOU.EQ.0)) THEN
512 : VART(KD+1) =EQSTR3
513 : VART(KD+2) =EPSEQ3
514 : VART(KD+3) =EPST3
515 : VART(KD+4) =EPSC3
516 : VART(KD+5) =DEFR3
517 : VART(KD+6) =RTM3
518 : VART(KD+7) =EDC3
519 : VART(KD+8) =ETS3
520 : VART(KD+9) =EDT3
521 : VART(KD+10)=TANG3
522 : VART(KD+11)=OUV3
523 : VART(KD+12)=IPLA3
524 : VART(KD+13)=IFISU3
525 : VART(KD+14)=JFISU3
526 : KD = 46
527 : ENDIF
528 : C
529 : VART(KD+1)=EPSR(1)
530 : VART(KD+2)=EPSR(2)
531 : VART(KD+3)=EPSR(3)
532 : K1=KD+3
533 : C
534 : GO TO (11,11,31,41,51,61),NSTRS
535 : 11 CONTINUE
536 : KERRE=437
537 : WRITE(*,*) '!! ATTENTION DANS BETON NSTRS=',NSTRS
538 : STOP
539 : c+mdj
540 : C GO TO 1000
541 : c+mdj
542 : C
543 : 61 CONTINUE
544 : C
545 : C COMPORTEMENT TRIDIMENSIONNEL
546 : C
547 : VART(KD+4)=EPSR(4)
548 : VART(KD+5)=EPSR(5)
549 : VART(KD+6)=EPSR(6)
550 : K1=KD+6
551 : GO TO 31
552 : C
553 : 51 CONTINUE
554 : C
555 : C CONTRAINTES OU DEFORMATIONS PLANES
556 : C AVEC TOUS LES CISAILLEMENTS
557 : C
558 : VART(KD+4)=EPSR(4)
559 : VART(KD+5)=EPSR(5)
560 : K1=KD+5
561 : GO TO 31
562 : C
563 : 41 CONTINUE
564 : C
565 : C COMPORTEMENT AXISYMETRIQUE
566 : C
567 : VART(KD+4)=EPSR(4)
568 : K1=KD+4
569 : GO TO 31
570 : 31 CONTINUE
571 : C
572 : C CONTRAINTES OU DEFORMATIONS PLANES
573 : C SANS CISAILLEMENTS TRANSVERSAL
574 : C
575 : L1=0
576 : DO 66 I1 = 1,6
577 : DO 67 J1 = 1,6
578 : L1=L1+1
579 : VV(L1) = D(J1,I1)
580 : 67 CONTINUE
581 : 66 CONTINUE
582 : *
583 : L1=0
584 : DO 62 I1 = 1,NSTRS
585 : DO 63 J1 = 1,NSTRS
586 : K=K1+(I1-1)*NSTRS+J1
587 : L1=L1+1
588 : VART(K)=VV(L1)
589 : 63 CONTINUE
590 : 62 CONTINUE
591 : C
592 : DO 23 II=1,NVARI
593 : VARF(II)=VART(II)
594 : 23 CONTINUE
595 : C---------------------------------------------------------------------
596 : IF(IECR.EQ.10) THEN
597 : WRITE(*,66770) IBB,IGAU
598 : 66770 format(///'SORTIE DE BEINSA : element ',i6,2x,'point ',i3//)
599 : WRITE(*,66771) MELE,IFOU,MFR,NSTRS
600 : 66771 format('0 MELE=',i4,2x,'ifour=',i4,'mfr=',i2,'nstrs=',i2/)
601 : WRITE(*,66774) (DEPST(I),I=1,NSTRSS)
602 : 66774 format(2x,' depst '/(6(1x,1pe12.5)))
603 : WRITE(*,66772) (SIGF(I),I=1,NSTRSS)
604 : 66772 format(2x,' sigf '/(6(1x,1pe12.5)))
605 : WRITE(*,66773) (VARF(I),I=1,NVARI)
606 : 66773 format(2x,' varf '/(6(1x,1pe12.5)))
607 : ENDIF
608 : C
609 : 1991 FORMAT(18(1X,E12.5))
610 : 1992 FORMAT(18(1X,I6))
611 : C---------------------------------------------------------------------
612 : C
613 : 1000 CONTINUE
614 : C
615 : RETURN
616 : END
617 :
618 :
619 :
620 :
621 :
622 :
623 :
© Cast3M 2003 - All rights reserved.
Disclaimer