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

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

!||    r13ke3         ../engine/source/elements/spring/r13ke3.F

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

!||    imp_glob_k     ../engine/source/implicit/imp_glob_k.F

!||    imp_glob_k0    ../engine/source/implicit/imp_glob_k.F

!||--- calls      -----------------------------------------------------

!||    assem_p        ../engine/source/implicit/assem_p.F

!||    peoff          ../engine/source/elements/beam/peoff.F

!||    pkeg3          ../engine/source/elements/beam/pke3.F

!||    r13mat3        ../engine/source/elements/spring/r13mat3.F

!||    r13sumg3       ../engine/source/elements/spring/r13sumg3.F

!||    r4coork3       ../engine/source/elements/spring/r4coork3.F

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

!||    elbufdef_mod   ../common_source/modules/mat_elem/elbufdef_mod.F90

!||    element_mod    ../common_source/modules/elements/element_mod.F90

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


      SUBROUTINE r13ke3(JFT    ,JLT   ,NEL   ,MTN    ,PM    ,

     1                  GEO    ,IXR   ,X     ,ELBUF_TAB,NPF   ,

     2                  TF  ,SKEW  ,OFFSET,FR_WAVE,IKGEO  ,IGEO,

     3                  ETAG  , IDDL  ,NDOF  ,K_DIAG ,K_LT  ,

     4                  IADK  ,JDIK  )


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

C   M o d u l e s

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

      USE elbufdef_mod

      use element_mod , only : nixr

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

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

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

#include      "implicit_f.inc"

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

C   G l o b a l   P a r a m e t e r s

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

#include      "mvsiz_p.inc"

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

C   C o m m o n   B l o c k s

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

#include      "com04_c.inc"

#include      "param_c.inc"

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

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

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

      INTEGER IXR(NIXR,*), NPF(*),JFT,JLT,OFFSET,NEL,MTN,IKGEO

      INTEGER IGEO(*),

     .         ETAG(*),IDDL(*)  ,NDOF(*)  ,IADK(*) ,JDIK(*)


C     REAL

      my_real

     .   geo(npropg,*),x(*),tf(*),skew(lskew,*),

     .   fr_wave(*),pm(*),ke11(36,mvsiz),ke12(36,mvsiz),ke22(36,mvsiz),

     .   off(mvsiz),k_diag(*) ,k_lt(*)

C

      TYPE (ELBUF_STRUCT_), TARGET:: ELBUF_TAB

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

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

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

C     REAL

      my_real

     .     al(mvsiz),kx(mvsiz),ky(mvsiz) ,kz(mvsiz),

     .     mx(mvsiz)  ,my(mvsiz)     ,mz(mvsiz),

     .    r11(mvsiz),r12(mvsiz),r13(mvsiz),

     .    r21(mvsiz),r22(mvsiz),r23(mvsiz),

     .    r31(mvsiz),r32(mvsiz),r33(mvsiz)

      INTEGER  PID(MVSIZ),NGL(MVSIZ)

      INTEGER

     .

     .

     .

     .

     .    i,ii(3)

C

      TYPE(g_bufel_)  ,POINTER :: GBUF

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

      GBUF => elbuf_tab%GBUF

C

C

      DO i=1,3

        ii(i) = (i-1)*nel + 1

      ENDDO

C

C      I0 = IXR(1,1)

C      IGTYP =  NINT(GEO(12,I0))

c      NB1 =1

c      NB2 =NB1  + NEL

c      NB1 =NB1  + OFFSET

c      NEL3 = 3 * NEL

c      NEL4 = 4 * NEL

c      OFFSET3 = 3 * OFFSET

c      OFFSET4 = 4 * OFFSET

C

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

C      IF(IGTYP==13)THEN

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

c       NB3 =NB2  + NEL3

c       NB4 =NB3  + NEL

c       NB5 =NB4  + NEL3

c       IF (SCODVER<101) THEN

c         NB6 =NB5  + NEL3

c       ELSE

c         NB6 =NB5  + 2*NEL3

c       END IF

c       NB7 =NB6  + NEL3

c       NB8 =NB7  + NEL3

c       NB9 =NB8  + NEL3

