46
47
48
50
51
52
53#include "implicit_f.inc"
54
55
56
57#include "mvsiz_p.inc"
58
59
60
61 INTEGER, INTENT(IN) :: NEL
63 . voln(*),deltax(*),
64 . x1(*), x2(*), x3(*), x4(*), x5(*), x6(*), x7(*), x8(*),
65 . y1(*), y2(*), y3(*), y4(*), y5(*), y6(*), y7(*), y8(*),
66 . z1(*), z2(*), z3(*), z4(*), z5(*), z6(*), z7(*), z8(*),
67 . n1x(*), n2x(*), n3x(*), n4x(*), n5x(*), n6x(*),
68 . n1y(*), n2y(*), n3y(*), n4y(*), n5y(*), n6y(*),
69 . n1z(*), n2z(*), n3z(*), n4z(*), n5z(*), n6z(*)
70
71
72
73 INTEGER I, NFAC
74 my_real xioff(mvsiz), aream(mvsiz), atest(mvsiz),
area(6,mvsiz)
75
76 DO i=1,nel
77 xioff(i)=one
78 aream(i)=em20
79 ENDDO
80
81 IF(
ale%GLOBAL%ICAA==0)
THEN
83 1 x1, x2, x3, x4,
84 2 y1, y2, y3, y4,
85 3 z1, z2, z3, z4,
86 4 1,
area, aream, nel)
88 1 x5, x6, x7, x8,
89 2 y5, y6, y7, y8,
90 3 z5, z6, z7, z8,
91 4 2,
area, aream, nel)
93 1 x1, x2, x6, x5,
94 2 y1, y2, y6, y5,
95 3 z1, z2, z6, z5,
96 4 3,
area, aream, nel)
98 1 x2, x3, x7, x6,
99 2 y2, y3, y7, y6,
100 3 z2, z3, z7, z6,
101 4 4,
area, aream, nel)
103 1 x3, x4, x8, x7,
104 2 y3, y4, y8, y7,
105 3 z3, z4, z8, z7,
106 4 5,
area, aream, nel)
108 1 x4, x1, x5, x8,
109 2 y4, y1, y5, y8,
110 3 z4, z1, z5, z8,
111 4 6,
area, aream, nel)
112
113 DO i=1,nel
114 deltax(i)=four*voln(i)*xioff(i)/sqrt(aream(i))
115 ENDDO
116
117 ELSE
119 1 n1x, n2x, n3x, n4x,
120 2 n5x, n6x, n1y, n2y,
121 3 n3y, n4y, n5y, n6y,
122 4 n1z, n2z, n3z, n4z,
123 5 n5z, n6z, x1, x2,
124 6 x3, x4, x5, x6,
125 7 x7, x8, y1, y2,
126 8 y3, y4, y5, y6,
127 9 y7, y8, z1, z2,
128 a z3, z4, z5, z6,
129 b z7, z8, aream, nel)
130 DO i=1,nel
131 deltax(i)=two*voln(i)*xioff(i)/sqrt(aream(i))
132 ENDDO
133 ENDIF
134 RETURN
subroutine area(d1, x, x2, y, y2, eint, stif0)
subroutine slena(n1x, n2x, n3x, n4x, n5x, n6x, n1y, n2y, n3y, n4y, n5y, n6y, n1z, n2z, n3z, n4z, n5z, n6z, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, aream, nel)
subroutine slen(x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, j, area, aream)
1.15.0