ede112.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/.
!||====================================================================
!||    ede112                 ../engine/source/ale/euler2d/ede112.f
!||--- called by ------------------------------------------------------
!||    qforc2                 ../engine/source/elements/solid_2d/quad/qforc2.F
!||--- uses       -----------------------------------------------------
!||    ale_connectivity_mod   ../common_source/modules/ale/ale_connectivity_mod.F
!||====================================================================
      SUBROUTINE ede112(PM  ,V   ,X   ,IXQ ,ALE_CONNECT ,
     .                  WYZ ,DYZ ,DZY  ,EYY ,EZZ )
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE ale_connectivity_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "param_c.inc"
#include      "com04_c.inc"
#include      "com08_c.inc"
#include      "vect01_c.inc"
#include      "tabsiz_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
! SPMD CASE : SIXQ >= NIXQ*NUMELQ    (SIXQ = NIXQ*NUMELQ_L+NIXQ*NQVOIS_L)
! IXQ(1:NIXQ, 1:NUMELQ) local elems
!    (1:NIXQ, NUMELQ+1:) additional elems (also on adjacent domains but connected to the boundary of the current domain)
      INTEGER,INTENT(IN) :: IXQ(NIXQ,SIXQ/NIXQ)
! SPMD CASE : SX >= 3*NUMNOD    (SX = 3*(NUMNOD_L+NRCVVOIS_L))
! X(1:3,1:NUMNOD) : local nodes
!  (1:3, NUMNOD+1:) additional nodes (also on adjacent domains but connected to the boundary of the current domain)
      my_real,INTENT(IN) :: x(3,sx/3)
!
      my_real,INTENT(IN) :: v(3,numnod), pm(npropm,nummat)
      my_real,INTENT(INOUT) :: eyy(*),ezz(*),dyz(*), dzy(*),wyz(*)
      TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, II,IADJ,ML,NC1,NC2,NC3,NC4,J,IAD2
      my_real vy13, vy24, vz13,vz24,py1,py2,pz1,pz2,aire
C-----------------------------------------------
C   S o u r c e   L i n e s 
C-----------------------------------------------
      DO i=lft,llt
       ii=i+nft
       iad2 = ale_connect%ee_connect%iad_connect(ii)
       DO j=1,4
         iadj=ale_connect%ee_connect%connected(iad2 + j - 1)
         IF(iadj /= 0) THEN
           ml=nint(pm(19,ixq(1,iadj)))
           IF(ml /= 11) THEN
             nc1=ixq(2,iadj)
             nc2=ixq(3,iadj)
             nc3=ixq(4,iadj)
             nc4=ixq(5,iadj)
             py1=half*(x(3,nc2)-x(3,nc4))
             py2=half*(x(3,nc3)-x(3,nc1))
             pz1=half*(x(2,nc4)-x(2,nc2))
             pz2=half*(x(2,nc1)-x(2,nc3))    
             aire=two*(py1*pz2-pz1*py2)
             vy13=v(2,nc1)-v(2,nc3)
             vy24=v(2,nc2)-v(2,nc4)
             vz13=v(3,nc1)-v(3,nc3)
             vz24=v(3,nc2)-v(3,nc4)
             eyy(i)=(py1*vy13+py2*vy24)/aire
             ezz(i)=(pz1*vz13+pz2*vz24)/aire
             dzy(i)=(py1*vz13+py2*vz24)/aire
             dyz(i)=(pz1*vy13+pz2*vy24)/aire
             ezz(i)=zero
             wyz(i)=half*dt1*(dzy(i)-dyz(i))
           ENDIF
         ENDIF
       enddo!next J
      enddo!next I
C-----------------------------------------------
      RETURN
      END