cmain3pinch.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

cmain3pinch (elbuf_str, jft, jlt, nft, iparg, nel, mtn, ipla, ithk, group_param, pm, geo, npf, tf, bufmat, ssp, rho, viscmx, dt1c, sigy, area, exx, eyy, exy, exz, eyz, kxx, kyy, kxy, nu, off, thk0, mat, pid, for, mom, gstr, failwave, fwave_el, thk, eint, iofc, g, a11, a12, vol0, indxdel, ngl, zcfac, shf, gs, epsp, kfts, jhbe, alpe, dir_a, dir_b, igeo, ipm, ifailure, npg, tempel, die, jthe, iexpan, tempel0, ishplyxfem, ply_exx, ply_eyy, ply_exy, ply_exz, ply_eyz, ply_f, del_ply, th_iply, sig_iply, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, ng, table, ixfem, offi, a11_iply, elcrkini, dir1_crk, dir2_crk, aldt, ismstr, ir, is, nlay, npt, ixlay, ixel, isubstack, stack, f_def, itask, drape, varnl, pinch_local, forp, momp, ezzavg, areapinch)

Function/Subroutine Documentation

cmain3pinch()

subroutine cmain3pinch	(	type(elbuf_struct_), target	elbuf_str,
		integer	jft,
		integer	jlt,
		integer	nft,
		integer, dimension(*)	iparg,
		integer	nel,
		integer	mtn,
		integer	ipla,
		integer	ithk,
		type (group_param_)	group_param,
			pm,
			geo,
		integer, dimension(*)	npf,
			tf,
			bufmat,
			ssp,
			rho,
			viscmx,
			dt1c,
			sigy,
			area,
			exx,
			eyy,
			exy,
			exz,
			eyz,
			kxx,
			kyy,
			kxy,
			nu,
			off,
			thk0,
		integer, dimension(mvsiz)	mat,
		integer, dimension(mvsiz)	pid,
			for,
			mom,
			gstr,
		type (failwave_str_), target	failwave,
		integer, dimension(nel)	fwave_el,
			thk,
			eint,
		integer	iofc,
			g,
			a11,
			a12,
			vol0,
		integer, dimension(mvsiz)	indxdel,
		integer, dimension(mvsiz)	ngl,
			zcfac,
			shf,
			gs,
			epsp,
		integer	kfts,
		integer	jhbe,
			alpe,
			dir_a,
			dir_b,
		integer, dimension(npropgi,*)	igeo,
		integer, dimension(npropmi,*)	ipm,
		integer	ifailure,
		integer	npg,
			tempel,
			die,
		integer	jthe,
		integer	iexpan,
			tempel0,
		integer	ishplyxfem,
			ply_exx,
			ply_eyy,
			ply_exy,
			ply_exz,
			ply_eyz,
			ply_f,
			del_ply,
			th_iply,
			sig_iply,
			e1x,
			e1y,
			e1z,
			e2x,
			e2y,
			e2z,
			e3x,
			e3y,
			e3z,
		integer	ng,
		type(ttable), dimension(*)	table,
		integer	ixfem,
			offi,
		integer, dimension(mvsiz,*)	a11_iply,
		integer, dimension(*)	elcrkini,
			dir1_crk,
			dir2_crk,
			aldt,
		integer	ismstr,
		integer	ir,
		integer	is,
		integer	nlay,
		integer	npt,
		integer	ixlay,
		integer	ixel,
		integer	isubstack,
		type (stack_ply)	stack,
			f_def,
		integer	itask,
		type (drape_)	drape,
			varnl,
		type (pinch_local_struct_)	pinch_local,
			forp,
			momp,
			ezzavg,
			areapinch )

Definition at line 40 of file cmain3pinch.F.

