cpxpy3.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	cpxpy3 (jft, jlt, pm, sti, smstr, px1, px2, py1, py2, stir, ixc, area, x2, x3, x4, y2, y3, y4, z2, thk0, mat, nel)

Function/Subroutine Documentation

cpxpy3()

subroutine cpxpy3	(	integer	jft,
		integer	jlt,
			pm,
			sti,
			smstr,
			px1,
			px2,
			py1,
			py2,
			stir,
		integer, dimension(nixc,*)	ixc,
			area,
			x2,
			x3,
			x4,
			y2,
			y3,
			y4,
			z2,
			thk0,
		integer, dimension(mvsiz)	mat,
		integer	nel )

Definition at line 31 of file cpxpy3.F.

      use element_mod , only : nixc
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 MAT(MVSIZ), IXC(NIXC,*), JFT, JLT,NEL
C     REAL
      my_real
     .   pm(npropm,*), sti(*), stir(*),smstr(*),
     .   px1(*), px2(*), py1(*), py2(*),
     .   x2(mvsiz), x3(mvsiz), x4(mvsiz),area(mvsiz),
     .   y2(mvsiz), y3(mvsiz), y4(mvsiz), z2(mvsiz),thk0(mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, II(4), MX
C     REAL
      my_real
     .     young , 
     .     pxx, al2
C=======================================================================
      mx = mat(jft)  
      young=pm(20,mx)
C
      DO i=1,4
        ii(i) = nel*(i-1)
      ENDDO
C
      DO i=jft,jlt
        z2(i)  = zero
        px1(i) = smstr(ii(1)+i)
        px2(i) = smstr(ii(2)+i)
        py1(i) = smstr(ii(3)+i)
        py2(i) = smstr(ii(4)+i)
        x2(i)  =-two*py1(i)
        y2(i)  = zero
        x3(i)  =-two*py2(i)
        y3(i)  = two*px2(i)
        x4(i)  = zero
        y4(i)  =-two*px1(i)
        z2(i)  = zero
        area(i)= two*(py2(i)*px1(i)-py1(i)*px2(i))
      ENDDO
C
      DO i=jft,jlt
        pxx = max(px1(i)**2+py1(i)**2,px2(i)**2+py2(i)**2)
        sti(i) = young * pxx  * thk0(i) / area(i)
        stir(i)= zero
      ENDDO
C
      DO i=jft,jlt
        IF (ixc(4,i) == ixc(5,i)) THEN
          al2= (py1(i)+py2(i))**2
          al2= min(al2, px1(i)**2+py1(i)**2 )
          al2= four * min(al2, px2(i)**2+py2(i)**2 )
          sti(i) = half * area(i) * thk0(i) * young / al2
        ENDIF
      ENDDO
C-----------
      RETURN