princ_u1.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	princ_u1 (nel, es1, es2, es3, es4, es5, es6, ev, dirprv)

Function/Subroutine Documentation

princ_u1()

subroutine princ_u1	(	integer	nel,
			es1,
			es2,
			es3,
			es4,
			es5,
			es6,
			ev,
			dirprv )

Definition at line 31 of file princ_u1.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 
C-----------------------------------------------
C----------------------------------------------------------------
C  I N P U T   A R G U M E N T S
C----------------------------------------------------------------
      INTEGER       NEL     
      my_real
     .   es1(mvsiz), es2(mvsiz), es3(mvsiz), 
     .   es4(mvsiz), es5(mvsiz), es6(mvsiz),
     .   ev(mvsiz,3),dirprv(mvsiz,3,3)
C----------------------------------------------------------------
C  L O C A L  V A R I B L E S
C----------------------------------------------------------------
      INTEGER    I,J,K
      my_real
     .            av(mvsiz,6),evv(mvsiz,3)
C----------------------------------------------------------------
C--------- [B]=[F][F]^t -1 strain-----
      DO i=1,nel
       av(i,1) = es1(i)
       av(i,2) = es2(i)
       av(i,3) = es3(i)
       av(i,4) = es4(i)
       av(i,5) = es5(i)
       av(i,6) = es6(i)
      ENDDO
      CALL valpvecop_v(av,evv,dirprv,nel)
 
      ev(1:nel,1:3)=sqrt(max(zero,(evv(1:nel,1:3)+ one))) - one
C----------------
C      
      RETURN