czbe3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	czbe3 (jft, jlt, area, a_i, x13, x24, y13, y24, mx13, mx23, mx34, my13, my23, my34, z1, px1, px2, py1, py2, rx, ry, sx, sy, rx2, ry2, sx2, sy2, rhx, rhy, shx, shy, ph1, ph2, hxx, hyy, hxy, nplat, iplat)
subroutine	czber3 (jft, jlt, area, a_i, x13, x24, y13, y24, mx13, mx23, mx34, my13, my23, my34, z1, rx, ry, sx, sy, prx, pry, prxy, prz, phkrx, phkry, phkrxy, pherx, phery, pherxy, phkrz, pherz, phkx, phky, phex, phey, iplat)

Function/Subroutine Documentation

czbe3()

subroutine czbe3	(	integer	jft,
		integer	jlt,
			area,
			a_i,
			x13,
			x24,
			y13,
			y24,
			mx13,
			mx23,
			mx34,
			my13,
			my23,
			my34,
			z1,
			px1,
			px2,
			py1,
			py2,
			rx,
			ry,
			sx,
			sy,
			rx2,
			ry2,
			sx2,
			sy2,
			rhx,
			rhy,
			shx,
			shy,
			ph1,
			ph2,
			hxx,
			hyy,
			hxy,
		integer	nplat,
		integer, dimension(*)	iplat )

Definition at line 28 of file czbe3.F.

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
#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-----------------------------------------------
      INTEGER JFT,JLT,NPLAT,IPLAT(*)
      my_real 
     .   x13(*),x24(*),y13(*),y24(*),area(*),z1(*),a_i(*),
     .   mx13(*),mx23(*),mx34(*),my13(*),my23(*),my34(*)
      my_real 
     .   px1(*),px2(*),py1(*),py2(*),ph1(*),ph2(*),
     .   rx(*),ry(*), sx(*),sy(*),rx2(*),ry2(*), sx2(*),sy2(*),
     .   rhx(mvsiz,4),rhy(mvsiz,4), shx(mvsiz,4),shy(mvsiz,4),hxx(*),hyy(*),hxy(*)
C-----------------------------------------------
C   L O C A L   V A R I A B L E S
C-----------------------------------------------
      INTEGER I,J,EP
C-----------gama(I)=hI/4-PH(I), PH:anti-sym comme bxI------
#include "vectorize.inc"
       DO i=jft,jlt
        ep=iplat(i)
        px1(i) = y24(ep)*a_i(ep)
        px2(i) = -y13(ep)*a_i(ep)
        py1(i) = -x24(ep)*a_i(ep)
        py2(i) =  x13(ep)*a_i(ep)
        ph1(i) =  mx13(ep)*px1(i)+ my13(ep)*py1(i)
        ph2(i) =  mx13(ep)*px2(i)+ my13(ep)*py2(i)
       ENDDO 
C-----------Bc = Bcr0(1+r*rI)+Bcs0(1+s*sI)------
C-----------rhx(I)=rkxk+hkxk*SI---shx(I)=skxk+hkxk*RI-
C----------shear----R : -1 1 1 -1; S: -1 -1 1 1--
#include "vectorize.inc"
       DO i=jft,jlt
        ep=iplat(i)
        sx(i) = mx34(ep)*a_i(ep) 
        sy(i) = my34(ep)*a_i(ep) 
        rx(i) = mx23(ep)*a_i(ep) 
        ry(i) = my23(ep)*a_i(ep) 
        rhx(i,1) = mx23(ep)-mx13(ep) 
        rhx(i,3) = mx23(ep)+mx13(ep) 
        rhy(i,1) = my23(ep)-my13(ep) 
        rhy(i,3) = my23(ep)+my13(ep) 
        shx(i,1) = mx34(ep)-mx13(ep)
        shx(i,2) = mx34(ep)+mx13(ep)
        shy(i,1) = my34(ep)-my13(ep)
        shy(i,2) = my34(ep)+my13(ep)
       ENDDO