c       NB10=NB9  + NEL3

c       NB11=NB10 + NEL3

c       NB12=NB11 + NEL3

c       NB13=NB12 + NEL3

c       NB14=NB13 + NEL3

c       IF (SCODVER<101) THEN

c         NB15 = NB14 + NEL3

c       ELSE

c         NB15 = NB14 + 2*NEL3

c       END IF

c       NBFI = NB15 + NEL4*6

c       NB2A=NB2  + NEL

c       NB2B=NB2A + NEL

c       NB4A=NB4  + NEL

c       NB4B=NB4A + NEL

c       NB5A=NB5  + NEL

c       NB5B=NB5A + NEL

c       NB6A=NB6  + NEL

c       NB6B=NB6A + NEL

c       NB7A=NB7  + NEL

c       NB7B=NB7A + NEL

c       NB8A=NB8  + NEL

c       NB8B=NB8A + NEL

c       NB9A=NB9  + NEL

c       NB9B=NB9A + NEL

c       NB10A=NB10  + NEL

c       NB10B=NB10A + NEL

c       NB11A=NB11 + NEL

c       NB11B=NB11A+ NEL

c       NB12A=NB12 + NEL

c       NB12B=NB12A+ NEL

c       NB13A=NB13 + NEL

c       NB13B=NB13A+ NEL

c       NB15A=NB15 + NEL4

c       NB15B=NB15A+ NEL4

c       NB15C=NB15B+ NEL4

c       NB15D=NB15C+ NEL4

c       NB15E=NB15D+ NEL4

c       NB16 =NB15E+ NEL4

c       NBFI =NB16 + 6*NEL

C

c       NB2 =NB2  + OFFSET

c       NB3 =NB3  + OFFSET

c       NB4 =NB4  + OFFSET

c       NB5 =NB5  + OFFSET

c       NB6 =NB6  + OFFSET

c       NB7 =NB7  + OFFSET

c       NB8 =NB8  + OFFSET

c       NB9 =NB9  + OFFSET

c       NB10=NB10 + OFFSET

c       NB11=NB11 + OFFSET

c       NB12=NB12 + OFFSET

c       NB13=NB13 + OFFSET

c       NB14=NB14 + OFFSET3

c       NB2A=NB2A + OFFSET

c       NB2B=NB2B + OFFSET

c       NB4A=NB4A + OFFSET

c       NB4B=NB4B + OFFSET

c       NB5A=NB5A + OFFSET

c       NB5B=NB5B + OFFSET

c       NB6A=NB6A + OFFSET

c       NB6B=NB6B + OFFSET

c       NB7A=NB7A + OFFSET

c       NB7B=NB7B + OFFSET

c       NB8A=NB8A + OFFSET

c       NB8B=NB8B + OFFSET

c       NB9A=NB9A + OFFSET

c       NB9B=NB9B + OFFSET

c       NB10A=NB10A + OFFSET

c       NB10B=NB10B + OFFSET

c       NB11A=NB11A+ OFFSET

c       NB11B=NB11B+ OFFSET

c       NB12A=NB12A+ OFFSET

c       NB12B=NB12B+ OFFSET

c       NB13A=NB13A+ OFFSET

c       NB13B=NB13B+ OFFSET

c       NB15 =NB15  + OFFSET4

c       NB15A=NB15A + OFFSET4

c       NB15B=NB15B + OFFSET4

c       NB15C=NB15C + OFFSET4

c       NB15D=NB15D + OFFSET4

c       NB15E=NB15E + OFFSET4

C    FX,FY,FZ       -> NB2,NB2A,NB2B        E               -> NB3

C    DX,DY,DZ       -> NB4,NB4A,NB4B        OFF             -> NB1

C    DPX,DPY,DPZ    -> NB6,NB6A,NB6B        DPX2,DPY2,DPZ2  -> NB7,NB7A,NB7B

C    FXEP,FYEP,FZEP -> NB12,NB12A,NB12B     X0,Y0,Z0        -> NB5,NB5A,NB5B

C    XMOM,YMOM,ZMOM -> NB8,NB8A,NB8B        RX,RY,RZ        -> NB9,NB9A,NB9B

