Lines Matching refs:gderiv
51 real(kind=kreal) :: gderiv(nn+3, 3), stress(6) local
91 call getGlobalDeriv( etype, nn, naturalcoord, elem, det, gderiv(1:nn, 1:3) )
119 gderiv(nn+1, :) = -2.0D0*naturalcoord(1)*inverse(1, :)/det
120 gderiv(nn+2, :) = -2.0D0*naturalcoord(2)*inverse(2, :)/det
121 gderiv(nn+3, :) = -2.0D0*naturalcoord(3)*inverse(3, :)/det
126 B(1, 3*j-2) = gderiv(j, 1)
127 B(2, 3*j-1) = gderiv(j, 2)
128 B(3, 3*j ) = gderiv(j, 3)
129 B(4, 3*j-2) = gderiv(j, 2)
130 B(4, 3*j-1) = gderiv(j, 1)
131 B(5, 3*j-1) = gderiv(j, 3)
132 B(5, 3*j ) = gderiv(j, 2)
133 B(6, 3*j-2) = gderiv(j, 3)
134 B(6, 3*j ) = gderiv(j, 1)
140 gdispderiv(1:ndof,1:ndof) = matmul( unode(1:ndof, 1:nn+3), gderiv(1:nn+3, 1:ndof) )
143 B1(1, 3*j-2) = gdispderiv(1, 1)*gderiv(j, 1)
144 B1(1, 3*j-1) = gdispderiv(2, 1)*gderiv(j, 1)
145 B1(1, 3*j ) = gdispderiv(3, 1)*gderiv(j, 1)
146 B1(2, 3*j-2) = gdispderiv(1, 2)*gderiv(j, 2)
147 B1(2, 3*j-1) = gdispderiv(2, 2)*gderiv(j, 2)
148 B1(2, 3*j ) = gdispderiv(3, 2)*gderiv(j, 2)
149 B1(3, 3*j-2) = gdispderiv(1, 3)*gderiv(j, 3)
150 B1(3, 3*j-1) = gdispderiv(2, 3)*gderiv(j, 3)
151 B1(3, 3*j ) = gdispderiv(3, 3)*gderiv(j, 3)
152 B1(4, 3*j-2) = gdispderiv(1, 2)*gderiv(j, 1)+gdispderiv(1, 1)*gderiv(j, 2)
153 B1(4, 3*j-1) = gdispderiv(2, 2)*gderiv(j, 1)+gdispderiv(2, 1)*gderiv(j, 2)
154 B1(4, 3*j ) = gdispderiv(3, 2)*gderiv(j, 1)+gdispderiv(3, 1)*gderiv(j, 2)
155 B1(5, 3*j-2) = gdispderiv(1, 2)*gderiv(j, 3)+gdispderiv(1, 3)*gderiv(j, 2)
156 B1(5, 3*j-1) = gdispderiv(2, 2)*gderiv(j, 3)+gdispderiv(2, 3)*gderiv(j, 2)
157 B1(5, 3*j ) = gdispderiv(3, 2)*gderiv(j, 3)+gdispderiv(3, 3)*gderiv(j, 2)
158 B1(6, 3*j-2) = gdispderiv(1, 3)*gderiv(j, 1)+gdispderiv(1, 1)*gderiv(j, 3)
159 B1(6, 3*j-1) = gdispderiv(2, 3)*gderiv(j, 1)+gdispderiv(2, 1)*gderiv(j, 3)
160 B1(6, 3*j ) = gdispderiv(3, 3)*gderiv(j, 1)+gdispderiv(3, 1)*gderiv(j, 3)
179 BN(1, 3*j-2) = gderiv(j, 1)
180 BN(2, 3*j-1) = gderiv(j, 1)
181 BN(3, 3*j ) = gderiv(j, 1)
182 BN(4, 3*j-2) = gderiv(j, 2)
183 BN(5, 3*j-1) = gderiv(j, 2)
184 BN(6, 3*j ) = gderiv(j, 2)
185 BN(7, 3*j-2) = gderiv(j, 3)
186 BN(8, 3*j-1) = gderiv(j, 3)
187 BN(9, 3*j ) = gderiv(j, 3)
257 …real(kind=kreal) :: gderiv(nn+3,3), gderiv0(nn+3,3), gdispderiv(3,3), F(3,3), det, det0, WG, ttc… local
294 call getGlobalDeriv( etype, nn, naturalcoord, elem, det, gderiv(1:nn, 1:3) )
322 gderiv(nn+1, :) = -2.0D0*naturalcoord(1)*inverse(1, :)/det
323 gderiv(nn+2, :) = -2.0D0*naturalcoord(2)*inverse(2, :)/det
324 gderiv(nn+3, :) = -2.0D0*naturalcoord(3)*inverse(3, :)/det
329 B(1, 3*j-2) = gderiv(j, 1)
330 B(2, 3*j-1) = gderiv(j, 2)
331 B(3, 3*j ) = gderiv(j, 3)
332 B(4, 3*j-2) = gderiv(j, 2)
333 B(4, 3*j-1) = gderiv(j, 1)
334 B(5, 3*j-1) = gderiv(j, 3)
335 B(5, 3*j ) = gderiv(j, 2)
336 B(6, 3*j-2) = gderiv(j, 3)
337 B(6, 3*j ) = gderiv(j, 1)
343 gdispderiv(1:ndof,1:ndof) = matmul( totaldisp(1:ndof, 1:nn+3), gderiv(1:nn+3, 1:ndof) )
346 B1(1, 3*j-2) = gdispderiv(1, 1)*gderiv(j, 1)
347 B1(1, 3*j-1) = gdispderiv(2, 1)*gderiv(j, 1)
348 B1(1, 3*j ) = gdispderiv(3, 1)*gderiv(j, 1)
349 B1(2, 3*j-2) = gdispderiv(1, 2)*gderiv(j, 2)
350 B1(2, 3*j-1) = gdispderiv(2, 2)*gderiv(j, 2)
351 B1(2, 3*j ) = gdispderiv(3, 2)*gderiv(j, 2)
352 B1(3, 3*j-2) = gdispderiv(1, 3)*gderiv(j, 3)
353 B1(3, 3*j-1) = gdispderiv(2, 3)*gderiv(j, 3)
354 B1(3, 3*j ) = gdispderiv(3, 3)*gderiv(j, 3)
355 B1(4, 3*j-2) = gdispderiv(1, 2)*gderiv(j, 1)+gdispderiv(1, 1)*gderiv(j, 2)
356 B1(4, 3*j-1) = gdispderiv(2, 2)*gderiv(j, 1)+gdispderiv(2, 1)*gderiv(j, 2)
357 B1(4, 3*j ) = gdispderiv(3, 2)*gderiv(j, 1)+gdispderiv(3, 1)*gderiv(j, 2)
358 B1(5, 3*j-2) = gdispderiv(1, 2)*gderiv(j, 3)+gdispderiv(1, 3)*gderiv(j, 2)
359 B1(5, 3*j-1) = gdispderiv(2, 2)*gderiv(j, 3)+gdispderiv(2, 3)*gderiv(j, 2)
360 B1(5, 3*j ) = gdispderiv(3, 2)*gderiv(j, 3)+gdispderiv(3, 3)*gderiv(j, 2)
361 B1(6, 3*j-2) = gdispderiv(1, 3)*gderiv(j, 1)+gdispderiv(1, 1)*gderiv(j, 3)
362 B1(6, 3*j-1) = gdispderiv(2, 3)*gderiv(j, 1)+gdispderiv(2, 1)*gderiv(j, 3)
363 B1(6, 3*j ) = gdispderiv(3, 3)*gderiv(j, 1)+gdispderiv(3, 1)*gderiv(j, 3)
385 BN(1, 3*j-2) = gderiv(j, 1)
386 BN(2, 3*j-1) = gderiv(j, 1)
387 BN(3, 3*j ) = gderiv(j, 1)
388 BN(4, 3*j-2) = gderiv(j, 2)
389 BN(5, 3*j-1) = gderiv(j, 2)
390 BN(6, 3*j ) = gderiv(j, 2)
391 BN(7, 3*j-2) = gderiv(j, 3)
392 BN(8, 3*j-1) = gderiv(j, 3)
393 BN(9, 3*j ) = gderiv(j, 3)
480 call getGlobalDeriv(etype, nn, naturalcoord, elem, det, gderiv)
509 gderiv(nn+1, :) = -2.0D0*naturalcoord(1)*inverse(1, :)/det
510 gderiv(nn+2, :) = -2.0D0*naturalcoord(2)*inverse(2, :)/det
511 gderiv(nn+3, :) = -2.0D0*naturalcoord(3)*inverse(3, :)/det
514 gdispderiv(1:ndof, 1:ndof) = matmul( totaldisp(1:ndof, 1:nn+3), gderiv(1:nn+3, 1:ndof) )
572 B(1,3*j-2) = gderiv(j, 1)
573 B(2,3*j-1) = gderiv(j, 2)
574 B(3,3*j ) = gderiv(j, 3)
575 B(4,3*j-2) = gderiv(j, 2)
576 B(4,3*j-1) = gderiv(j, 1)
577 B(5,3*j-1) = gderiv(j, 3)
578 B(5,3*j ) = gderiv(j, 2)
579 B(6,3*j-2) = gderiv(j, 3)
580 B(6,3*j ) = gderiv(j, 1)
588 gdispderiv(1:ndof, 1:ndof) = matmul( totaldisp(1:ndof, 1:nn+3), gderiv(1:nn+3, 1:ndof) )
591 B1(1, 3*j-2) = gdispderiv(1, 1)*gderiv(j, 1)
592 B1(1, 3*j-1) = gdispderiv(2, 1)*gderiv(j, 1)
593 B1(1, 3*j ) = gdispderiv(3, 1)*gderiv(j, 1)
594 B1(2, 3*j-2) = gdispderiv(1, 2)*gderiv(j, 2)
595 B1(2, 3*j-1) = gdispderiv(2, 2)*gderiv(j, 2)
596 B1(2, 3*j ) = gdispderiv(3, 2)*gderiv(j, 2)
597 B1(3, 3*j-2) = gdispderiv(1, 3)*gderiv(j, 3)
598 B1(3, 3*j-1) = gdispderiv(2, 3)*gderiv(j, 3)
599 B1(3, 3*j ) = gdispderiv(3, 3)*gderiv(j, 3)
600 B1(4, 3*j-2) = gdispderiv(1, 2)*gderiv(j, 1)+gdispderiv(1, 1)*gderiv(j, 2)
601 B1(4, 3*j-1) = gdispderiv(2, 2)*gderiv(j, 1)+gdispderiv(2, 1)*gderiv(j, 2)
602 B1(4, 3*j ) = gdispderiv(3, 2)*gderiv(j, 1)+gdispderiv(3, 1)*gderiv(j, 2)
603 B1(5, 3*j-2) = gdispderiv(1, 2)*gderiv(j, 3)+gdispderiv(1, 3)*gderiv(j, 2)
604 B1(5, 3*j-1) = gdispderiv(2, 2)*gderiv(j, 3)+gdispderiv(2, 3)*gderiv(j, 2)
605 B1(5, 3*j ) = gdispderiv(3, 2)*gderiv(j, 3)+gdispderiv(3, 3)*gderiv(j, 2)
606 B1(6, 3*j-2) = gdispderiv(1, 3)*gderiv(j, 1)+gdispderiv(1, 1)*gderiv(j, 3)
607 B1(6, 3*j-1) = gdispderiv(2, 3)*gderiv(j, 1)+gdispderiv(2, 1)*gderiv(j, 3)
608 B1(6, 3*j ) = gdispderiv(3, 3)*gderiv(j, 1)+gdispderiv(3, 1)*gderiv(j, 3)
617 call getGlobalDeriv(etype, nn, naturalcoord, elem1, det, gderiv(1:nn,1:3))
624 gderiv(nn+1, :) = -2.0D0*naturalcoord(1)*inverse1(1, :)/det
625 gderiv(nn+2, :) = -2.0D0*naturalcoord(2)*inverse1(2, :)/det
626 gderiv(nn+3, :) = -2.0D0*naturalcoord(3)*inverse1(3, :)/det
630 B(1, 3*J-2) = gderiv(j, 1)
631 B(2, 3*J-1) = gderiv(j, 2)
632 B(3, 3*J ) = gderiv(j, 3)
633 B(4, 3*J-2) = gderiv(j, 2)
634 B(4, 3*J-1) = gderiv(j, 1)
635 B(5, 3*J-1) = gderiv(j, 3)
636 B(5, 3*J ) = gderiv(j, 2)
637 B(6, 3*J-2) = gderiv(j, 3)
638 B(6, 3*J ) = gderiv(j, 1)
701 real(kind=kreal) :: naturalcoord(3), gderiv(nn+3, 3) local
738 call getGlobalDeriv( etype, nn, naturalcoord, ecoord, det, gderiv(1:nn, 1:3) )
758 gderiv(nn+1, :) = -2.0D0*naturalcoord(1)*inverse(1, :)/det
759 gderiv(nn+2, :) = -2.0D0*naturalcoord(2)*inverse(2, :)/det
760 gderiv(nn+3, :) = -2.0D0*naturalcoord(3)*inverse(3, :)/det
764 B(1, 3*j-2) = gderiv(j, 1)
765 B(2, 3*j-1) = gderiv(j, 2)
766 B(3, 3*j ) = gderiv(j, 3)
767 B(4, 3*j-2) = gderiv(j, 2)
768 B(4, 3*j-1) = gderiv(j, 1)
769 B(5, 3*j-1) = gderiv(j, 3)
770 B(5, 3*j ) = gderiv(j, 2)
771 B(6, 3*j-2) = gderiv(j, 3)
772 B(6, 3*j ) = gderiv(j, 1)