ccoef3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	ccoef3 (jft, jlt, pm, mat, geo, pid, off, area, sti, stir, shf, thk0, thk02, nu, g, ym, a11, a12, thk, ssp, rho, h1, h2, h3, vol0, vol00, alpe, gs, mtn, ithk, ismstr, npt, kfts, srh1, srh2, srh3, igeo, a11r, isubstack, pm_stack)

Function/Subroutine Documentation

ccoef3()

subroutine ccoef3	(	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,
			h1,
			h2,
			h3,
			vol0,
			vol00,
			alpe,
			gs,
		integer	mtn,
		integer	ithk,
		integer	ismstr,
		integer	npt,
		integer	kfts,
			srh1,
			srh2,
			srh3,
		integer, dimension(npropgi,*)	igeo,
			a11r,
		integer	isubstack,
			pm_stack )

Definition at line 29 of file ccoef3.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"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT, JLT, MTN,ITHK,
     .            ISMSTR,NPT ,KFTS ,ISUBSTACK
      INTEGER MAT(MVSIZ), PID(MVSIZ), IGEO(NPROPGI,*)
C     REAL
      my_real
     .   geo(npropg,*), pm(npropm,*), off(*), area(*),
     .   sti(*),stir(*),shf(*),thk0(*),thk02(*),thk(*),
     .   nu(*),g(*),ym(*),a11(*),a12(*),gs(*),
     .   h1(*),h2(*),h3(*),vol0(*),ssp(*),rho(*),alpe(*),vol00(*),
     .   srh1(*)    ,srh2(*)   ,srh3(*),a11r(*),pm_stack(20,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,ISH,MX,IPID,J,IGTYP,IPGMAT,IGMAT
C     REAL
      my_real
     .   fac1(mvsiz),fsh, fac1tmp,
     .   rho_1, ym_1, nu_1, g_1, a11_1, a12_1, ssp_1      
C=======================================================================
      IF(ithk>0.AND.ismstr/=3)THEN
        DO i=jft,jlt
          alpe(i)=one
C --- A CORRIGER : VOL00(I) = THK0(I)*AREA0(I)
          vol00(i) = thk0(i)*area(i)
          thk0(i)=thk(i)
          vol0(i) = thk0(i)*area(i)
          thk02(i) = thk0(i)*thk0(i)
        ENDDO
      ELSE
        DO i=jft,jlt
          alpe(i)=one
C --- A CORRIGER : VOL00(I) = THK0(I)*AREA0(I)
          vol00(i) = thk0(i)*area(i)
          vol0(i) = vol00(i)
          thk02(i) = thk0(i)*thk0(i)
        ENDDO
      ENDIF
C
      igtyp = igeo(11,pid(1))
      igmat = igeo(98,pid(1))
      ipgmat = 700
 
      IF(igtyp == 11 .AND. igmat > 0) THEN
        mx = mat(jft)
        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)
         rho_1=pm(1,mx)
         ym_1  =pm(20,mx)
         nu_1  =pm(21,mx)
         g_1   =pm(22,mx)
         a11_1 =pm(24,mx)
         a12_1 =pm(25,mx)
         ssp_1 =pm(27,mx)
         DO i=jft,jlt
           rho(i)=rho_1
           ym(i)  =ym_1
           nu(i)  =nu_1
           g(i)   =g_1
           a11(i) =a11_1
           a12(i) =a12_1
           ssp(i) =ssp_1
         ENDDO
      ENDIF !igtyp 
C
      DO i=jft,jlt
       ipid=pid(i)
       h1(i) =geo(13,ipid)
       h2(i) =geo(14,ipid)
       h3(i) =geo(15,ipid)
       srh1(i)=geo(18,ipid)
       srh2(i)=geo(19,ipid)
       srh3(i)=geo(20,ipid)
      ENDDO
C    
      IF(npt==1) THEN
        DO i=jft,jlt
         shf(i)=zero
        ENDDO
      ELSE
        DO i=jft,jlt
C ds021  41i 26/9/00  +++
          fac1tmp  =2.*(1.+nu(i))*thk02(i)
          ish   =nint(geo(37,pid(i)))
          fsh   =geo(38,pid(i))
          shf(i)=fsh*(1. - ish +
     .      ish*fac1tmp / (fsh*area(i)+fac1tmp) )
          fac1(i) = fac1tmp
C ds021  41i 26/9/00  ---
        ENDDO
      ENDIF
      DO i=jft,jlt
        gs(i)=g(i)*shf(i)
      ENDDO
C
      RETURN