s4rcoor12.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	s4rcoor12 (off, nc1, nc2, nc3, nc4, x, xdp, d, r11, r12, r13, r21, r22, r23, r31, r32, r33, nel, jcvt)

Function/Subroutine Documentation

s4rcoor12()

subroutine s4rcoor12	(		off,
		integer, dimension(*)	nc1,
		integer, dimension(*)	nc2,
		integer, dimension(*)	nc3,
		integer, dimension(*)	nc4,
			x,
		double precision, dimension(3,*)	xdp,
			d,
			r11,
			r12,
			r13,
			r21,
			r22,
			r23,
			r31,
			r32,
			r33,
		integer, intent(in)	nel,
		integer, intent(in)	jcvt )

Definition at line 30 of file s4rcoor12.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      "scr05_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER, INTENT(IN) :: JCVT
      INTEGER NC1(*), NC2(*), NC3(*), NC4(*)
C     REAL
      my_real
     .   off(*),x(3,*),d(3,*),
     .   r11(mvsiz),r12(mvsiz),r13(mvsiz),
     .   r21(mvsiz),r22(mvsiz),r23(mvsiz),
     .   r31(mvsiz),r32(mvsiz),r33(mvsiz)
      DOUBLE PRECISION 
     .   XDP(3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      
      double precision
     .   x1(mvsiz), x2(mvsiz), x3(mvsiz), x4(mvsiz), 
     .   y1(mvsiz), y2(mvsiz), y3(mvsiz), y4(mvsiz), 
     .   z1(mvsiz), z2(mvsiz), z3(mvsiz), z4(mvsiz)
    
C     REAL
      my_real
     .   rx(mvsiz) , ry(mvsiz) , rz(mvsiz) ,
     .   sx(mvsiz) , sy(mvsiz) , sz(mvsiz) ,
     .   tx(mvsiz) , ty(mvsiz) , tz(mvsiz) 
      my_real
     .   x41, y41, z41, x42, y42, z42, x43, y43, z43
C=======================================================================
C---for case of Ismstr=10-> Ismstr11 :
C----1) config is fixed at t= t_dtmin but will not be updated like Ismstr11
C----2) [H] = [B](t=0)*{u}, [B](t=0)<-JAC_I, {u}={d}; 
C----3) in case of JCVT=0, [H] and [D] are transformed to local system and
C-      return to global system for stress
C ---------local system compute----------     
      IF(iresp==1)THEN 
        DO i=1,nel
         IF (off(i) <= one ) cycle 
          x1(i)=xdp(1,nc1(i))  
          y1(i)=xdp(2,nc1(i))  
          z1(i)=xdp(3,nc1(i))  
          x2(i)=xdp(1,nc2(i))  
          y2(i)=xdp(2,nc2(i))  
          z2(i)=xdp(3,nc2(i))  
          x3(i)=xdp(1,nc3(i))  
          y3(i)=xdp(2,nc3(i))  
          z3(i)=xdp(3,nc3(i))  
          x4(i)=xdp(1,nc4(i))  
          y4(i)=xdp(2,nc4(i))  
          z4(i)=xdp(3,nc4(i))  
        ENDDO
      ELSE  
        DO i=1,nel
         IF (off(i) <= one ) cycle 
          x1(i)=x(1,nc1(i))  
          y1(i)=x(2,nc1(i))  
          z1(i)=x(3,nc1(i))  
          x2(i)=x(1,nc2(i))  
          y2(i)=x(2,nc2(i))  
          z2(i)=x(3,nc2(i))  
          x3(i)=x(1,nc3(i))  
          y3(i)=x(2,nc3(i))  
          z3(i)=x(3,nc3(i))  
          x4(i)=x(1,nc4(i))  
          y4(i)=x(2,nc4(i))  
          z4(i)=x(3,nc4(i))    
        ENDDO      
      ENDIF     
C-----------
C     REPERE CONVECTE (ITERATIONS).
C-----------
      DO i=1,nel
       IF (off(i) <= one ) cycle 
       x43 = x4(i) - x3(i)
       y43 = y4(i) - y3(i)
       z43 = z4(i) - z3(i)
       x41 = x4(i) - x1(i)
       y41 = y4(i) - y1(i)
       z41 = z4(i) - z1(i)
       x42 = x4(i) - x2(i)
       y42 = y4(i) - y2(i)
       z42 = z4(i) - z2(i)
C
       rx(i) =  -x41
       ry(i) =  -y41
       rz(i) =  -z41
       sx(i) =  -x42
       sy(i) =  -y42
       sz(i) =  -z42
C
       tx(i) =  -x43
       ty(i) =  -y43
       tz(i) =  -z43
      END DO 
      CALL smortho3(
     1   off,     rx,      ry,      rz,
     2   sx,      sy,      sz,      tx,
     3   ty,      tz,      r11,     r12,
     4   r13,     r21,     r22,     r23,
     5   r31,     r32,     r33,     nel)
     
C-----------
C     PASSAGE AU REPERE CONVECTE.
C-----------
C       CALL S4RROTA3(OFF,
C     .   R11, R12, R13, R21, R22, R23, R31, R32, R33,
C     .   X0(I,1), X0(I,2), X0(I,3), X0(I,4), 
C     .   Y0(I,1), Y0(I,2), Y0(I,3), Y0(I,4), 
C     .   Z0(I,1), Z0(I,2), Z0(I,3), Z0(I,4))
C       CALL S4RROTA3(OFF,
C     .   R11, R12, R13, R21, R22, R23, R31, R32, R33,
C     .   VX0(I,1), VX0(I,2), VX0(I,3), VX0(I,4), 
C     .   VY0(I,1), VY0(I,2), VY0(I,3), VY0(I,4), 
C     .   VZ0(I,1), VZ0(I,2), VZ0(I,3), VZ0(I,4))
C    
C-----------
      RETURN