qrota_vect.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	qrota_vect (x, ixq, kcvt, vect, gama, nel)

Function/Subroutine Documentation

qrota_vect()

subroutine qrota_vect	(		x,
		integer, dimension(nixq,*)	ixq,
		integer	kcvt,
			vect,
			gama,
		integer	nel )

Definition at line 30 of file qrota_vect.F.

      use element_mod , only : nixq
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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real
     .   x(3,*),vect(3,*),gama(6,*)
      INTEGER IXQ(NIXQ,*), KCVT, NEL
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real
     .   y1(mvsiz), y2(mvsiz), y3(mvsiz), y4(mvsiz),
     .   z1(mvsiz), z2(mvsiz), z3(mvsiz), z4(mvsiz),
     .   r22(mvsiz),r23(mvsiz),r32(mvsiz),r33(mvsiz),
     .   sy(mvsiz),sz(mvsiz),ty(mvsiz),tz(mvsiz),
     .   ct,cs,suma,s1,s2,s3
      INTEGER NC1(MVSIZ), NC2(MVSIZ), NC3(MVSIZ), NC4(MVSIZ)
      INTEGER I
C-----------------------------------------------
      IF (kcvt == 0) THEN
        DO i=1,nel
          nc1(i)=ixq(2,i)
          nc2(i)=ixq(3,i)
          nc3(i)=ixq(4,i)
          nc4(i)=ixq(5,i)
C----------------------------
C     NODE COORDINATES
C----------------------------
          y1(i)=x(2,nc1(i))
          z1(i)=x(3,nc1(i))
          y2(i)=x(2,nc2(i))
          z2(i)=x(3,nc2(i))
          y3(i)=x(2,nc3(i))
          z3(i)=x(3,nc3(i))
          y4(i)=x(2,nc4(i))
          z4(i)=x(3,nc4(i))
C---------------------------------------
C     LOCAL REFERENCE FRAME (ISOPARAMETRIC)
C---------------------------------------
          sy(i)=half*(y2(i)+y3(i)-y1(i)-y4(i))
          sz(i)=half*(z2(i)+z3(i)-z1(i)-z4(i))
          ty(i)=half*(y3(i)+y4(i)-y1(i)-y2(i))
          tz(i)=half*(z3(i)+z4(i)-z1(i)-z2(i))
        ENDDO
C-----------
C       convected frame
C-----------
        DO i=1,nel
          ct = ty(i)*ty(i)+tz(i)*tz(i)
          cs = sy(i)*sy(i)+sz(i)*sz(i)
          IF(cs /= zero) THEN
            suma = sqrt(ct/max(em20,cs))
            sy(i) = sy(i)*suma + tz(i)
            sz(i) = sz(i)*suma - ty(i)
          ELSEIF(ct /= zero)THEN
            suma = sqrt(cs/max(em20,ct))
            sy(i) = sy(i) + tz(i)*suma
            sz(i) = sz(i) - ty(i)*suma
          END IF 
          suma=one/max(sqrt(sy(i)**2+sz(i)**2),em20)
          sy(i)=sy(i)*suma
          sz(i)=sz(i)*suma
C-----------
C       TRANSFORMATION MATRIX GLOBAL -> CONVECTED
C-----------
          r22(i)= sy(i)
          r32(i)=-sz(i)
          r23(i)= sz(i)
          r33(i)= sy(i)
        ENDDO
      END IF
c
      DO i=1,nel
C-----------
C       change of reference frame
C-----------
        s1=vect(1,i)
        s2=vect(2,i)
        s3=vect(3,i)
        vect(1,i)=s2*r22(i)+s3*r23(i)
        vect(2,i)=s2*r32(i)+s3*r33(i)
        vect(3,i)=zero
      ENDDO
C-----------
      RETURN