1 : C NORMAL SOURCE CHAT 05/01/13 02:00:11 5004
2 : SUBROUTINE NORMAL(A,B,C,SIGMA,ANOR,ALPHA,SELAS,IBOU,ITYP)
3 : C
4 : C CALCUL DE LA NORMNALE A LA SURFACE DE PLASTICITE
5 : C
6 : IMPLICIT INTEGER(I-N)
7 : IMPLICIT REAL*8(A-H,O-Z)
8 : DIMENSION SIGMA(*),ANOR(*)
9 : C CALCUL DE M*SIGMA
10 : DO 1001 IA=1,IBOU
11 : 1001 ANOR(IA)=0.
12 : GO TO (1,2,2,4,5,6,7,8,9,10,99,99,99,14),ITYP
13 : 1 CONTINUE
14 : ANOR(1)=A*SIGMA(1)-B*(SIGMA(2)+SIGMA(3))
15 : ANOR(2)=A*SIGMA(2)-B*(SIGMA(3)+SIGMA(1))
16 : ANOR(3)=A*SIGMA(3)-B*(SIGMA(1)+SIGMA(2))
17 : DO 1002 IA=1,3
18 : 1002 ANOR(IA+3)=C*SIGMA(3+IA)
19 : GO TO 1000
20 : 2 CONTINUE
21 : ANOR(1)=ALPHA*(SIGMA(1)*A-B*SIGMA(2))
22 : ANOR(2)=ALPHA*(SIGMA(2)*A-B*SIGMA(1))
23 : ANOR(4)= SIGMA(4)*A-B*SIGMA(5)
24 : ANOR(5)= SIGMA(5)*A-B*SIGMA(4)
25 : ANOR(3)=C*SIGMA(3)*ALPHA
26 : ANOR(6)=C*SIGMA(6)
27 : GO TO 1000
28 : 4 ANOR(3)=SIGMA(3)
29 : GO TO 1000
30 : 5 CONTINUE
31 : GO TO 1000
32 : 6 ANOR(1)=A*SIGMA(1)-B*SIGMA(2)
33 : ANOR(2)=A*SIGMA(2)-B*SIGMA(1)
34 : ANOR(4)=C*SIGMA(4)
35 : GO TO 1000
36 : 7 CONTINUE
37 : ANOR(1)=ALPHA*A*SIGMA(1)
38 : ANOR(4)= A*SIGMA(4)
39 : GO TO 1000
40 : 8 ANOR(1)=SIGMA(1)
41 : GO TO 1000
42 : C= Modes de calcul UNIDIMENSIONNELs (1D)
43 : 14 ANOR(1)=A*SIGMA(1)-B*(SIGMA(2)+SIGMA(3))
44 : ANOR(2)=A*SIGMA(2)-B*(SIGMA(3)+SIGMA(1))
45 : ANOR(3)=A*SIGMA(3)-B*(SIGMA(1)+SIGMA(2))
46 : C* GOTO 1000
47 : 9 CONTINUE
48 : 10 CONTINUE
49 : C NORMALISATION
50 : 1000 CONTINUE
51 : X=1./SELAS
52 : Y=0.
53 : DO 1003 IA=1,IBOU
54 : ANOR(IA)=ANOR(IA)*X
55 : 1003 Y=Y+ANOR(IA)*ANOR(IA)
56 : Y=1./SQRT(Y)
57 : DO 1004 IA=1,IBOU
58 : 1004 ANOR(IA)=ANOR(IA)*Y
59 : RETURN
60 : 99 WRITE(6,7999)
61 : 7999 FORMAT('0 NORMAL - CAS NON IMPLEMENTE ' /)
62 : RETURN
63 : END
64 :
65 :
66 :
67 :
© Cast3M 2003 - All rights reserved.
Disclaimer