Download eligeh.eso
Back to the list
1 : C ELIGEH SOURCE CHAT 12/04/06 21:15:07 7348
2 : C
3 : SUBROUTINE ELIGEH(EPSO,G12,S1T,S2T,S3T,T1T,T3T,ALFA,EPPP,EPPM
4 : . ,EPPPF,EPPMF,XKTAN,TAUF)
5 : IMPLICIT INTEGER(I-N)
6 : IMPLICIT REAL*8(A-H,O-Z)
7 : C----------------------------------------------------------------------
8 : C
9 : C LOI ELIGEHAUSEN 1D EXPRIMEE EN CONTRAINTE - GLISSEMENT
10 : C UTILISEE SUR DES ELEMENTS JOINTS 2D
11 : C ET POUR LES LOIS D'ANCRAGE (BARRE et FIBRE)
12 : C----------------------------------------------------------------------
13 : C Didier COMBESCURE - EMSI - MAI 2000
14 : C----------------------------------------------------------------------
15 : C
16 : C EPSO: Glissement total
17 : C G12: Module elastique (charge et décharge)
18 : C EPPP: Glissement plastique positif
19 : C EPPM: Glissement plastique négatif
20 : C EPPPF: Glissement plastique positif final
21 : C EPPMF: Glissement plastique négatif final
22 : C TAUF: Contrainte de cisaillement finale
23 : C KTAN: Module tangent
24 : C
25 : C S1T, S2T et S3T: Glissements caractéristiques
26 : C T1T et T3T: Contraintes de cisaillement caractéristiques
27 : C ALFA: Paramètre définissant la première partie de la courbe
28 : C
29 : C -----------------------------------
30 : C
31 : C VERIFICATION DES PARAMETRES
32 : C
33 : C
34 :
35 : IF ((S1T.LT.0.D0).OR.(S2T.LT.0.D0).OR.(S3T.LT.0.D0).OR.
36 : . (T1T.LT.0.D0).OR.(T3T.LT.0.D0).OR.(S2T.LT.S1T).OR.
37 : . (S3T.LT.S2T).OR.(T1T.LT.T3T).OR.(ALFA.GT.1.D0).OR.
38 : . (ALFA.LT.0.D0)) THEN
39 : WRITE (*,*) 'Vos paramètres Eligehausen ne sont pas corrects'
40 : ENDIF
41 : C
42 : C
43 : C
44 : IF (EPSO.GE.EPPP) THEN
45 : C
46 : C CAS EPSO>EPPP : ELASTIQUE OU PLASTIQUE
47 : C
48 : C PREDICTION ELASTIQUE
49 : C
50 : TAUFELA = G12*(EPSO - EPPP)
51 : C
52 : C COURBE NON LINEAIRE
53 : C
54 : IF (EPSO.LT.S1T) THEN
55 : TAUFPLA = T1T*((EPSO/S1T)**ALFA)
56 : XKTANP = ALFA*(T1T/S1T)*((EPSO/S1T)**(ALFA - 1.D0))
57 : ELSEIF (EPSO.LT.S2T) THEN
58 : TAUFPLA = T1T
59 : XKTANP = 0.D0
60 : ELSEIF (EPSO.LT.S3T) THEN
61 : TAUFPLA = T1T - ((T1T - T3T)*(EPSO - S2T)/(S3T - S2T))
62 : XKTANP = (T3T - T1T)/(S3T - S2T)
63 : ELSE
64 : TAUFPLA = T3T
65 : XKTANP = 0.D0
66 : ENDIF
67 : C
68 : C
69 : C
70 : IF (TAUFELA.GT.TAUFPLA) THEN
71 : TAUF = TAUFPLA
72 : EPPPF = EPSO - (TAUFPLA/G12)
73 : EPPMF = EPPM
74 : XKTAN = XKTANP
75 : ELSE
76 : TAUF = TAUFELA
77 : EPPPF = EPPP
78 : EPPMF = EPPM
79 : XKTAN = G12
80 : ENDIF
81 :
82 : C
83 : ELSEIF (EPSO.GT.EPPM) THEN
84 : C
85 : C CAS EPPM