cfint3.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "com01_c.inc"
#include "com08_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	cfint3 (jft, jlt, for, mom, thk0, thk02, px1, px2, py1, py2, area, z2, f11, f12, f13, f14, f21, f22, f23, f24, f31, f32, f33, f34, h11, h12, h13, h21, h22, h23, h31, h32, h33, b11, b12, b13, b14, b21, b22, b23, b24, nel, m11, m12, m13, m14, m21, m22, m23, m24, m31, m32, m33, m34, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, ihbe, npt, fzero)

Function/Subroutine Documentation

cfint3()

subroutine cfint3	(	integer	jft,
		integer	jlt,
			for,
			mom,
			thk0,
			thk02,
			px1,
			px2,
			py1,
			py2,
			area,
			z2,
			f11,
			f12,
			f13,
			f14,
			f21,
			f22,
			f23,
			f24,
			f31,
			f32,
			f33,
			f34,
			h11,
			h12,
			h13,
			h21,
			h22,
			h23,
			h31,
			h32,
			h33,
			b11,
			b12,
			b13,
			b14,
			b21,
			b22,
			b23,
			b24,
		integer	nel,
			m11,
			m12,
			m13,
			m14,
			m21,
			m22,
			m23,
			m24,
			m31,
			m32,
			m33,
			m34,
			e1x,
			e1y,
			e1z,
			e2x,
			e2y,
			e2z,
			e3x,
			e3y,
			e3z,
		integer	ihbe,
		integer	npt,
			fzero )

