38
39
40
41#include "implicit_f.inc"
42
43
44
45#include "mvsiz_p.inc"
46
47
48
49
50
51
52 INTEGER NEL
53
55 . px1(*), px2(*), px3(*), px4(*),
56 . py1(*), py2(*), py3(*), py4(*),
57 . pz1(*), pz2(*), pz3(*), pz4(*),
58 . k11(3,3,*),k12(3,3,*),k13(3,3,*),k14(3,3,*) ,
59 . k22(3,3,*) ,k23(3,3,*) ,k24(3,3,*),
60 . k33(3,3,*),k34(3,3,*),
61 . k44(3,3,*),vol(*),sig(nel,6)
62
63
64
65 INTEGER I,J
66
68 . fxx(mvsiz),fyy(mvsiz),fzz(mvsiz),fxy(mvsiz),fxz(mvsiz),
69 . fyz(mvsiz)
70
71 DO i=1,nel
72 fxx(i)=sig(i,1)*vol(i)
73 fyy(i)=sig(i,2)*vol(i)
74 fzz(i)=sig(i,3)*vol(i)
75 fxy(i)=sig(i,4)*vol(i)
76 fyz(i)=sig(i,5)*vol(i)
77 fxz(i)=sig(i,6)*vol(i)
78 ENDDO
79
81 . px1, py1 ,pz1 ,px1, py1, pz1,
82 . fxx, fyy ,fzz ,fxy, fyz, fxz,k11 )
84 . px2, py2 ,pz2 ,px2, py2, pz2,
85 . fxx, fyy ,fzz ,fxy, fyz, fxz,k22 )
87 . px3, py3 ,pz3 ,px3, py3, pz3,
88 . fxx, fyy ,fzz ,fxy, fyz, fxz,k33 )
90 . px4, py4 ,pz4 ,px4, py4, pz4,
91 . fxx, fyy ,fzz ,fxy, fyz, fxz,k44 )
93 . px1, py1 ,pz1 ,px2, py2, pz2,
94 . fxx, fyy ,fzz ,fxy, fyz, fxz,k12 )
96 . px1, py1 ,pz1 ,px3, py3, pz3,
97 . fxx, fyy ,fzz ,fxy, fyz, fxz,k13 )
99 . px1, py1 ,pz1 ,px4, py4, pz4,
100 . fxx, fyy ,fzz ,fxy, fyz, fxz,k14 )
102 . px2, py2 ,pz2 ,px3, py3, pz3,
103 . fxx, fyy ,fzz ,fxy, fyz, fxz,k23 )
105 . px2, py2 ,pz2 ,px4, py4, pz4,
106 . fxx, fyy ,fzz ,fxy, fyz, fxz,k24 )
108 . px3, py3 ,pz3 ,px4, py4, pz4,
109 . fxx, fyy ,fzz ,fxy, fyz, fxz,k34 )
110
111 RETURN
subroutine s4keh3(lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, fxx, fyy, fzz, fxy, fyz, fxz, kij)