c3coork3.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	c3coork3 (jft, jlt, x, ixtg, offg, geo, pid, mat, ngl, area, irep, npt, ismstr, nlay, elbuf_str, smstr, dir_a, dir_b, xl2, xl3, yl3, r11, r12, r13, r21, r22, r23, r31, r32, r33, k11, k12, k13, k22, k23, k33, m11, m12, m13, m22, m23, m33, mf11, mf12, mf13, mf22, mf23, mf33, fm12, fm13, fm23, off, nel)

Function/Subroutine Documentation

c3coork3()

subroutine c3coork3	(	integer	jft,
		integer	jlt,
			x,
		integer, dimension(nixtg,*)	ixtg,
			offg,
			geo,
		integer, dimension(*)	pid,
		integer, dimension(*)	mat,
		integer, dimension(*)	ngl,
			area,
		integer	irep,
		integer	npt,
		integer	ismstr,
		integer	nlay,
		type(elbuf_struct_)	elbuf_str,
			smstr,
			dir_a,
			dir_b,
			xl2,
			xl3,
			yl3,
			r11,
			r12,
			r13,
			r21,
			r22,
			r23,
			r31,
			r32,
			r33,
			k11,
			k12,
			k13,
			k22,
			k23,
			k33,
			m11,
			m12,
			m13,
			m22,
			m23,
			m33,
			mf11,
			mf12,
			mf13,
			mf22,
			mf23,
			mf33,
			fm12,
			fm13,
			fm23,
			off,
		integer	nel )

Definition at line 34 of file c3coork3.F.

C-----------------------------------------------
C   M o d u l e s
C----------------------------------------------- 
      USE elbufdef_mod
      use element_mod , only : nixtg
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 ,IREP,NPT,ISMSTR,NLAY,NEL
      INTEGER IXTG(NIXTG,*),PID(*),MAT(*),NGL(*)
C     REAL
      my_real
     .   x(3,*),offg(*),area(*),
     .   geo(npropg,*),xl2(*),xl3(*),yl3(*),
     .   r11(*),r12(*),r13(*),
     .   r21(*),r22(*),r23(*),
     .   r31(*),r32(*),r33(*),
     .    k11(9,*),k12(9,*),k13(9,*),
     .    k22(9,*),k23(9,*),k33(9,*),
     .    m11(9,*),m12(9,*),m13(9,*),
     .    m22(9,*),m23(9,*),m33(9,*),
     .    mf11(9,*),mf12(9,*),mf13(9,*),
     .    mf22(9,*),mf23(9,*),mf33(9,*),
     .    fm12(9,*),fm13(9,*),fm23(9,*),off(*),smstr(*),
     .    dir_a(nel,*),dir_b(nel,*)
      TYPE(ELBUF_STRUCT_) :: ELBUF_STR
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NC1, NC2, NC3,I,II(3),J,K, MAT_1
C     REAL
      my_real
     .     x1(mvsiz), y1(mvsiz), z1(mvsiz),deta1(mvsiz),
     .     x2(mvsiz), y2(mvsiz), z2(mvsiz),
     .     x3(mvsiz), y3(mvsiz), z3(mvsiz),
     .     rx(mvsiz), ry(mvsiz), rz(mvsiz),
     .     sx(mvsiz), sy(mvsiz), sz(mvsiz)
C-----------------------------------------------
      DO i=1,3
        ii(i) = nel*(i-1)
      ENDDO
C
      mat_1 = ixtg(1,jft)
      DO i=jft,jlt
       mat(i)  = mat_1
       nc1     = ixtg(2,i)
       nc2     = ixtg(3,i)
       nc3     = ixtg(4,i)
       pid(i)  = ixtg(5,i)
       ngl(i)  = ixtg(6,i)
C----------------------------
C     COORDONNEES
C----------------------------
       x1(i)=x(1,nc1)
       y1(i)=x(2,nc1)
       z1(i)=x(3,nc1)
       x2(i)=x(1,nc2)
       y2(i)=x(2,nc2)
       z2(i)=x(3,nc2)
       x3(i)=x(1,nc3)
       y3(i)=x(2,nc3)
       z3(i)=x(3,nc3)
      ENDDO
C
      DO i=jft,jlt
        rx(i)=x2(i)-x1(i)
        sx(i)=x3(i)-x1(i)
        ry(i)=y2(i)-y1(i)
        sy(i)=y3(i)-y1(i)
        rz(i)=z2(i)-z1(i)
        sz(i)=z3(i)-z1(i)
      ENDDO 
C----------------------------
C     LOCAL SYSTEM
C----------------------------
      k = 2
      CALL clskew3(jft,jlt,k,
     .   rx, ry, rz, 
     .   sx, sy, sz, 
     .   r11,r12,r13,r21,r22,r23,r31,r32,r33,deta1,offg )
C----------------------------
C     CHANGE TO LOCAL SYSTEM
C----------------------------
      DO i=jft,jlt
       area(i)=half*deta1(i)
       xl2(i)=r11(i)*rx(i)+r21(i)*ry(i)+r31(i)*rz(i)
       xl3(i)=r11(i)*sx(i)+r21(i)*sy(i)+r31(i)*sz(i)
       yl3(i)=r12(i)*sx(i)+r22(i)*sy(i)+r32(i)*sz(i)
      ENDDO
C----------------------------
C     SMALL STRAIN
C----------------------------
      IF(ismstr==1.OR.ismstr==2)THEN
        DO i=jft,jlt
          IF(abs(offg(i))==2.)THEN
            xl2(i)=smstr(ii(1)+i)
            xl3(i)=smstr(ii(2)+i)
            yl3(i)=smstr(ii(3)+i)
          ELSE
            smstr(ii(1)+i)=xl2(i)
            smstr(ii(2)+i)=xl3(i)
            smstr(ii(3)+i)=yl3(i)
          ENDIF
        ENDDO
      ENDIF
      IF(ismstr==1)THEN
        DO i=jft,jlt
          IF(offg(i) == one) offg(i)=two
        ENDDO
      ENDIF
C----------------------------
C     Orthotropy later
C----------------------------
      IF (irep > 0) THEN
       CALL cortdir3(elbuf_str,dir_a,dir_b ,jft    ,jlt    ,
     .               nlay   ,irep   ,rx    ,ry     ,rz     , 
     .               sx     ,sy     ,sz    ,r11    ,r21    ,
     .               r31    ,r12    ,r22   ,r32    ,nel    )
      ENDIF
c
      DO j=1,9
       DO i=jft,jlt
          k11(j,i)=zero
          k12(j,i)=zero
          k13(j,i)=zero
          k22(j,i)=zero
          k23(j,i)=zero
          k33(j,i)=zero
          m11(j,i)=zero
          m12(j,i)=zero
          m13(j,i)=zero
          m22(j,i)=zero
          m23(j,i)=zero
          m33(j,i)=zero
       ENDDO 
      ENDDO 
      DO i=jft,jlt
       off(i)=offg(i)
      ENDDO
C
      RETURN