c3coef3.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"
#include "com04_c.inc"
#include "impl1_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	c3coef3 (jft, jlt, pm, mat, geo, pid, off, area, sti, stir, shf, thk0, thk02, nu, g, ym, a11, a12, thk, ssp, rho, vol0, gs, mtn, ithk, npt, ismstr, vol00, igeo, a11r, isubstack, pm_stack, nel, zoffset)

Function/Subroutine Documentation

c3coef3()

subroutine c3coef3	(	integer	jft,
		integer	jlt,
			pm,
		integer, dimension(mvsiz)	mat,
			geo,
		integer, dimension(mvsiz)	pid,
			off,
			area,
			sti,
			stir,
			shf,
			thk0,
			thk02,
			nu,
			g,
			ym,
			a11,
			a12,
			thk,
			ssp,
			rho,
			vol0,
			gs,
		integer	mtn,
		integer	ithk,
		integer	npt,
		integer	ismstr,
			vol00,
		integer, dimension(npropgi,numgeo)	igeo,
			a11r,
		integer	isubstack,
			pm_stack,
		integer, intent(in)	nel,
		intent(out)	zoffset )

Definition at line 29 of file c3coef3.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      "param_c.inc"
#include      "com04_c.inc"
#include      "impl1_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT, JLT,MTN,ITHK,NPT,ISMSTR,ISUBSTACK
      INTEGER MAT(MVSIZ), PID(MVSIZ), IGEO(NPROPGI,NUMGEO)
      INTEGER , INTENT(IN) :: NEL
C     REAL
      my_real
     .   geo(npropg,numgeo), pm(npropm,*), off(*), area(*),
     .   sti(*),stir(*),shf(*),thk0(*),thk02(*),thk(*),
     .   nu(*),g(*),ym(*),a11(*),a12(*),
     .   vol0(*),vol00(*),ssp(*),rho(*),gs(*),
     .   a11r(*),pm_stack(20,*)
      my_real , DIMENSION(NEL), INTENT(OUT) :: zoffset
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,ISH,MX,IPID,J,IGTYP,IGMAT,IPGMAT,IPOS
C     REAL
      my_real
     .   fac1,fsh,viscdef,z0  
C-----------------------------------------------
      IF (ithk>0.AND.ismstr/=3.AND.ismdisp == 0) THEN
        DO i=jft,jlt
          thk0(i)=thk(i)
        ENDDO
      ENDIF
C
      DO i=jft,jlt
        thk02(i) = thk0(i)*thk0(i)
        vol0(i) = thk0(i)*area(i)
C --- A CORRIGER : VOL00(I) = THK0(I)*AREA0(I)
        vol00(i) = thk0(i)*area(i)
      ENDDO
C
      IF (mtn == 19) THEN
        viscdef=fourth
      ELSEIF (mtn == 25.OR.mtn == 27) THEN
        viscdef=fiveem2
      ELSE
        viscdef=zero
      ENDIF
C
      igtyp = igeo(11,pid(1))
      igmat = igeo(98,pid(1))
      ipgmat = 700
      IF(igtyp == 11 .AND. igmat > 0) THEN
          DO i=jft,jlt
            mx = pid(i)
            rho(i) = geo(ipgmat +1 ,mx)   
            ym(i)  = geo(ipgmat +2 ,mx) 
            nu(i)  = geo(ipgmat +3 ,mx)
            g(i)   = geo(ipgmat +4 ,mx) 
            a11(i) = geo(ipgmat +5 ,mx)
            a12(i) = geo(ipgmat +6 ,mx)  
            a11r(i)= geo(ipgmat +7 ,mx)
            ssp(i) = geo(ipgmat +9 ,mx)
        ENDDO
      ELSEIF(igtyp == 52 .OR.
     .       ((igtyp == 17 .OR. igtyp == 51) .AND. igmat > 0)) THEN
          DO i=jft,jlt
            rho(i) = pm_stack(1 ,isubstack)   
            ym(i)  = pm_stack(2 ,isubstack) 
            nu(i)  = pm_stack(3 ,isubstack)
            g(i)   = pm_stack(4 ,isubstack) 
            a11(i) = pm_stack(5 ,isubstack)
            a12(i) = pm_stack(6 ,isubstack)  
            a11r(i)= pm_stack(7 ,isubstack)
            ssp(i) = pm_stack(9 ,isubstack)
        ENDDO 
      ELSE
      mx  =mat(jft)
        DO i=jft,jlt
             rho(i)=pm(1,mx)
             ym(i) =pm(20,mx)
             nu(i) =pm(21,mx)
             g(i)  =pm(22,mx)
             a11(i)=pm(24,mx)
             a12(i)=pm(25,mx)
             ssp(i)=pm(27,mx)
        ENDDO
      ENDIF
C
      IF (npt == 1) THEN
        DO i=jft,jlt
          shf(i) = zero
        ENDDO
      ELSE
        DO i=jft,jlt
          fac1 = two*(one+nu(i))*thk02(i)
          ish  = nint(geo(37,pid(i)))
          fsh  = geo(38,pid(i))
          shf(i)=fsh*(one - ish + ish*fac1 / (fsh*area(i)+fac1) )
        ENDDO
      ENDIF
      DO i=jft,jlt
         gs(i)=g(i)*shf(i)
      ENDDO
      z0 = geo(199,pid(1))
      zoffset(jft:jlt)  = zero
      SELECT CASE(igtyp)
        CASE (1,9,10,11,16)
           DO i=jft,jlt
            zoffset(i) = z0
           ENDDO 
        CASE (17,51,52)
            ipos   = igeo(99,pid(1))
             IF(ipos == 2) THEN
               DO i=jft,jlt
                 zoffset(i) = z0 - half*thk0(i)
               ENDDO 
              ELSEIF (ipos== 3 .OR. ipos == 4) THEN
               DO i=jft,jlt
                 z0= half*thk0(i)
                 zoffset(i) = z0
               ENDDO  
             ENDIF    
        CASE DEFAULT 
            zoffset(jft:jlt)  = zero
      END SELECT
C-----------
      RETURN