C-----------------------------------------------
C   M o d u l e s
C----------------------------------------------- 
      USE table_mod
      USE elbufdef_mod
      USE stack_mod
      USE failwave_mod
      USE group_param_mod 
      USE pinchtype_mod
      USE pinch_local_mod
      USE drape_mod
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,NFT,IR,IS,NPT,NG,NPG,NLAY,NEL,MTN,IPLA,IOFC,ITHK,
     .   KFTS,JHBE,IFAILURE,JTHE,IEXPAN,ISHPLYXFEM,ISMSTR,
     .   IXLAY,IXEL,IXFEM,ISUBSTACK,ITASK
      INTEGER MAT(MVSIZ), PID(MVSIZ), INDXDEL(MVSIZ), NGL(MVSIZ),NPF(*), 
     .   IGEO(NPROPGI,*),IPM(NPROPMI,*),IPARG(*),
     .   A11_IPLY(MVSIZ,*),ELCRKINI(*),FWAVE_EL(NEL)
      my_real
     .   pm(npropm,*), geo(npropg,*),tf(*),thk(*),eint(jlt,*),for(nel,5), 
     .   mom(nel,3),epsp(*),gstr(nel,8),bufmat(*),g(*),a11(*),a12(*),vol0(*),
     .   shf(*),gs(*),sigy(mvsiz),rho(mvsiz),ssp(mvsiz),viscmx(mvsiz),
     .   off(mvsiz) ,thk0(mvsiz) ,nu(mvsiz),dt1c(mvsiz),
     .   exx(mvsiz), eyy(mvsiz), exy(mvsiz), exz(mvsiz), eyz(mvsiz),
     .   kxx(mvsiz), kyy(mvsiz), kxy(mvsiz), area(mvsiz),
     .   zcfac(mvsiz,2),alpe(mvsiz),dir_a(*),dir_b(*),tempel(*),die(*),
     .   tempel0(mvsiz), ply_f(mvsiz,5,*),
     .   ply_exx(mvsiz,*),ply_eyy(mvsiz,*),ply_exy(mvsiz,*), 
     .   ply_exz(mvsiz,*),ply_eyz(mvsiz,*),del_ply(mvsiz,3,*) ,
     .   th_iply(mvsiz,*), sig_iply(mvsiz,3,*),
     .   e1x(mvsiz), e1y(mvsiz), e1z(mvsiz),
     .   e2x(mvsiz), e2y(mvsiz), e2z(mvsiz),
     .   e3x(mvsiz), e3y(mvsiz), e3z(mvsiz),offi(mvsiz,*),
     .   dir1_crk(*),dir2_crk(*),dmg_scale(mvsiz),
     .   f_def(mvsiz,*),aldt(mvsiz),varnl(nel,*),
     .   forp(mvsiz),momp(mvsiz,2),ezzavg(mvsiz),areapinch(mvsiz)
      TYPE(TTABLE) TABLE(*)
      TYPE(ELBUF_STRUCT_), TARGET :: ELBUF_STR
      TYPE (STACK_PLY) :: STACK
      TYPE (FAILWAVE_STR_) ,TARGET :: FAILWAVE 
      TYPE (GROUP_PARAM_) :: GROUP_PARAM
      TYPE (PINCH_LOCAL_STRUCT_) :: PINCH_LOCAL
      TYPE (DRAPE_) :: DRAPE
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J,IT,MX,IGTYP,IFUNC_ALPHA,IPMAT_IPLY,ILAY,LAYNPT_MAX,LAY_MAX
      parameter(laynpt_max = 10)
      parameter(lay_max = 100)
      INTEGER MATLY(MVSIZ*LAY_MAX),MAT_IPLY(MVSIZ,NPT)
      my_real
     .   thkly(mvsiz*lay_max*laynpt_max),posly(mvsiz,lay_max*laynpt_max),
     .   thk_ly(nel,lay_max*laynpt_max),tstar(mvsiz),dtemp(mvsiz),
     .   sxx(mvsiz),syy(mvsiz), sxy(mvsiz),syz(mvsiz),szx(mvsiz),
     .   epsxx(mvsiz),epsyy(mvsiz),epsxy(mvsiz),epsyz(mvsiz),epszx(mvsiz),
     .   depsxx(mvsiz),depsyy(mvsiz),depsxy(mvsiz),depsyz(mvsiz),
     .   depszx(mvsiz), epspxx(mvsiz),epspyy(mvsiz),epspxy(mvsiz),
     .   epspyz(mvsiz), epspzx(mvsiz),etimp(mvsiz),tens(nel,5),eth(mvsiz)
      my_real alpha,df,deintth,t0,tm,fscal_alpha,dm
      INTEGER, DIMENSION(:) ,POINTER :: FLD_IDX
C-------------------------------------
      TYPE(BUF_LAY_) ,POINTER :: BUFLY
      TYPE(BUF_FAIL_),POINTER :: FBUF
      TYPE(L_BUFEL_) ,POINTER :: LBUF     
      TYPE(G_BUFEL_) ,POINTER :: GBUF     
C-----------------------------------------------
       my_real finter 
      EXTERNAL finter
C=======================================================================
      gbuf => elbuf_str%GBUF
      igtyp=igeo(11,pid(1))
      mx = mat(jft) 
      dm = group_param%VISC_DM   ! membrane damping
C--------------------------------------------------
      die(jft:jlt)   = zero
 
      ilay = 1
      it = 1
      lbuf => elbuf_str%BUFLY(ilay)%LBUF(ir,is,it)
C
      IF(mtn == 1 .OR. mtn == 91) THEN
        CALL mulawglcpinch(elbuf_str,
     1         jft     ,jlt     ,pm      ,for     ,mom     ,thk     ,
     2         eint    ,off     ,gstr    ,lbuf%PLA,dir_a   ,shf     ,
     3         mat     ,area    ,exx     ,eyy     ,exy     ,nel     ,
     4         exz     ,eyz     ,kxx     ,kyy     ,kxy     ,dm      ,
     5         pid     ,tf      ,npf     ,mtn     ,dt1c    ,a11     ,
     6         bufmat  ,ssp     ,rho     ,viscmx  ,iofc    ,a12     ,
     7         indxdel ,ngl     ,zcfac   ,gs      ,sigy    ,g       ,
     8         thk0    ,epsp    ,ipla    ,igeo    ,ipm     ,table   ,
     9         ir      ,is      ,f_def   ,ismstr  ,nu      ,vol0    ,
     a         kfts    ,pinch_local%EPINCHZZ,pinch_local%EPINCHXZ   ,
     b         pinch_local%EPINCHYZ      ,forp    ,momp    ,aldt    ,
     c         ezzavg  ,areapinch)
      ENDIF 
 
      RETURN