s8sansb_8F_source.html

Copyright>        OpenRadioss

Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.

Copyright>

Copyright>        This program is free software: you can redistribute it and/or modify

Copyright>        it under the terms of the GNU Affero General Public License as published by

Copyright>        the Free Software Foundation, either version 3 of the License, or

Copyright>        (at your option) any later version.

Copyright>

Copyright>        This program is distributed in the hope that it will be useful,

Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of

Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

Copyright>        GNU Affero General Public License for more details.

Copyright>

Copyright>        You should have received a copy of the GNU Affero General Public License

Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.

Copyright>

Copyright>

Copyright>        Commercial Alternative: Altair Radioss Software

Copyright>

Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss

Copyright>        software under a commercial license.  Contact Altair to discuss further if the

Copyright>        commercial version may interest you: https://www.altair.com/radioss/.


!||====================================================================

!||    s8sansb   ../engine/source/elements/solid/solide8s/s8sansb.f

!||--- called by ------------------------------------------------------

!||    s8slke3   ../engine/source/elements/solid/solide8s/s8slke3.F

!||====================================================================


      SUBROUTINE s8sansb(

     .    X1,X2,X3,X4,X5,X6,X7,X8,Y1,Y2,Y3,Y4,Y5,Y6,Y7,Y8,

     .    Z1,Z2,Z3,Z4,Z5,Z6,Z7,Z8, INVJ,KSI,ETA,ZETA,

     .    A11,A12,A13,A21,A22,A23,A31,A32,A33,

     .    DN_R, DN_S, DN_T,BB)

C-----------------------------------------------

C   I m p l i c i t   T y p e s

C-----------------------------------------------

#include      "implicit_f.inc"

C-----------------------------------------------

C   D u m m y   A r g u m e n t s

C-----------------------------------------------

      INTEGER I

C     REAL

      my_real

     .   DN_R(8), DN_S(8), DN_T(8), INVJ(9),

     .   x1,x2,x3,x4,x5,x6,x7,x8,

     .   y1,y2,y3,y4,y5,y6,y7,y8,

     .   z1,z2,z3,z4,z5,z6,z7,z8,

     .   a11, a12, a13,

     .   a21, a22, a23,

     .   a31, a32, a33,

     .   ksi,eta,zeta, bb(6,24)

C-----------------------------------------------

C   L o c a l   V a r i a b l e s

C-----------------------------------------------

      INTEGER J,JJ,L,M,LM

C     REAL

      my_real

     .   B(6,24),C(6,6),UM,UP,

     .   a1a1,a1a2,a1a3,a2a1,a2a2,a2a3,

     .   a1b1,a1b2,a1b3,a2b1,a2b2,a2b3,

     .   b1a1,b1a2,b1a3,b2a1,b2a2,b2a3,

     .   b1b1,b1b2,b1b3,b2b1,b2b2,b2b3,

     .   c1a1,c1a2,c1a3,c2a1,c2a2,c2a3,

     .   c1b1,c1b2,c1b3,c2b1,c2b2,c2b3,

     .   ca11,ca21,ca31,ca12,ca22,ca32,ca13,ca23,ca33

C

C-----------------------------------------------

        DO j=1,8

         jj = 3*(j-1)

         b(1,jj+1) =  a11 *dn_t(j)

         b(1,jj+2) =  a12 *dn_t(j)

         b(1,jj+3) =  a13 *dn_t(j)

         b(2,jj+1) =  a21 *dn_r(j)

         b(2,jj+2) =  a22 *dn_r(j)

         b(2,jj+3) =  a23 *dn_r(j)

         b(3,jj+1) =  a31 *dn_s(j)

         b(3,jj+2) =  a32 *dn_s(j)

         b(3,jj+3) =  a33 *dn_s(j)

        ENDDO

! sortir cette partie car independante du PG

        a1a1 = (x2 +x3 -x1 -x4 )*fourth

        a1a2 = (y2 +y3 -y1 -y4 )*fourth

        a1a3 = (z2 +z3 -z1 -z4 )*fourth

        a2a1 = (x4 +x3 -x1 -x2 )*fourth

        a2a2 = (y4 +y3 -y1 -y2 )*fourth

        a2a3 = (z4 +z3 -z1 -z2 )*fourth

        a1b1 = (x6 +x7 -x5 -x8 )*fourth

        a1b2 = (y6 +y7 -y5 -y8 )*fourth

        a1b3 = (z6 +z7 -z5 -z8 )*fourth

        a2b1 = (x8 +x7 -x5 -x6 )*fourth

        a2b2 = (y8 +y7 -y5 -y6 )*fourth

        a2b3 = (z8 +z7 -z5 -z6 )*fourth

