fsigtrini.F

Go to the documentation of this file.

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/.
!||====================================================================
!||    fsigtrini   ../starter/source/constraints/fxbody/fsigtrini.F
!||--- called by ------------------------------------------------------
!||    fxbsini     ../starter/source/constraints/fxbody/fxbsini.F
!||====================================================================
      SUBROUTINE fsigtrini(FXBELM, IPARG , X     , PM  , IXT ,
     .                     GEO   , FXBMOD, FXBSIG, NELT)
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      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER FXBELM(*), IPARG(NPARG,*), IXT(NIXT,*), NELT
      my_real
     .        fxbsig(*), x(3,*), pm(npropm,*), fxbmod(*),
     .        geo(npropg,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IG, OFFSET, LAST, NFT, NFS, I, NG, IEL,
     .        mat, prop, n1, n2
      my_real
     .        area, e, x1, y1, z1, x2, y2, z2, l0, d11, d12, d13,
     .        d21, d22, d23, xc1, yc1, zc1, xc2, yc2, zc2, lc
C-----------------------------------------------
C
      DO ig=1,nelt,mvsiz
         offset=ig-1
         last=min(mvsiz,nelt-offset)
         nft=offset*7
         nfs=offset*2
         DO i=1,last
            ng=fxbelm(nft+7*(i-1)+1)
            iel=iparg(3,ng)+fxbelm(nft+7*(i-1)+2)
            mat=ixt(1,iel)
            prop=ixt(4,iel)
            area=geo(1,prop)
            e=pm(20,mat)
            x1=x(1,ixt(2,iel))
            y1=x(2,ixt(2,iel))
            z1=x(3,ixt(2,iel))
            x2=x(1,ixt(3,iel))
            y2=x(2,ixt(3,iel))
            z2=x(3,ixt(3,iel))
            l0=sqrt((x2-x1)**2+(y2-y1)**2+(z2-z1)**2)
            n1=fxbelm(nft+7*(i-1)+3)
            n2=fxbelm(nft+7*(i-1)+4)
            d11=fxbmod(6*(n1-1)+1)
            d12=fxbmod(6*(n1-1)+2)
            d13=fxbmod(6*(n1-1)+3)
            d21=fxbmod(6*(n2-1)+1)
            d22=fxbmod(6*(n2-1)+2)
            d23=fxbmod(6*(n2-1)+3)
            xc1=x1+d11
            yc1=y1+d12
            zc1=z1+d13
            xc2=x2+d21
            yc2=y2+d22
            zc2=z2+d23
            lc=sqrt((xc2-xc1)**2+(yc2-yc1)**2+(zc2-zc1)**2)
C
            fxbsig(nfs+2*(i-1)+1)=e*area*(lc-l0)
            fxbsig(nfs+2*(i-1)+2)=half*e*area*(lc-l0)**2
         ENDDO
      ENDDO
C
      RETURN
      END SUBROUTINE fsigtrini