fxbsini.F File Reference

#include "implicit_f.inc"
#include "param_c.inc"
#include "units_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	fxbsini (fxbelm, fxbsig, nels, nelc, neltg, x, iparg, pm, fxbmod, nml, nsn, ixs, ixc, ixtg, geo, fxbrpm, nfx, ifile, lvsig, nsni, nme, ircs, ircm0, nelt, nelp, ixt, ixp, ibeam_vector, rbeam_vector)

Function/Subroutine Documentation

fxbsini()

subroutine fxbsini	(	integer, dimension(*)	fxbelm,
			fxbsig,
		integer	nels,
		integer	nelc,
		integer	neltg,
			x,
		integer, dimension(nparg,*)	iparg,
			pm,
			fxbmod,
		integer	nml,
		integer	nsn,
		integer, dimension(nixs,*)	ixs,
		integer, dimension(nixc,*)	ixc,
		integer, dimension(nixtg,*)	ixtg,
			geo,
			fxbrpm,
		integer	nfx,
		integer	ifile,
		integer	lvsig,
		integer	nsni,
		integer	nme,
		integer	ircs,
		integer	ircm0,
		integer	nelt,
		integer	nelp,
		integer, dimension(nixt,*)	ixt,
		integer, dimension(nixp,*)	ixp,
		integer, dimension(nelp), intent(in)	ibeam_vector,
		dimension(3,nelp), intent(in)	rbeam_vector )

Definition at line 35 of file fxbsini.F.

      use element_mod , only : nixs,nixc,nixtg,nixt,nixp
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      "param_c.inc"
#include      "units_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER FXBELM(*), NELS, NELC, NELTG, IPARG(NPARG,*), NML, 
     .        NSN, IXS(NIXS,*), IXC(NIXC,*), IXTG(NIXTG,*), NFX, IFILE,
     .        LVSIG, NSNI, NME, IRCS, IRCM0, NELT, NELP, IXT(NIXT,*),
     .        IXP(NIXP,*)
      INTEGER, INTENT (IN   ) :: IBEAM_VECTOR(NELP)
      my_real
     .        fxbsig(*), x(3,*), pm(npropm,*), fxbmod(*),
     .        geo(npropg,*), fxbrpm(*)
      my_real, INTENT (IN   ) :: rbeam_vector(3,nelp)     
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,II,IM,IADMOD,IADEL,IADSIG, IADM, IADS, IRCM, IAD, J
      my_real
     .        rini(3,3), vmod(nsn*6), vsig(lvsig), vv(6)
C=======================================================================
      ircm=ircm0
      ircm=ircm+nme*(nsn-nsni)
C      
      DO i=1,3
         DO ii=1,3
            rini(i,ii)=fxbrpm(1+(i-1)*3+ii)
         ENDDO
      ENDDO
C
      iadsig=1
      DO im=1,nml
        IF (ifile==0) THEN
           iadmod=nsn*6*(im-1)+1
           DO i=1,nsn*6
              vmod(i)=fxbmod(iadmod+i-1)
           ENDDO
        ELSEIF (ifile==1) THEN
           iadmod=nsni*6*(im-1)+1
           DO i=1,nsni*6
              vmod(i)=fxbmod(iadmod+i-1)
           ENDDO
           iadm=nsni*6
           DO i=1,nsn-nsni
              ircm=ircm+1
              READ(ifxm,rec=ircm) (vv(ii),ii=1,6)
              DO ii=1,6
                 vmod(iadm+ii)=vv(ii)
              ENDDO
              iadm=iadm+6
           ENDDO
        ENDIF
C   1. Solid elements
        iadel=1
        iads=1
        IF (nels>0) CALL fsigsini(
     .     fxbelm(iadel), iparg, x,          pm,   ixs, 
     .     geo,           vmod,  vsig(iads), rini, nels)
C   2. 4-nodes shell elements
        iadel=iadel+nels*13
        iads=iads+nels*7             
        IF (nelc>0) CALL fsigcini(
     .     fxbelm(iadel), iparg, x,          pm,   ixc, 
     .     geo,           vmod,  vsig(iads), rini, nelc)
C   3. Truss elements
        iadel=iadel+nelc*10
        iads=iads+nelc*10
        IF (nelt>0) CALL fsigtrini(
     .     fxbelm(iadel), iparg, x,         pm,   ixt,
     .     geo,           vmod,  vsig(iads),nelt)
C   4. Beam elements
        iadel=iadel+nelt*7
        iads=iads+nelt*2
        IF (nelp>0) CALL fsigpini(
     .     fxbelm(iadel), iparg, x,          pm,   ixp, 
     .     geo,           vmod,  vsig(iads), rini, nelp,
     .     ibeam_vector ,rbeam_vector)
C   5. 3-nodes shell elements
        iadel=iadel+nelp*9
        iads=iads+nelp*8             
        IF (neltg>0) CALL fsigtini(
     .     fxbelm(iadel), iparg, x,          pm,   ixtg, 
     .     geo,           vmod,  vsig(iads), rini, neltg)
        IF (ifile==0) THEN
           DO i=1,lvsig
              fxbsig(iadsig+i-1)=vsig(i)
           ENDDO
        ELSE
           iad=0
           DO i=1,lvsig/6
              ircs=ircs+1
              WRITE(ifxs,rec=ircs) (vsig(iad+j),j=1,6)
              iad=iad+6
           ENDDO
           ii=lvsig-(lvsig/6)*6
           IF (ii/=0) THEN
              ircs=ircs+1
              WRITE(ifxs,rec=ircs) (vsig(iad+j),j=1,ii),(zero,j=ii+1,6)
           ENDIF
        ENDIF
        iadsig=iadsig+nels*7+nelc*10+nelt*2+nelp*8+neltg*10
      ENDDO
C-----------
      RETURN