qfint2.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    qfint2   ../engine/source/elements/solid_2d/quad/qfint2.F
!||--- called by ------------------------------------------------------
!||    bforc2   ../engine/source/ale/bimat/bforc2.F
!||    qforc2   ../engine/source/elements/solid_2d/quad/qforc2.F
!||====================================================================
      SUBROUTINE qfint2(
     1   SIG,     PY1,     PY2,     PZ1,
     2   PZ2,     AREA,    VOL,     QVIS,
     3   F11,     F12,     F21,     F22,
     4   AX1,     AX2,     R22,     R23,
     5   R32,     R33,     NEL,     JCVT,
     6   SVIS)
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, INTENT(IN) :: JCVT
      INTEGER NEL
      my_real
     .   SIG(NEL,6),
     .   PY1(*), PY2(*), PZ1(*), PZ2(*), AREA(*), VOL(*), QVIS(*),
     .   F11(*), F12(*), F21(*), F22(*), AX1(*),AX2(*),
     .   r22(*), r23(*), r32(*), r33(*)
      my_real, DIMENSION(MVSIZ,6), INTENT(INOUT) :: svis
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     .   FAC(MVSIZ), S1(MVSIZ), S2(MVSIZ), S4(MVSIZ),
     .   T1,T2,T3,T4
C-----------------------------------------------
      IF (N2D==1) then
      DO i=1,nel
        fac(i)=area(i)*area(i)/ max(em15,vol(i))/four
      ENDDO
      ELSE
      DO i=1,nel
        fac(i)=zero
      ENDDO
      ENDIF
 
      IF(jcvt==0)THEN
       DO i=1,nel
         ax1(i)=(sig(i,3)+svis(i,3)-sig(i,1)-svis(i,1))*fac(i)
         ax2(i)=(sig(i,4)+svis(i,4))*fac(i)
       END DO
      ELSE
       DO i=1,nel
         s1(i)=sig(i,1)+svis(i,1)
         s2(i)=sig(i,2)+svis(i,2)
         s4(i)=sig(i,4)+svis(i,4)
         t1=s1(i)*r22(i)+s4(i)*r23(i)
         t2=s4(i)*r32(i)+s2(i)*r33(i)
         t3=s1(i)*r32(i)+s4(i)*r33(i)
         t4=s4(i)*r22(i)+s2(i)*r23(i)
         s1(i)=r22(i)*t1+r23(i)*t4
         s2(i)=r32(i)*t3+r33(i)*t2
         s4(i)=r22(i)*t3+r23(i)*t2
       END DO
       DO i=1,nel
         ax1(i)= (sig(i,3)+svis(i,3)-s1(i))*fac(i)
         ax2(i)= s4(i)*fac(i)
       END DO
      END IF
 
      DO i=1,nel
        s1(i) = sig(i,1)+svis(i,1)-qvis(i)
        s2(i) = sig(i,2)+svis(i,2)-qvis(i)
        s4(i) = sig(i,4)+svis(i,4)
        f11(i)=s1(i)*py1(i)+s4(i)*pz1(i)
        f21(i)=s2(i)*pz1(i)+s4(i)*py1(i)
        f12(i)=s1(i)*py2(i)+s4(i)*pz2(i)
        f22(i)=s2(i)*pz2(i)+s4(i)*py2(i)
      ENDDO
 
      RETURN
      END