31
32
33
34#include "implicit_f.inc"
35
36
37
38#include "mvsiz_p.inc"
39
40
41
42 INTEGER, INTENT(IN) :: JCVT
43 INTEGER
44
46 . sig(nel,6), wxx(mvsiz), wyy(mvsiz), wzz(mvsiz)
47
48
49
50
51
52
53 INTEGER I
54
56 . s1, s2, s3, s4, s5, s6,
57 . q1, q2, q3
58
59 IF (jcvt == 0) THEN
60 DO i=1,nel
61 s1=sig(i,1)
62 s2=sig(i,2)
63 s3=sig(i,3)
64 s4=sig(i,4)
65 s5=sig(i,5)
66 s6=sig(i,6)
67
68 q1=two*s4*wzz(i)
69 q2=two*s6*wyy(i)
70 q3=two*s5*wxx(i)
71
72
73 sig(i,1)=s1-q1+q2
74 sig(i,2)=s2+q1-q3
75 sig(i,3)=s3-q2+q3
76 sig(i,4)=s4+wzz(i)*(s1-s2)+wyy(i)*s5-wxx(i)*s6
77 sig(i,5)=s5+wxx(i)*(s2-s3)+wzz(i)*s6-wyy(i)*s4
78 sig(i,6)=s6+wyy(i)*(s3-s1)+wxx(i)*s4-wzz(i)*s5
79 ENDDO
80 ENDIF
81
82 RETURN