!

        um = (one-zeta)/eight

        up = (one+zeta)/eight

        b(4,1) = -um*(a1a1+a2a1)

        b(4,2) = -um*(a1a2+a2a2)

        b(4,3) = -um*(a1a3+a2a3)

        b(4,4) =  um*(a2a1-a1a1)

        b(4,5) =  um*(a2a2-a1a2)

        b(4,6) =  um*(a2a3-a1a3)

        b(4,7) =  um*(a2a1+a1a1)

        b(4,8) =  um*(a2a2+a1a2)

        b(4,9) =  um*(a2a3+a1a3)

        b(4,10)=  um*(a1a1-a2a1)

        b(4,11)=  um*(a1a2-a2a2)

        b(4,12)=  um*(a1a3-a2a3)

        b(4,13) =-up*(a1b1+a2b1)

        b(4,14) =-up*(a1b2+a2b2)

        b(4,15) =-up*(a1b3+a2b3)

        b(4,16) = up*(a2b1-a1b1)

        b(4,17) = up*(a2b2-a1b2)

        b(4,18) = up*(a2b3-a1b3)

        b(4,19) = up*(a2b1+a1b1)

        b(4,20) = up*(a2b2+a1b2)

        b(4,21) = up*(a2b3+a1b3)

        b(4,22) = up*(a1b1-a2b1)

        b(4,23) = up*(a1b2-a2b2)

        b(4,24) = up*(a1b3-a2b3)


! sortir cette partie car independante du PG

        b1a1 = (x2 +x6 -x1 -x5 )*fourth

        b1a2 = (y2 +y6 -y1 -y5 )*fourth

        b1a3 = (z2 +z6 -z1 -z5 )*fourth

        b2a1 = (x5 +x6 -x1 -x2 )*fourth

        b2a2 = (y5 +y6 -y1 -y2 )*fourth

        b2a3 = (z5 +z6 -z1 -z2 )*fourth

        b1b1 = (x3 +x7 -x4 -x8 )*fourth

        b1b2 = (y3 +y7 -y4 -y8 )*fourth

        b1b3 = (z3 +z7 -z4 -z8 )*fourth

        b2b1 = (x7 +x8 -x3 -x4 )*fourth

        b2b2 = (y7 +y8 -y3 -y4 )*fourth

        b2b3 = (z7 +z8 -z3 -z4 )*fourth

!

        um = (one-eta)/eight

        up = (one+eta)/eight

        b(5,1) = -um*(b1a1+b2a1)

        b(5,2) = -um*(b1a2+b2a2)

        b(5,3) = -um*(b1a3+b2a3)

        b(5,4) =  um*(b2a1-b1a1)

        b(5,5) =  um*(b2a2-b1a2)

        b(5,6) =  um*(b2a3-b1a3)

        b(5,7) =  up*(b2b1-b1b1)

        b(5,8) =  up*(b2b2-b1b2)

        b(5,9) =  up*(b2b3-b1b3)

        b(5,10)= -up*(b1b1+b2b1)

        b(5,11)= -up*(b1b2+b2b2)

        b(5,12)= -up*(b1b3+b2b3)

        b(5,13)=  um*(b1a1-b2a1)

        b(5,14)=  um*(b1a2-b2a2)

        b(5,15)=  um*(b1a3-b2a3)

        b(5,16)=  um*(b1a1+b2a1)

        b(5,17)=  um*(b1a2+b2a2)

        b(5,18)=  um*(b1a3+b2a3)

        b(5,19)=  up*(b1b1+b2b1)

        b(5,20)=  up*(b1b2+b2b2)

        b(5,21)=  up*(b1b3+b2b3)

        b(5,22)=  up*(b1b1-b2b1)

        b(5,23)=  up*(b1b2-b2b2)

        b(5,24)=  up*(b1b3-b2b3)


