pfint3.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	pfint3 (for, mom, geo, offg, off, al, f1, f2, f3, m1, m2, m3, sti, stir, nel, pid, f11, f12, f21, f22, f31, f32, m11, m12, m21, m22, m31, m32)

Function/Subroutine Documentation

pfint3()

subroutine pfint3	(		for,
			mom,
			geo,
			offg,
			off,
			al,
			f1,
			f2,
			f3,
			m1,
			m2,
			m3,
			sti,
			stir,
		integer	nel,
		integer, dimension(mvsiz)	pid,
			f11,
			f12,
			f21,
			f22,
			f31,
			f32,
			m11,
			m12,
			m21,
			m22,
			m31,
			m32 )

Definition at line 28 of file pfint3.F.

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 NEL,PID(MVSIZ)
C     REAL
      my_real
     .   for(nel,3),mom(nel,3),geo(npropg,*),offg(*),off(*),al(mvsiz), 
     .   f1(mvsiz),f2(mvsiz),f3(mvsiz),m1(mvsiz),m2(mvsiz),m3(mvsiz),
     .   sti(mvsiz), stir(mvsiz),
     .   f11(mvsiz), f12(mvsiz), f21(mvsiz),  
     .   f22(mvsiz), f31(mvsiz), f32(mvsiz), 
     .   m11(mvsiz), m12(mvsiz), m21(mvsiz),
     .   m22(mvsiz), m31(mvsiz), m32(mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J, IG, IR1Y, IR1Z, IR2Y, IR2Z
C     REAL
      my_real
     .   off_l
C---------------------------------------------------
      DO i=1,nel
        f11(i)= -f1(i)
        f21(i)= -f2(i)
        f31(i)= -f3(i)
C
        f12(i)=-f11(i)
        f22(i)=-f21(i)
        f32(i)=-f31(i)
      ENDDO
C
C     CONTRIBUTION AU MOMENT DU CISAILLEMENT
C
      DO i=1,nel
        m11(i)=-m1(i)
        m21(i)=-m2(i)+half*al(i)*f3(i)
        m31(i)=-m3(i)-half*al(i)*f2(i)
C
        m12(i)=-m11(i)
        m22(i)=m2(i)+half*al(i)*f3(i)
        m32(i)=m3(i)-half*al(i)*f2(i)
      ENDDO
C---------------------------------------------------
C     LIBERATION DES ROTATIONS
C---------------------------------------------------
      DO i=1,nel
        ig=pid(i)
        ir1y=nint(geo(8 ,ig))
        ir1z=nint(geo(9 ,ig))
        ir2y=nint(geo(10,ig))
        ir2z=nint(geo(11,ig))
        m21(i)=ir1y*m21(i)+(1-ir2y)*m22(i)
        m22(i)=ir2y*m22(i)+(1-ir1y)*m21(i)
        m31(i)=ir1z*m31(i)+(1-ir2z)*m32(i)
        m32(i)=ir2z*m32(i)+(1-ir1z)*m31(i)
      ENDDO
C
      off_l = zero
      DO i=1,nel
        IF (off(i) < one) offg(i) = off(i)
        off_l = min(off_l,offg(i))
      ENDDO
      IF (off_l < zero) THEN
        DO i=1,nel
          IF (offg(i) < zero) THEN
            f11(i)=zero
            f21(i)=zero
            f31(i)=zero
            m11(i)=zero
            m21(i)=zero
            m31(i)=zero
            f12(i)=zero
            f22(i)=zero
            f32(i)=zero
            m12(i)=zero
            m22(i)=zero
            m32(i)=zero
            sti(i)=zero
            stir(i)=zero
          ENDIF
        ENDDO
      ENDIF
C---------------------------------------------------
      RETURN