checksvolume_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/.


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

!||    checkvolume_8n   ../starter/source/elements/solid/solide/checksvolume.F

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

!||    hm_read_solid    ../starter/source/elements/reader/hm_read_solid.F

!||    s16coor3         ../starter/source/elements/thickshell/solide16/s16coor3.F

!||    s20coor3         ../starter/source/elements/solid/solide20/s20coor3.F

!||    sccoor3          ../starter/source/elements/thickshell/solidec/sccoor3.f

!||    scoor3           ../starter/source/elements/solid/solide/scoor3.F

!||    spcoor3          ../starter/source/elements/solid/sconnect/spcoor3.f

!||    srcoor3          ../starter/source/elements/solid/solide/srcoor3.f

!||--- uses       -----------------------------------------------------

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


      FUNCTION checkvolume_8n(X,IXS)

      use element_mod , only : nixs

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

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

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

#include      "implicit_f.inc"

      my_real :: checkvolume_8n

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

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

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

      INTEGER :: ixs(nixs)

      my_real :: x(3,*)

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

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

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

      my_real

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

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

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

     .   x17, x28, x35, x46, y17, y28, y35, y46, z17, z28, z35, z46,

     .   jac1,jac2,jac3,jac4,jac5,jac6,jac7,jac8,jac9,

     .   a17,a28,b17,b28,c17,c28,

     .   jac_59_68, jac_67_49, jac_48_57

C=======================================================================

      checkvolume_8n = zero

C

      x1=x(1,ixs(2))

      y1=x(2,ixs(2))

      z1=x(3,ixs(2))

      x2=x(1,ixs(3))

      y2=x(2,ixs(3))

      z2=x(3,ixs(3))

      x3=x(1,ixs(4))

      y3=x(2,ixs(4))

      z3=x(3,ixs(4))

      x4=x(1,ixs(5))

      y4=x(2,ixs(5))

      z4=x(3,ixs(5))

      x5=x(1,ixs(6))

      y5=x(2,ixs(6))

      z5=x(3,ixs(6))

      x6=x(1,ixs(7))

      y6=x(2,ixs(7))

      z6=x(3,ixs(7))

      x7=x(1,ixs(8))

      y7=x(2,ixs(8))

      z7=x(3,ixs(8))

      x8=x(1,ixs(9))

      y8=x(2,ixs(9))

      z8=x(3,ixs(9))

C

      x17=x7-x1

      x28=x8-x2

      x35=x5-x3

      x46=x6-x4

      y17=y7-y1

      y28=y8-y2

      y35=y5-y3

      y46=y6-y4

      z17=z7-z1

      z28=z8-z2

      z35=z5-z3

      z46=z6-z4

C

      jac1=x17+x28-x35-x46

      jac2=y17+y28-y35-y46

      jac3=z17+z28-z35-z46

      a17=x17+x46

      a28=x28+x35

      b17=y17+y46

      b28=y28+y35

      c17=z17+z46

      c28=z28+z35

      jac4=a17+a28

      jac5=b17+b28

      jac6=c17+c28

      jac7=a17-a28

      jac8=b17-b28

      jac9=c17-c28

      jac_59_68=jac5*jac9-jac6*jac8

      jac_67_49=jac6*jac7-jac4*jac9

      jac_48_57=jac4*jac8-jac5*jac7

      checkvolume_8n=one_over_64*(jac1*jac_59_68+jac2*jac_67_49+jac3*jac_48_57)


      RETURN


      END

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

!||    checkvolume_6n   ../starter/source/elements/solid/solide/checksvolume.F

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

!||    hm_read_solid    ../starter/source/elements/reader/hm_read_solid.F

!||    s6ccoor3         ../starter/source/elements/thickshell/solide6c/s6ccoor3.F

!||--- uses       -----------------------------------------------------

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


      FUNCTION checkvolume_6n(X,IXS)

      use element_mod , only : nixs

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

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

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

#include      "implicit_f.inc"

      my_real :: checkvolume_6n

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

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

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

      INTEGER :: ixs(nixs)

      my_real :: x(3,*)

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

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

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

      my_real

     .   x1,x2,x3,x4,x5,x6,y1,y2,y3,y4,y5,y6,z1,z2,z3,z4,z5,z6,

     .   x21,x31,x41,x54,x64,y21,y31,y41,y54,y64,z21,z31,z41,z54,z64,

     .   jac1,jac2,jac3,jac4,jac5,jac6,jac7,jac8,jac9, jac_59_68,jac_67_49,jac_48_57

C=======================================================================

      checkvolume_6n = zero

C

      x1=x(1,ixs(2))

      y1=x(2,ixs(2))

      z1=x(3,ixs(2))

      x2=x(1,ixs(3))

      y2=x(2,ixs(3))

      z2=x(3,ixs(3))

      x3=x(1,ixs(4))

      y3=x(2,ixs(4))

      z3=x(3,ixs(4))

      x4=x(1,ixs(6))

      y4=x(2,ixs(6))

      z4=x(3,ixs(6))

      x5=x(1,ixs(7))

      y5=x(2,ixs(7))

      z5=x(3,ixs(7))

      x6=x(1,ixs(8))

      y6=x(2,ixs(8))

      z6=x(3,ixs(8))

