set_elgroup_param.F File Reference

#include "implicit_f.inc"
#include "param_c.inc"
#include "com04_c.inc"
#include "tabsiz_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	set_elgroup_param (group_param_tab, iparg, ngroup, n2d, ipm, igeo, pm, geo, bufmat)

Function/Subroutine Documentation

set_elgroup_param()

subroutine set_elgroup_param	(	type(group_param_), dimension(ngroup)	group_param_tab,
		integer, dimension(nparg,ngroup)	iparg,
		integer	ngroup,
		integer	n2d,
		integer, dimension(npropmi,nummat)	ipm,
		integer, dimension(npropgi,numgeo)	igeo,
		dimension(npropm,nummat)	pm,
		dimension(npropg,numgeo)	geo,
		dimension(sbufmat)	bufmat )

Definition at line 29 of file set_elgroup_param.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE group_param_mod            
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "param_c.inc"
#include      "com04_c.inc"
#include      "tabsiz_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NGROUP,N2D
      INTEGER ,DIMENSION(NPARG,NGROUP)   :: IPARG
      INTEGER ,DIMENSION(NPROPMI,NUMMAT) :: IPM
      INTEGER ,DIMENSION(NPROPGI,NUMGEO) :: IGEO
      my_real ,DIMENSION(NPROPM,NUMMAT)  :: pm
      my_real ,DIMENSION(NPROPG,NUMGEO)  :: geo
      my_real ,DIMENSION(SBUFMAT)        :: bufmat
      TYPE(GROUP_PARAM_) ,DIMENSION(NGROUP) :: GROUP_PARAM_TAB
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NG,ITY,IGTYP,IDEF,IGMAT,IHBE,ISH3N,IEL,PID,IMAT,ILAW,IADD,
     .        NIP,ISMSTR,ICOAT
      my_real dm,tcoat
C=======================================================================
      DO ng=1,ngroup
        ity   = iparg(5,ng)
        IF (n2d == 0 .and. (ity == 3 .or. ity == 7)) THEN
          ilaw  = iparg(1,ng)
          imat  = iparg(18,ng)
          pid   = iparg(62,ng)
          iel   = iparg(3,ng)+1
          igtyp = iparg(38,ng)
          dm    = geo(16,pid)
          idef  = igeo(31,pid)
          igmat = igeo(98,pid)
c          IHBE  = IGEO(10,PID)   ! should be the same, GEO(171) is obsolete, check !!!
          ihbe  = nint(geo(171,pid))
          ish3n = igeo(18,pid)
c
c         Set default membrane shell viscosity parameter (DM)         
c
          IF (dm == zero) THEN
            SELECT CASE (ilaw)
              CASE (65)
                dm = fiveem2
              CASE (19, 119)
                IF (ity == 3 .AND. ihbe == 23) THEN
                  dm = zep015
                ELSE
                  dm = fourth
                END IF
              CASE (58, 158)
                IF (ihbe == 11) THEN
                  dm = one_fifth
                ELSE
                  dm = three_over_4 * em01
                END IF
              CASE (25, 125, 127)
                IF (ihbe == 23) dm = em01
              CASE DEFAULT
                IF (ihbe == 23 .AND. (igtyp == 1 .or. igtyp == 9)) dm = zep015
            END SELECT
          END IF
c------   specific composite IGMAT > 0
          IF (igmat > 0) THEN  ! .and. (ILAW /= 15 .AND. ILAW /= 27)) THEN
            IF (dm == zero) dm = fiveem2
          ELSEIF ( ilaw == 15 .or. ilaw == 27 .or. ilaw == 25 .or. ilaw==125 
     .                                                        .or. ilaw == 127) THEN
            IF (dm == zero) dm = fiveem2
          ENDIF !(IGMAT > 0
c          
          group_param_tab(ng)%VISC_DM = dm  
c
c------   special settings for law119
c          
          IF (ilaw == 119) THEN
            iadd = max(1,ipm(7,imat))-1
            icoat = nint(bufmat(iadd+14))
            tcoat = bufmat(iadd+16)
            ismstr = 11   ! force small strain formulation for law 119 regardless of input
            ihbe   = 11   ! force full shell integration (Batoz) for law 119 regardless of input
            IF (icoat > zero .AND. tcoat > zero) THEN 
              nip = 3
              ! update shell thickness in the property to take coating into account
              geo(1,pid) = geo(1,pid) + two*tcoat
            ELSE
              nip = 1
            END IF
            iparg(6,ng) = nip
            iparg(9,ng) = ismstr
            iparg(23,ng)= ihbe
          END IF
           
 
        ELSE
c
          group_param_tab(ng)%VISC_DM = zero  
c
        ENDIF
c
        group_param_tab(ng)%VISC_DN = zero  
      
      ENDDO
c-----------
      RETURN