cm1inif.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	cm1inif (jft, jlt, for, mom, ithk, thk, eint, nu, g, a1, a2, gs, kxx, kyy, kxy, exx, eyy, exy, eyz, exz, area)

Function/Subroutine Documentation

cm1inif()

subroutine cm1inif	(	integer	jft,
		integer	jlt,
			for,
			mom,
		integer	ithk,
			thk,
			eint,
			nu,
			g,
			a1,
			a2,
			gs,
			kxx,
			kyy,
			kxy,
			exx,
			eyy,
			exy,
			eyz,
			exz,
			area )

Definition at line 29 of file cm1inif.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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT, JLT, ITHK
      my_real
     .    exx(mvsiz) , eyy(mvsiz) , exy(mvsiz) , 
     .    exz(mvsiz) , eyz(mvsiz) , kxx(mvsiz) ,
     .    kyy(mvsiz) , kxy(mvsiz) , area(mvsiz)
      my_real
     .   for(5,*) ,mom(3,*), thk(*), eint(2,*), nu(*),
     .   g(*), a1(*), a2(*), gs(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real thk08, ezz,
     .    degmb(mvsiz), degfx(mvsiz), b1(mvsiz), b2(mvsiz), b3(mvsiz)
C=======================================================================
      DO i=jft,jlt
        thk08  =thk(i)*one_over_12
        b1(i)  =a1(i)*thk08
        b2(i)  =a2(i)*thk08
        b3(i)  =g(i) *thk08
      ENDDO
C
      DO i=jft,jlt
        for(1,i)=a1(i)*exx(i)+a2(i)*eyy(i)
        for(2,i)=a1(i)*eyy(i)+a2(i)*exx(i)
        for(3,i)=g(i) *exy(i)
        for(4,i)=gs(i)*eyz(i)
        for(5,i)=gs(i)*exz(i)
C
        mom(1,i)=b1(i)*kxx(i)+b2(i)*kyy(i)
        mom(2,i)=b1(i)*kyy(i)+b2(i)*kxx(i)
        mom(3,i)=b3(i)*kxy(i)
      ENDDO
C
      DO i=jft,jlt
        degmb(i) = for(1,i)*exx(i)+for(2,i)*eyy(i)+for(3,i)*exy(i)
        degfx(i) = mom(1,i)*kxx(i)+mom(2,i)*kyy(i)
     +                   +mom(3,i)*kxy(i)
        eint(1,i) = degmb(i)*half*thk(i)*area(i)
        eint(2,i) = degfx(i)*half*thk(i)*thk(i)*area(i)
      ENDDO
      IF(ithk==1) THEN
        DO i=jft,jlt
          ezz = -nu(i) * (exx(i) + eyy(i)) / (one-nu(i))
          thk(i) = thk(i) * (one + ezz)
        ENDDO
      ENDIF
C-----------
      RETURN