i9sti3.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	i9sti3 (x, irect, stf, ixs, nrt, nint, nsn, nsv, noint, iele, knod2els, nod2els, igrsurf, isu, ixs10, ixs16, ixs20, id, titr)

Function/Subroutine Documentation

i9sti3()

subroutine i9sti3	(		x,
		integer, dimension(4,*)	irect,
			stf,
		integer, dimension(nixs,*)	ixs,
		integer	nrt,
		integer	nint,
		integer	nsn,
		integer, dimension(*)	nsv,
		integer	noint,
		integer, dimension(*)	iele,
		integer, dimension(*)	knod2els,
		integer, dimension(*)	nod2els,
		type (surf_), dimension(nsurf)	igrsurf,
		integer	isu,
		integer, dimension(*)	ixs10,
		integer, dimension(*)	ixs16,
		integer, dimension(*)	ixs20,
		integer	id,
		character(len=nchartitle)	titr )

Definition at line 35 of file i9sti3.F.

C-----------------------------------------------
C   D e s c r i p t i o n
C-----------------------------------------------
C     NODE AND SEGMENT STIFFNESS
C       specific case : segment belongs to both solid and shell elem ; shell stiffness is used (if > 0)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
      USE groupdef_mod
      USE names_and_titles_mod , ONLY : nchartitle
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      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NRT, NINT, NSN, NOINT,ISU
      INTEGER IRECT(4,*), IXS(NIXS,*),
     .   NSV(*),  IELE(*),
     .   KNOD2ELS(*), NOD2ELS(*), 
     .   IXS10(*), IXS16(*), IXS20(*)
      my_real x(3,*), stf(*)
      INTEGER ID
      CHARACTER(LEN=NCHARTITLE) :: TITR
      TYPE (SURF_)   , DIMENSION(NSURF)   :: IGRSURF
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, INRT, NELS, J, IG, IL
      my_real area
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
 
      DO i=1,nrt
        stf(i)=zero
        inrt=i
C
        IF (isu > 0) THEN
          CALL inelts(x       ,irect  ,ixs     ,nint    ,nels              ,
     .                inrt    ,area   ,noint   ,0       ,igrsurf(isu)%ELTYP,
     .                igrsurf(isu)%ELEM)
          IF(nels/=0) THEN
            iele(i) = nels
            cycle
          END IF
        ENDIF
C----------------------
C     SOLID ELEMS
C----------------------
        CALL insol3(x,irect,ixs,nint,nels,inrt,
     .              area,noint,knod2els ,nod2els ,0 ,ixs10,
     .              ixs16,ixs20)
        IF(nels/=0) THEN
         iele(i) = nels
        ELSE
           IF(nint>0) THEN
              CALL ancmsg(msgid=92,msgtype=msgwarning,anmode=aninfo_blind_2,i1=id,c1=titr,i2=i)
           ENDIF            
           IF(nint<0) THEN
              CALL ancmsg(msgid=93,msgtype=msgwarning,anmode=aninfo_blind_2,i1=id,c1=titr,i2=i)
           ENDIF
        ENDIF
C
      ENDDO
C---------------------------------------------
C
      DO i=1,nrt
        DO j=1,4
          ig=irect(j,i)
          CALL local_index(il,ig,nsv,nsn)
          irect(j,i)=il
        ENDDO
      ENDDO
C
      RETURN