lecm85_void.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	lecm85_void (iin, iout, uparam, maxuparam, nuparam, nuvar, ifunc, maxfunc, nfunc, parmat)

Function/Subroutine Documentation

lecm85_void()

subroutine lecm85_void	(	integer	iin,
		integer	iout,
			uparam,
		integer	maxuparam,
		integer	nuparam,
		integer	nuvar,
		integer, dimension(maxfunc)	ifunc,
		integer	maxfunc,
		integer	nfunc,
			parmat )

Definition at line 26 of file lecm85_void.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C----------+---------+---+---+--------------------------------------------
C VAR      | SIZE    |TYP| RW| DEFINITION
C----------+---------+---+---+--------------------------------------------
C IIN      |  1      | I | R | INPUT FILE UNIT (D00 file) 
C IOUT     |  1      | I | R | OUTPUT FILE UNIT (L00 file)
C UPARAM   | NUPARAM | F | W | USER MATERIAL PARAMETER ARRAY
C MAXUPARAM|  1      | I | R | MAXIMUM SIZE OF UPARAM 
C NUPARAM  |  1      | I | W | SIZE OF UPARAM =< MAXUPARAM
C NUVAR    |  1      | I | W | NUMBER OF USER ELEMENT VARIABLES
C----------+---------+---+---+--------------------------------------------
C IFUNC    | NFUNC   | I | W | FUNCTION NUMBER ARRAY
C MAXFUNC  |  1      | I | R | MAXIMUM SIZE OF IFUNC
C NFUNC    |  1      | I | W | SIZE OF IFUNC =< MAXFUNC
C----------+---------+---+---+--------------------------------------------
C (STIFINT)|                 | Replaced by PARMAT(1)
C PARMAT   |  3      | F | W | 1)STIFNESS MODULUS FOR INTERFACE
C          |(in v31) |   |   |       (AND INITIAL TIME STEP COMPUTATION)
C          |         |   |   | 2)YOUNG MODULUS (only needed for shells)
C          |         |   |   | 3)POISSON'S RATIO (only needed for shells)
C          |         |   |   | 4... Possible future use
C----------+---------+---+---+--------------------------------------------
C
C     This subroutine read the user material parameters.
C     The material cards that are common for all materials
C     (card 1 to 7 in version 2.2) have previously been read.
C     The NUPARAM material datas have to bee stored in UPARAM array.
C     If some standard radioss functions (time function or 
C     x,y function) are needed, this NFUNC function numbers have to 
C     bee stored in IFUNC array.  
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IIN,IOUT,MAXUPARAM,NUPARAM,NUVAR,MAXFUNC,NFUNC,
     .        IFUNC(MAXFUNC)
      my_real   uparam(maxuparam),parmat(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real e,nu,a11,a12,a44
C=======================================================================
C
C     EXAMPLE : ELASTIC LAW 
C
C=======================================================================
C
      READ(iin,'(2f16.0)')E,NU
      A11 = E * (ONE - NU) / (ONE + NU) / (ONE - TWO*NU)
      A12 = E * NU / (ONE + NU) / (ONE - TWO*NU)
      A44 = E / TWO / (ONE+NU)
C
.OR.      IF(NU<ZERONU>=HALF)THEN
         WRITE(IOUT,*)' ** error : wrong nu VALUE'
      ENDIF
      NUPARAM = 7
      IF(NUPARAM>MAXUPARAM)THEN
        WRITE(IOUT,*)' ** error : nuparam gt maxuparam'
        WRITE(IOUT,*)'      nuparam =',NUPARAM,
     .                ' maxuparam =',MAXUPARAM
      ELSE
C used in sigeps29 (solid 2d,3d)
        UPARAM(1) = A11 
        UPARAM(2) = A12
        UPARAM(3) = A44
        UPARAM(4) = E/(ONE-NU*NU)
        UPARAM(5) = NU*E/(ONE - NU*NU)
        UPARAM(6) = HALF*E/(ONE + NU)
 
        UPARAM(7) = E
      ENDIF
C
      NUVAR = 0
      NFUNC = 0
C used for interface (solid+shell)
      PARMAT(1) = A11
C used for shell
      PARMAT(2) = E
      PARMAT(3) = NU
C
      WRITE(IOUT,1000)
      WRITE(IOUT,1100)E,NU
C
 1000 FORMAT(
     & 5X,'  elastic user law ',/,
     & 5X,'  ---------------- ',//)
 1100 FORMAT(
     & 5X,'e . . . . . . . . . . . . . . . . . . .=',E12.4/
     & 5X,'nu. . . . . . . . . . . . . . . . . . .=',E12.4//)
C
      RETURN