i3sti2.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	i3sti2 (x, irect, stf, ixq, pm, nrt, stfn, nseg, lnsv, nint, nsn, nsv, slsfac, noint, ipm, id, titr, areas, knod2elq, nod2elq, nty, nsns, nsvs, segquadfr)

Function/Subroutine Documentation

i3sti2()

subroutine i3sti2	(		x,
		integer, dimension(4,*)	irect,
			stf,
		integer, dimension(7,*)	ixq,
			pm,
		integer	nrt,
			stfn,
		integer, dimension(*)	nseg,
		integer, dimension(*)	lnsv,
		integer	nint,
		integer	nsn,
		integer, dimension(*)	nsv,
			slsfac,
		integer	noint,
		integer, dimension(npropmi,*)	ipm,
		integer	id,
		character(len=nchartitle)	titr,
			areas,
		integer, dimension(*)	knod2elq,
		integer, dimension(*)	nod2elq,
		integer	nty,
		integer	nsns,
		integer, dimension(*)	nsvs,
		integer, dimension(2,*)	segquadfr )

Definition at line 37 of file i3sti2.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_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      "com01_c.inc"
#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, NTY, NSNS, IPM(NPROPMI,*)
      my_real
     .   slsfac
      INTEGER IRECT(4,*), IXQ(7,*), NSEG(*), LNSV(*), NSV(*),
     .        KNOD2ELQ(*),NOD2ELQ(*), NSVS(*), SEGQUADFR(2,*)
      my_real
     .   x(3,*), stf(*), pm(npropm,*), stfn(*),areas(*)
      INTEGER ID
      CHARACTER(LEN=NCHARTITLE) :: TITR
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, II, NEL, MT, J, NUM, NPT, JJ, LL, IG, IL, IE, INRT,
     .   N1 ,N2 ,STAT, ILINE, LIN, L, N, K
      INTEGER LINES(2,4)
      INTEGER, DIMENSION(:),ALLOCATABLE ::INRTIE
C     REAL
      my_real
     .   area, xl2, ym1, ym2, zm1, zm2,ye(4) ,ze(4),
     .   y1 ,y2 ,z1 ,z2
      DATA lines/1,2,
     .           2,3,
     .           3,4,
     .           4,1/
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
C
       ALLOCATE(inrtie(numelq),stat=stat)
       IF (stat /= 0) CALL ancmsg(msgid=268,anmode=aninfo,
     .                             msgtype=msgerror,
     .                         c1='INRTIE')
        inrtie(1:numelq)=0
C
      DO i=1,nrt
        ii=i
        CALL inrch2(x   , irect,ixq,  ii , nel,
     .              nint, ym1  ,ym2,  zm1,
     .              zm2 , ye   ,ze ,  id , titr)
        IF(nel/=0) THEN
         inrtie(nel) = ii
         CALL inare2(area,ye   ,ze)
         xl2=(ym2-ym1)**2+(zm2-zm1)**2
         mt=ixq(1,nel)
         IF(mt>0)THEN
          stf(i)=slsfac*xl2*pm(32,mt)/area
         ELSE
          stf(i)=zero
          CALL ancmsg(msgid=347,
     .                msgtype=msgwarning,
     .                anmode=aninfo_blind_2,
     .                i1=id,
     .                c1=titr,
     .                i2=ipm(1,mt),
     .                i3=nel,
     .                i4=i)
C
         ENDIF
        ELSE
          stf(i)=zero
        ENDIF
C
        CALL inori2(irect,ii,nel,nint,
     .              noint, ym1, ym2, zm1,zm2 ,
     .              ye   ,ze)
      ENDDO
C---------------------------------------------
C     CALCUL DES RIGIDITES NODALES
C---------------------------------------------
      DO j=1,nsn
        num=nseg(j+1)-nseg(j)
        npt=nseg(j)-1
        DO jj=1,num
          ll=lnsv(npt+jj)
          stfn(j)=stfn(j) + half*stf(ll)
        ENDDO
      ENDDO    
 
C---------------------------------------------
C     SECND NODAL SURFACE COMPUTATION
C---------------------------------------------
C
      IF(nty == 3) THEN
       DO i = 1,nsn    
         areas(i) = zero
         DO j= knod2elq(nsv(i))+1,knod2elq(nsv(i)+1)
            ie = nod2elq(j)
            inrt = inrtie(ie)
            IF(inrt/=0)THEN
               n1=irect(1,inrt)
               n2=irect(2,inrt)
               y1=x(2,n1)
               z1=x(3,n1)
               y2=x(2,n2)
               z2=x(3,n2)
c
               area = sqrt((y2-y1)*(y2-y1)+(z2-z1)*(z2-z1))
               area = area*half
c
              areas(i) = areas(i) + area
           ENDIF
          ENDDO
       ENDDO
      ELSEIF(nty == 5) THEN
       DO i = 1,nsns   
         areas(i) = zero
          DO j= knod2elq(nsvs(i))+1,knod2elq(nsvs(i)+1)
            ie = nod2elq(j)
            lin = -huge(lin)
            DO l=1,4 
              IF(ixq(lines(1,l)+1,ie) ==nsvs(i)) THEN
                 lin = l
                 EXIT
              ENDIF   
            ENDDO
 
           DO k=1,nsegquadfr
              n  =segquadfr(1,k)
              iline=segquadfr(2,k)
              
              IF(n==ie.AND.iline==lin) THEN
 
                n1=ixq(lines(1,iline)+1,n)
                n2=ixq(lines(2,iline)+1,n)
 
               y1=x(2,n1)
               z1=x(3,n1)
               y2=x(2,n2)
               z2=x(3,n2)
 
              area = sqrt((y2-y1)*(y2-y1)+(z2-z1)*(z2-z1))
              area = area*half
 
              areas(i) = areas(i) + area
             ENDIF
           ENDDO
          ENDDO
        ENDDO
      ENDIF
      DEALLOCATE(inrtie)
 
C-----------------------------------------------------
C     MISE DANS IRECT DU NUMERO LOCAL DU NOEUD
C-----------------------------------------------------
      DO i=1,nrt
        ig=irect(1,i)
        CALL local_index(il,ig,nsv,nsn)
        irect(1,i)=il
        ig=irect(2,i)
        CALL local_index(il,ig,nsv,nsn)
        irect(2,i)=il
      ENDDO
C
      RETURN