Definition at line 29 of file cfint3.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,NPT,IHBE,NEL
C     REAL
      my_real
     .   for(nel,5), mom(nel,3), thk0(*), thk02(*),
     .   area(*),   px1(*), px2(*), py1(*), py2(*), z2(*),
     .   f11(mvsiz), f12(mvsiz), f13(mvsiz), f14(mvsiz),
     .   f21(mvsiz), f22(mvsiz), f23(mvsiz), f24(mvsiz),
     .   f31(mvsiz), f32(mvsiz), f33(mvsiz), f34(mvsiz),
     .     b11(mvsiz), b12(mvsiz), b13(mvsiz), b14(mvsiz),
     .     b21(mvsiz), b22(mvsiz), b23(mvsiz), b24(mvsiz),
     .     h11(mvsiz), h12(mvsiz), h13(mvsiz), h14(mvsiz),
     .     h21(mvsiz), h22(mvsiz), h23(mvsiz), h24(mvsiz),
     .     h31(mvsiz), h32(mvsiz), h33(mvsiz), h34(mvsiz)
      my_real
     .     m11(mvsiz), m12(mvsiz), m13(mvsiz), m14(mvsiz),
     .     m21(mvsiz), m22(mvsiz), m23(mvsiz), m24(mvsiz),
     .     m31(mvsiz), m32(mvsiz), m33(mvsiz), m34(mvsiz)
      my_real
     .     e1x(mvsiz), e1y(mvsiz), e1z(mvsiz),
     .     e2x(mvsiz), e2y(mvsiz), e2z(mvsiz),
     .     e3x(mvsiz), e3y(mvsiz), e3z(mvsiz),fzero(3,4,*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com08_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J
C     REAL
      my_real
     .     m4(mvsiz), m5(mvsiz), f1(mvsiz), f2(mvsiz), f3(mvsiz), 
     .     f4(mvsiz), f5(mvsiz), m1(mvsiz), m2(mvsiz), m3(mvsiz),
     .     g11t(mvsiz), g21t(mvsiz), g12t(mvsiz), g22t(mvsiz), 
     .     g13t(mvsiz), g23t(mvsiz), g14t(mvsiz), g24t(mvsiz),
     .     g31t(mvsiz), g32t(mvsiz), g33t(mvsiz), g34t(mvsiz),
     .     f1a,f2a,f3a,f4a,f5a, g11ta, g21ta, g12ta, g22ta, 
     .     g13ta, g23ta, g14ta, g24ta, g31ta, g32ta, g33ta, g34ta,
     .     m1a,m2a,m3a,temp1,temp2,temp3
C=======================================================================
C               POUR ISIGI = 5 OU -5
C     CHANGEMENT DE REPERE DE FZERO POUR TT=DT1
C          REPERE GLOBAL ---> REPERE LOCAL
C-----------------------------------------------
      IF(iabs(isigi)==5.AND.ncycle==1.AND.irun==1) THEN
        DO i=jft,jlt
          DO j=1,4
            temp1=fzero(1,j,i)
            temp2=fzero(2,j,i)
            temp3=fzero(3,j,i)
            fzero(1,j,i)=temp1*e1x(i)+temp2*e1y(i)+temp3*e1z(i)
            fzero(2,j,i)=temp1*e2x(i)+temp2*e2y(i)+temp3*e2z(i)
            fzero(3,j,i)=temp1*e3x(i)+temp2*e3y(i)+temp3*e3z(i)
          ENDDO
        ENDDO
      ENDIF
C-----------------------------------------------
C               POUR ISIGI = 5 OU -5
C     PRISE EN COMPTE DES FORCES D EQUILIBRAGE
C      DANS L EVALUATION DES EFFORTS INTERNES
C-----------------------------------------------
      IF((iabs(isigi)==5.AND.ncycle>=1.AND.irun==1).OR.
     .     (iabs(isigi)==5.AND.irun>1)) THEN
          DO i=jft,jlt
            f1a= for(i,1)*thk0(i)
            f2a= for(i,2)*thk0(i)
            f3a= for(i,3)*thk0(i)
            f4a= for(i,4)*thk0(i)
            f5a= for(i,5)*thk0(i)
            m4(i)= f4a*area(i)
            m5(i)= f5a*area(i)
C
            f12(i)= f1a*px2(i)+f3a*py2(i)
            f22(i)= f2a*py2(i)+f3a*px2(i)
            f32(i)= f5a*px2(i)+f4a*py2(i)
            f11(i)= f1a*px1(i)+f3a*py1(i)
            f21(i)= f2a*py1(i)+f3a*px1(i)
            f31(i)= f5a*px1(i)+f4a*py1(i)
C
            g11ta=f11(i)+h11(i)
            g13ta=h13(i)-f11(i)
            g21ta=f21(i)+h21(i)
            g23ta=h23(i)-f21(i)
            g31ta=f31(i)+h31(i)
            g33ta=h33(i)-f31(i)
            g12ta=f12(i)+h12(i)
            g22ta=f22(i)+h22(i)
            g32ta=f32(i)+h32(i)
            g14ta=-g11ta-g12ta-g13ta
            g24ta=-g21ta-g22ta-g23ta
            g34ta=-g31ta-g32ta-g33ta
C
            g11ta=g11ta+fzero(1,1,i)
            g13ta=g13ta+fzero(1,3,i)
            g21ta=g21ta+fzero(2,1,i)
            g23ta=g23ta+fzero(2,3,i)
            g31ta=g31ta+fzero(3,1,i)
            g33ta=g33ta+fzero(3,3,i)
            g12ta=g12ta+fzero(1,2,i)
            g22ta=g22ta+fzero(2,2,i)
            g32ta=g32ta+fzero(3,2,i)
            g14ta=g14ta+fzero(1,4,i)
            g24ta=g24ta+fzero(2,4,i)
            g34ta=g34ta+fzero(3,4,i)
C
            f11(i)=e1x(i)*g11ta+e2x(i)*g21ta+e3x(i)*g31ta
            f12(i)=e1x(i)*g12ta+e2x(i)*g22ta+e3x(i)*g32ta
            f13(i)=e1x(i)*g13ta+e2x(i)*g23ta+e3x(i)*g33ta
            f14(i)=e1x(i)*g14ta+e2x(i)*g24ta+e3x(i)*g34ta
C
            f21(i)=e1y(i)*g11ta+e2y(i)*g21ta+e3y(i)*g31ta
            f22(i)=e1y(i)*g12ta+e2y(i)*g22ta+e3y(i)*g32ta
            f23(i)=e1y(i)*g13ta+e2y(i)*g23ta+e3y(i)*g33ta
            f24(i)=e1y(i)*g14ta+e2y(i)*g24ta+e3y(i)*g34ta
C
            f31(i)=e1z(i)*g11ta+e2z(i)*g21ta+e3z(i)*g31ta
            f32(i)=e1z(i)*g12ta+e2z(i)*g22ta+e3z(i)*g32ta
            f33(i)=e1z(i)*g13ta+e2z(i)*g23ta+e3z(i)*g33ta
            f34(i)=e1z(i)*g14ta+e2z(i)*g24ta+e3z(i)*g34ta
          ENDDO
      ENDIF
      IF(iabs(isigi)/=5.OR.tt==0) THEN
          DO i=jft,jlt
            f1a= for(i,1)*thk0(i)
            f2a= for(i,2)*thk0(i)
            f3a= for(i,3)*thk0(i)
            f4a= for(i,4)*thk0(i)
            f5a= for(i,5)*thk0(i)
            m4(i)= f4a*area(i)
            m5(i)= f5a*area(i)
C
            f12(i)= f1a*px2(i)+f3a*py2(i)
            f22(i)= f2a*py2(i)+f3a*px2(i)
            f32(i)= f5a*px2(i)+f4a*py2(i)
            f11(i)= f1a*px1(i)+f3a*py1(i)
            f21(i)= f2a*py1(i)+f3a*px1(i)
            f31(i)= f5a*px1(i)+f4a*py1(i)
C
            g11ta=f11(i)+h11(i)
            g13ta=h13(i)-f11(i)
            g21ta=f21(i)+h21(i)
            g23ta=h23(i)-f21(i)
            g31ta=f31(i)+h31(i)
            g33ta=h33(i)-f31(i)
            g12ta=f12(i)+h12(i)
            g22ta=f22(i)+h22(i)
            g32ta=f32(i)+h32(i)
            f11(i)=e1x(i)*g11ta+e2x(i)*g21ta+e3x(i)*g31ta
            f12(i)=e1x(i)*g12ta+e2x(i)*g22ta+e3x(i)*g32ta
            f13(i)=e1x(i)*g13ta+e2x(i)*g23ta+e3x(i)*g33ta
            f21(i)=e1y(i)*g11ta+e2y(i)*g21ta+e3y(i)*g31ta
            f22(i)=e1y(i)*g12ta+e2y(i)*g22ta+e3y(i)*g32ta
            f23(i)=e1y(i)*g13ta+e2y(i)*g23ta+e3y(i)*g33ta
            f31(i)=e1z(i)*g11ta+e2z(i)*g21ta+e3z(i)*g31ta
            f32(i)=e1z(i)*g12ta+e2z(i)*g22ta+e3z(i)*g32ta
            f33(i)=e1z(i)*g13ta+e2z(i)*g23ta+e3z(i)*g33ta
            f14(i)=-f11(i)-f12(i)-f13(i)
            f24(i)=-f21(i)-f22(i)-f23(i)
            f34(i)=-f31(i)-f32(i)-f33(i)
          ENDDO
      ENDIF
C
      IF((ihbe>=2).AND.iabs(npt)/=1)THEN
        DO i=jft,jlt
          m4(i)= m4(i)+(h21(i)+h23(i))*z2(i)
          m5(i)= m5(i)+(h11(i)+h13(i))*z2(i)
        ENDDO
C
      ENDIF
C
      DO i=jft,jlt
        m1a= mom(i,1)*thk02(i)
        m2a= mom(i,2)*thk02(i)
        m3a= mom(i,3)*thk02(i)
        m4(i)= m4(i)*fourth
        m5(i)= m5(i)*fourth
C
        m11(i)=-m2a*py1(i)-m3a*px1(i)
        m21(i)= m1a*px1(i)+m3a*py1(i)
        m12(i)=-m2a*py2(i)-m3a*px2(i)
        m22(i)= m1a*px2(i)+m3a*py2(i)
C
        g11ta = m11(i)-m4(i)+b11(i)
        g13ta =-m11(i)-m4(i)+b13(i)
        g12ta = m12(i)-m4(i)+b12(i)
        g14ta =-m12(i)-m4(i)+b14(i)
        g21ta = m21(i)+m5(i)+b21(i)
        g23ta =-m21(i)+m5(i)+b23(i)
        g22ta = m22(i)+m5(i)+b22(i)
        g24ta =-m22(i)+m5(i)+b24(i)
C
        m11(i)=e1x(i)*g11ta+e2x(i)*g21ta
        m12(i)=e1x(i)*g12ta+e2x(i)*g22ta
        m13(i)=e1x(i)*g13ta+e2x(i)*g23ta
        m14(i)=e1x(i)*g14ta+e2x(i)*g24ta
        m21(i)=e1y(i)*g11ta+e2y(i)*g21ta
        m22(i)=e1y(i)*g12ta+e2y(i)*g22ta
        m23(i)=e1y(i)*g13ta+e2y(i)*g23ta
        m24(i)=e1y(i)*g14ta+e2y(i)*g24ta
        m31(i)=e1z(i)*g11ta+e2z(i)*g21ta
        m32(i)=e1z(i)*g12ta+e2z(i)*g22ta
        m33(i)=e1z(i)*g13ta+e2z(i)*g23ta
        m34(i)=e1z(i)*g14ta+e2z(i)*g24ta
      ENDDO
c-----------
      RETURN