C

      x21=x2-x1

      x31=x3-x1

      x41=x4-x1

      x54=x5-x4

      x64=x6-x4

      y21=y2-y1

      y31=y3-y1

      y41=y4-y1

      y54=y5-y4

      y64=y6-y4

      z21=z2-z1

      z31=z3-z1

      z41=z4-z1

      z54=z5-z4

      z64=z6-z4

C

      jac1=x21+x54

      jac2=y21+y54

      jac3=z21+z54

      jac4=x31+x64

      jac5=y31+y64

      jac6=z31+z64

      jac7=third*(x41+x5-x2+x6-x3)

      jac8=third*(y41+y5-y2+y6-y3)

      jac9=third*(z41+z5-z2+z6-z3)

C

      jac_59_68=jac5*jac9-jac6*jac8

      jac_67_49=jac6*jac7-jac4*jac9

      jac_48_57=jac4*jac8-jac5*jac7

C

      checkvolume_6n=one_over_8*(jac1*jac_59_68+jac2*jac_67_49+jac3*jac_48_57)

      RETURN


      END

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

!||    checkvolume_4n   ../starter/source/elements/solid/solide/checksvolume.F

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

!||    hm_read_solid    ../starter/source/elements/reader/hm_read_solid.f

!||    s10coor3         ../starter/source/elements/solid/solide10/s10coor3.F

!||    s10coor3_old     ../starter/source/elements/solid/solide10/s10coor3.F

!||    s4coor3          ../starter/source/elements/solid/solide4/s4coor3.F

!||    sgsavini         ../starter/source/elements/solid/solide/scoor3.F

!||    soltosphx4       ../starter/source/elements/sph/soltosph.F

!||--- uses       -----------------------------------------------------

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


      FUNCTION checkvolume_4n(X,IXS)

      use element_mod , only : nixs

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

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

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

#include      "implicit_f.inc"

      my_real :: checkvolume_4n

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

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

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

      INTEGER :: ixs(nixs)

      my_real :: x(3,*)

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

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

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

      my_real

     .   x1,x2,x3,x4,y1,y2,y3,y4,z1,z2,z3,z4,

     .   jac1 ,jac2 ,jac3 ,

     .   jac4 ,jac5 ,jac6 ,

     .   jac7 ,jac8 ,jac9 ,

     .   jac_59_68, jac_67_49, jac_48_57

C=======================================================================

      checkvolume_4n = zero

C

      x1=x(1,ixs(2))

      y1=x(2,ixs(2))

      z1=x(3,ixs(2))

      x2=x(1,ixs(4))

      y2=x(2,ixs(4))

      z2=x(3,ixs(4))

      x3=x(1,ixs(7))

      y3=x(2,ixs(7))

      z3=x(3,ixs(7))

      x4=x(1,ixs(6))

      y4=x(2,ixs(6))

      z4=x(3,ixs(6))

C

      jac1=x1-x4

      jac2=y1-y4

      jac3=z1-z4

C

      jac4=x2-x4

      jac5=y2-y4

      jac6=z2-z4

C

      jac7=x3-x4

      jac8=y3-y4

      jac9=z3-z4

C

      jac_59_68=jac5*jac9-jac6*jac8

      jac_67_49=jac6*jac7-jac4*jac9

      jac_48_57=jac4*jac8-jac5*jac7

C

      checkvolume_4n=(jac1*jac_59_68+jac2*jac_67_49+jac3*jac_48_57)*one_over_6

C

      RETURN


      END

checkvolume_6n
function checkvolume_6n(x, ixs)
Definition checksvolume.F:132

checkvolume_8n
function checkvolume_8n(x, ixs)
Definition checksvolume.F:39

checkvolume_4n
function checkvolume_4n(x, ixs)
Definition checksvolume.F:217

my_real
#define my_real
Definition cppsort.cpp:32

hm_read_solid
subroutine hm_read_solid(ixs, pm, itab, itabm1, ipart, iparts, isolnod, ixs10, ixs20, ixs16, igeo, lsubmodel, is_dyna, x)
Definition hm_read_solid.F:42

sccoor3
subroutine sccoor3(x, ixs, geo, mxt, ngeo, ngl, ix1, ix2, ix3, ix4, ix5, ix6, ix7, ix8, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, rx, ry, rz, sx, sy, sz, tx, ty, tz, r11, r21, r31, r12, r22, r32, r13, r23, r33, f1x, f1y, f1z, f2x, f2y, f2z, temp0, temp, nintemp)
Definition sccoor3.F:43

spcoor3
subroutine spcoor3(x, ixs, geo, nel, mxt, pid, ngl, ix1, ix2, ix3, ix4, ix5, ix6, ix7, ix8, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, volu, thick)
Definition spcoor3.F:41

srcoor3
subroutine srcoor3(x, xrefs, ixs, geo, mxt, ngeo, ngl, jhbe, ix1, ix2, ix3, ix4, ix5, ix6, ix7, ix8, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, f1x, f1y, f1z, f2x, f2y, f2z, temp0, temp, nintemp, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8)
Definition srcoor3.F:52

starter
program starter
Definition starter.F:39