C    RPX,RPY,RPZ    -> NB10,NB10A,NB10B     XMEP,YMEP,ZMEP  -> NB11,NB11A,NB11B

C    RPX2,RPY2,RPZ2 -> NB13,NB13A,NB13B

C    RLOC           -> NB14

      CALL r4coork3(jft    ,jlt    ,x      ,ixr   ,pid    ,

     2              ngl    ,al     ,r11    ,r12   ,r13    ,

     3              r21    ,r22    ,r23    ,r31   ,r32    ,

     4              r33   ,gbuf%SKEW,gbuf%OFF,off   )

      CALL r13mat3(jft    ,jlt    ,geo    ,kx      ,ky     ,

     2             kz     ,mx     ,my     ,mz      ,pid    ,

     3             gbuf%LENGTH(ii(1)),gbuf%FOR(ii(1)),gbuf%FOR(ii(2)),gbuf%FOR(ii(3)),gbuf%TOTDEPL(ii(1)),

     4             gbuf%TOTDEPL(ii(2)),gbuf%TOTDEPL(ii(3)),gbuf%MOM(ii(1)),gbuf%MOM(ii(2)),gbuf%MOM(ii(3)),

     5             gbuf%TOTROT(ii(1)),gbuf%TOTROT(ii(2)),gbuf%TOTROT(ii(3)),tf     ,npf     ,

     6             gbuf%POSX,gbuf%POSY,gbuf%POSZ,gbuf%POSXX,gbuf%POSYY,

     7             gbuf%POSZZ,igeo  )

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

C     STIFFNESS MATRIX

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

      CALL r13sumg3(jft    ,jlt    ,al      ,kx      ,ky     ,

     2              kz     ,mx     ,my      ,mz      ,r11    ,

     3              r12    ,r13    ,r21     ,r22     ,r23    ,

     4              r31    ,r32    ,r33     ,ke11    ,ke12   ,

     5              ke22   )

      IF (ikgeo==1) CALL pkeg3(jft    ,jlt    ,

     1                gbuf%FOR(ii(1))  ,al      ,ke11    ,ke12    ,ke22   )


              IF (neig>0) CALL peoff(

     1                   nixr, jft, jlt, ixr, etag,

     2                   off )

              CALL assem_p(nixr,6   ,

     1         ixr       ,nel       ,iddl      ,ndof      ,k_diag    ,

     2         k_lt      ,iadk      ,jdik      ,ke11      ,ke12      ,

     3         ke22      ,off       )


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

      RETURN


      END

assem_p
subroutine assem_p(nixpl, nd, ixp, nel, iddl, ndof, k_diag, k_lt, iadk, jdik, kc11, kc12, kc22, off)
Definition assem_p.F:40

my_real
#define my_real
Definition cppsort.cpp:32

peoff
subroutine peoff(nixpl, jft, jlt, ixp, etag, off)
Definition peoff.F:35

pkeg3
subroutine pkeg3(jft, jlt, fx, al, ke11, ke12, ke22)
Definition pke3.F:131

r13ke3
subroutine r13ke3(jft, jlt, nel, mtn, pm, geo, ixr, x, elbuf_tab, npf, tf, skew, offset, fr_wave, ikgeo, igeo, etag, iddl, ndof, k_diag, k_lt, iadk, jdik)
Definition r13ke3.F:44

r13mat3
subroutine r13mat3(jft, jlt, geo, kx, ky, kz, mx, my, mz, mgn, al0, fx, fy, fz, dx, dy, dz, xmom, ymom, zmom, rx, ry, rz, tf, npf, posx, posy, posz, posxx, posyy, poszz, igeo)
Definition r13mat3.F:38

r13sumg3
subroutine r13sumg3(jft, jlt, al, kx, ky, kz, mx, my, mz, r11, r12, r13, r21, r22, r23, r31, r32, r33, ke11, ke12, ke22)
Definition r13sumg3.F:33

r4coork3
subroutine r4coork3(jft, jlt, x, ncc, pid, ngl, al, r11, r12, r13, r21, r22, r23, r31, r32, r33, rloc, offg, off)
Definition r4coork3.F:34