54 xnu12, xnu13, xnu23, g12, g13, g23, vx1o, vy1o, vz1o, vx2o, vy2o, vz2o, &
60 INTEGER :: nummatvol, nummatortho
61 REAL*8,
DIMENSION(1:6,1:6,1:NumMatVol) :: ci_full
62 REAL*8,
DIMENSION(1:9,1:NumMatVol) :: ci
63 INTEGER,
DIMENSION(1:NumMatVol) :: matortho
64 REAL*8,
DIMENSION(1:NumMatOrtho) :: e11, e22, e33
65 REAL*8,
DIMENSION(1:NumMatOrtho) :: xnu12, xnu13, xnu23
66 REAL*8,
DIMENSION(1:NumMatOrtho) :: g12, g13, g23
67 REAL*8,
DIMENSION(1:NumMatOrtho) :: vx1o, vy1o, vz1o
68 REAL*8,
DIMENSION(1:NumMatOrtho) :: vx2o, vy2o, vz2o
69 REAL*8,
DIMENSION(1:NumMatOrtho) :: vx3o, vy3o, vz3o
72 INTEGER ::
i, ii,
j,
k
74 REAL*8,
DIMENSION(1:6,1:6) :: temp0, temp1, temp2, temp3, temp4
75 REAL*8,
DIMENSION(1:6,1:6) :: t, ti,
r, ri
76 REAL*8,
DIMENSION(1:3) ::
x,
y,
z, v1, v2, v3
77 REAL*8 :: m1, m2, m3, n1, n2, n3,
p1, p2, p3
81 IF ( matortho(
i) == 0 )
THEN
85 ci_full(1,1,
i) = ci(1,
i)
86 ci_full(1,2,
i) = ci(2,
i)
87 ci_full(1,3,
i) = ci(4,
i)
88 ci_full(2,1,
i) = ci(2,
i)
89 ci_full(2,2,
i) = ci(3,
i)
90 ci_full(2,3,
i) = ci(5,
i)
91 ci_full(3,1,
i) = ci(4,
i)
92 ci_full(3,2,
i) = ci(5,
i)
93 ci_full(3,3,
i) = ci(6,
i)
94 ci_full(4,4,
i) = ci(7,
i)
95 ci_full(5,5,
i) = ci(8,
i)
96 ci_full(6,6,
i) = ci(9,
i)
108 temp0(1,1) = 1.0/e11(ii)
109 temp0(1,2) = -xnu12(ii)/e11(ii)
110 temp0(1,3) = -xnu13(ii)/e11(ii)
111 temp0(2,1) = -xnu12(ii)/e11(ii)
112 temp0(2,2) = 1.0/e22(ii)
113 temp0(2,3) = -xnu23(ii)/e11(ii)
114 temp0(3,1) = -xnu13(ii)/e11(ii)
115 temp0(3,2) = -xnu23(ii)/e11(ii)
116 temp0(3,3) = 1.0/e33(ii)
117 temp0(4,4) = 1.0/g23(ii)
118 temp0(5,5) = 1.0/g13(ii)
119 temp0(6,6) = 1.0/g12(ii)
214 t(4,4) = n2*p3 + n3*p2
215 t(4,5) = p2*m3 + p3*m2
216 t(4,6) = m2*n3 + m3*n2
217 t(5,4) = n3*
p1 + n1*p3
218 t(5,5) = p3*m1 +
p1*m3
219 t(5,6) = m3*n1 + m1*n3
220 t(6,4) = n1*p2 + n2*
p1
221 t(6,5) =
p1*m2 + p2*m1
222 t(6,6) = m1*n2 + m2*n1
228 temp2 = matmul(ti,temp1)
229 temp3 = matmul(temp2,
r)
230 temp4 = matmul(temp3, t)
231 ci_full(:,:,
i) = matmul(temp4, ri)
void int int REAL REAL * y
subroutine invert6x6(T, Ti)
void int int int REAL REAL REAL * z
subroutine vol_elem_mat_ortho(ci_full, ci, NumMatVol, NumMatOrtho, MatOrtho, E11, E22, E33, xnu12, xnu13, xnu23, G12, G13, G23, vx1o, vy1o, vz1o, vx2o, vy2o, vz2o, vx3o, vy3o, vz3o)
long double dot_product(pnt vec1, pnt vec2)