! sortir cette partie car independante du PG

        c1a1 = (x4 +x8 -x1 -x5 )*fourth

        c1a2 = (y4 +y8 -y1 -y5 )*fourth

        c1a3 = (z4 +z8 -z1 -z5 )*fourth

        c2a1 = (x5 +x8 -x1 -x4 )*fourth

        c2a2 = (y5 +y8 -y1 -y4 )*fourth

        c2a3 = (z5 +z8 -z1 -z4 )*fourth

        c1b1 = (x3 +x7 -x2 -x6 )*fourth

        c1b2 = (y3 +y7 -y2 -y6 )*fourth

        c1b3 = (z3 +z7 -z2 -z6 )*fourth

        c2b1 = (x6 +x7 -x2 -x3 )*fourth

        c2b2 = (y6 +y7 -y2 -y3 )*fourth

        c2b3 = (z6 +z7 -z2 -z3 )*fourth

!

        um = (one-ksi)/eight

        up = (one+ksi)/eight

        b(6,1) = -um*(c1a1+c2a1)

        b(6,2) = -um*(c1a2+c2a2)

        b(6,3) = -um*(c1a3+c2a3)

        b(6,4) = -up*(c1b1+c2b1)

        b(6,5) = -up*(c1b2+c2b2)

        b(6,6) = -up*(c1b3+c2b3)

        b(6,7) =  up*(c2b1-c1b1)

        b(6,8) =  up*(c2b2-c1b2)

        b(6,9) =  up*(c2b3-c1b3)

        b(6,10)=  um*(c2a1-c1a1)

        b(6,11)=  um*(c2a2-c1a2)

        b(6,12)=  um*(c2a3-c1a3)

        b(6,13)=  um*(c1a1-c2a1)

        b(6,14)=  um*(c1a2-c2a2)

        b(6,15)=  um*(c1a3-c2a3)

        b(6,16)=  up*(c1b1-c2b1)

        b(6,17)=  up*(c1b2-c2b2)

        b(6,18)=  up*(c1b3-c2b3)

        b(6,19)=  up*(c1b1+c2b1)

        b(6,20)=  up*(c1b2+c2b2)

        b(6,21)=  up*(c1b3+c2b3)

        b(6,22)=  um*(c1a1+c2a1)

        b(6,23)=  um*(c1a2+c2a2)

        b(6,24)=  um*(c1a3+c2a3)


        ca11 = invj(3)

        ca21 = invj(1)

        ca31 = invj(2)

        ca12 = invj(6)

        ca22 = invj(4)

        ca32 = invj(5)

        ca13 = invj(9)

        ca23 = invj(7)

        ca33 = invj(8)


        c(1,1) = ca11*ca11

        c(2,1) = ca12*ca12

        c(3,1) = ca13*ca13

        c(4,1) = two*ca11*ca12

        c(5,1) = two*ca11*ca13

        c(6,1) = two*ca12*ca13

        c(1,2) = ca21*ca21

        c(2,2) = ca22*ca22

        c(3,2) = ca23*ca23

        c(4,2) = two*ca21*ca22

        c(5,2) = two*ca21*ca23

        c(6,2) = two*ca22*ca23

        c(1,3) = ca31*ca31

        c(2,3) = ca32*ca32

        c(3,3) = ca33*ca33

        c(4,3) = two*ca31*ca32

        c(5,3) = two*ca31*ca33

        c(6,3) = two*ca32*ca33

        c(1,4) = ca11*ca21

        c(2,4) = ca12*ca22

        c(3,4) = ca13*ca23

        c(4,4) = ca12*ca21+ca11*ca22

        c(5,4) = ca13*ca21+ca11*ca23

        c(6,4) = ca12*ca23+ca13*ca22

        c(1,5) = ca11*ca31

        c(2,5) = ca12*ca32

        c(3,5) = ca13*ca33

        c(4,5) = ca12*ca31+ca11*ca32

        c(5,5) = ca13*ca31+ca11*ca33

        c(6,5) = ca12*ca33+ca13*ca32

        c(1,6) = ca21*ca31

        c(2,6) = ca22*ca32

        c(3,6) = ca23*ca33

        c(4,6) = ca22*ca31+ca21*ca32

        c(5,6) = ca23*ca31+ca21*ca33

        c(6,6) = ca23*ca32+ca22*ca33


        DO l=1,6

         DO m=1,24

          bb(l,m) = zero

          DO lm=1,6

           bb(l,m) = bb(l,m)+c(l,lm)*b(lm,m)

          ENDDO

         ENDDO

        ENDDO

C

      RETURN


      END

s8sansb
subroutine s8sansb(x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, invj, ksi, eta, zeta, a11, a12, a13, a21, a22, a23, a31, a32, a33, dn_r, dn_s, dn_t, bb)
Definition s8sansb.F:35