fxbvini.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/.
!||====================================================================
!||    fxbvini   ../common_source/fxbody/fxbvini.F
!||--- called by ------------------------------------------------------
!||    initia    ../starter/source/elements/initia/initia.F
!||    lecinv    ../engine/source/input/lecinv.F
!||====================================================================
      SUBROUTINE fxbvini(FXBIPM, FXBVIT, FXBRPM, V, VR)
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
#include      "fxbcom.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        INTEGER, INTENT(IN) :: FXBIPM(NBIPM,NFXBODY)
        my_real, INTENT(IN) :: fxbrpm(lenrpm),v(3,numnod),vr(3,numnod)
        my_real, INTENT(INOUT) :: fxbvit(lenvar)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        INTEGER NFX,I,AVAR,ARPM,IDMAST
        my_real p11,p12,p13,p21,p22,p23,p31,p32,p33,w1,w2,w3,wl1,wl2,wl3
 
        DO nfx=1,nfxbody
          idmast=fxbipm(2,nfx)
          avar=fxbipm(13,nfx)
          arpm=fxbipm(14,nfx)
          DO i=1,3
            fxbvit(avar+9+i-1)=v(i,idmast)
          END DO
          p11=fxbrpm(arpm+1)
          p12=fxbrpm(arpm+2)
          p13=fxbrpm(arpm+3)
          p21=fxbrpm(arpm+4)
          p22=fxbrpm(arpm+5)
          p23=fxbrpm(arpm+6)
          p31=fxbrpm(arpm+7)
          p32=fxbrpm(arpm+8)
          p33=fxbrpm(arpm+9)
          w1=vr(1,idmast)
          w2=vr(2,idmast)
          w3=vr(3,idmast)
          wl1=p11*w1+p21*w2+p31*w3
          wl2=p12*w1+p22*w2+p32*w3
          wl3=p13*w1+p23*w2+p33*w3
          fxbvit(avar)  =v(1,idmast)+p12*wl3-p13*wl2
          fxbvit(avar+1)=v(2,idmast)+p22*wl3-p23*wl2
          fxbvit(avar+2)=v(3,idmast)+p32*wl3-p33*wl2
          fxbvit(avar+3)=v(1,idmast)-p11*wl3+p13*wl1
          fxbvit(avar+4)=v(2,idmast)-p21*wl3+p23*wl1
          fxbvit(avar+5)=v(3,idmast)-p31*wl3+p33*wl1
          fxbvit(avar+6)=v(1,idmast)+p11*wl2-p12*wl1
          fxbvit(avar+7)=v(2,idmast)+p21*wl2-p22*wl1
          fxbvit(avar+8)=v(3,idmast)+p31*wl2-p32*wl1
        ENDDO
C
        RETURN
      END