C
       DO i=jft,jlt
        sx2(i) = sx(i)*sx(i) 
        sy2(i) = sy(i)*sy(i) 
        rx2(i) = rx(i)*rx(i) 
        ry2(i) = ry(i)*ry(i)
        rhx(i,2) = rhx(i,1) 
        rhx(i,4) = rhx(i,3) 
        rhy(i,2) = rhy(i,1) 
        rhy(i,4) = rhy(i,3) 
        shx(i,3) = shx(i,2)
        shx(i,4) = shx(i,1)
        shy(i,3) = shy(i,2)
        shy(i,4) = shy(i,1)
        hxx(i)=fivep333*(sy2(i)+ry2(i))
        hyy(i)=fivep333*(sx2(i)+rx2(i))
        hxy(i)=-fivep333*(sx(i)*sy(i)+rx(i)*ry(i))
       ENDDO
C
      RETURN

czber3()

subroutine czber3	(	integer	jft,
		integer	jlt,
			area,
			a_i,
			x13,
			x24,
			y13,
			y24,
			mx13,
			mx23,
			mx34,
			my13,
			my23,
			my34,
			z1,
			rx,
			ry,
			sx,
			sy,
			prx,
			pry,
			prxy,
			prz,
			phkrx,
			phkry,
			phkrxy,
			pherx,
			phery,
			pherxy,
			phkrz,
			pherz,
			phkx,
			phky,
			phex,
			phey,
		integer, dimension(*)	iplat )

Definition at line 115 of file czbe3.F.

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
#include      "implicit_f.inc"
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   D U M M Y   A R G U M E N T S
C-----------------------------------------------
      INTEGER JFT,JLT,IPLAT(*)
      my_real 
     .   x13(*),x24(*),y13(*),y24(*),area(*),z1(*),a_i(*),
     .   mx13(*),mx23(*),mx34(*),my13(*),my23(*),my34(*),
     .   rx(*),ry(*), sx(*),sy(*)
      my_real 
     .   prx(4,*),pry(4,*),prxy(4,*),phkrx(4,*),phkry(4,*),
     .   phkrxy(4,*),pherx(4,*),phery(4,*),pherxy(4,*),prz(4,*),
     .   phkrz(4,*),pherz(4,*),phkx(*) ,phky(*) ,phex(*) ,phey(*) 
C-----------------------------------------------
C   L O C A L   V A R I A B L E S
C-----------------------------------------------
      INTEGER I,J,EP
      my_real 
     .   bm0rz(mvsiz,4,4),bmkrz(mvsiz,4,4),bmerz(mvsiz,4,4)
C---------------
        CALL czderirz(jft ,jlt  ,area ,x13 ,x24   ,
     2               y13  ,y24  ,mx13 ,mx23,mx34  ,
     3               my13 ,my23 ,my34 ,bm0rz,bmkrz,
     4               bmerz )
#include "vectorize.inc"
       DO i=jft,jlt
        ep=iplat(i)
        DO j=1,4
         prx(j,i)= bm0rz(ep,1,j)*a_i(ep)
         pry(j,i)= bm0rz(ep,2,j)*a_i(ep)
         prxy(j,i)=bm0rz(ep,3,j)*a_i(ep)
C         
         phkrx(j,i)= bmkrz(ep,1,j)*a_i(ep)
         phkry(j,i)= bmkrz(ep,2,j)*a_i(ep)
         phkrxy(j,i)=bmkrz(ep,3,j)*a_i(ep)
C         
         pherx(j,i)= bmerz(ep,1,j)*a_i(ep)
         phery(j,i)= bmerz(ep,2,j)*a_i(ep)
         pherxy(j,i)=bmerz(ep,3,j)*a_i(ep)
C         
         prz(j,i)=bm0rz(ep,4,j)*a_i(ep)
        ENDDO 
       ENDDO
C       
       DO j=1,4
       DO i=jft,jlt
        ep=iplat(i)
        phkrz(j,i)=bmkrz(ep,4,j)*a_i(ep)
        pherz(j,i)=bmerz(ep,4,j)*a_i(ep)
       ENDDO 
       ENDDO 
C       
       DO i=jft,jlt
        phex(i)= four*sy(i)
        phey(i) =-four*sx(i)
        phkx(i)=-four*ry(i)  
        phky(i)= four*rx(i)
       ENDDO 
C
      RETURN