preplyxfem.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/.
!||====================================================================
!||    preplyxfem   ../starter/source/properties/composite_options/stack/preplyxfem.F
!||--- called by ------------------------------------------------------
!||    lectur       ../starter/source/starter/lectur.F
!||====================================================================
      SUBROUTINE preplyxfem(
     1      MS_PLY0,ZI_PLY0,IEL, INOD , IXC ,MS_PLY,ZI_PLY,ADDCNE ,
     2      MSZ20,MSZ2)
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      "com04_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER  IXC(NIXC,*), ADDCNE(0:*),IEL(*),INOD(*)
      my_real 
     .        ms_ply0(numnod,*),  zi_ply0(numnod,*),
     .        ms_ply(nplyxfe,*),  zi_ply(nplyxfe,*),
     .        msz20(*),msz2(*) 
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,II,K,N1,J,N
      my_real
     .  msz
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
         addcne(0)=0
         addcne(1)= 0
         DO i=1,numnod
            ii = inod(i)
            addcne(i + 1) = 0
            msz = zero
            IF( ii > 0 )THEN
              msz2(ii) = zero
              DO j=1,nplymax
                 ms_ply(ii,j) = ms_ply0(i,j)
                 zi_ply(ii,j) = zi_ply0(i,j)
C  ementary compute
cc                 MSZ2(II) = MSZ20(I)  
c nodale compute                 
                 msz2(ii) = msz2(ii) + 
     .                      ms_ply(ii,j)*zi_ply(ii,j)*zi_ply(ii,j)
              ENDDO
               IF(msz2(ii) == zero) msz2(ii) = ep30
            ENDIF
         ENDDO    
C
         DO  i=1,numelc
            IF(iel(i) > 0 ) THEN
              DO  k=2,5
                n = inod(ixc(k,i)) + 1
                addcne(n) = addcne(n) + 1
              ENDDO
            ENDIF
          ENDDO 
          addcne(1) = 1
          DO i=2,nplyxfe + 1
             addcne(i) = addcne(i) + addcne(i-1)
          END DO        
C          
          RETURN
          END
!||====================================================================
!||    fillcne_pxfem   ../starter/source/properties/composite_options/stack/preplyxfem.F
!||--- called by ------------------------------------------------------
!||    lectur          ../starter/source/starter/lectur.F
!||--- calls      -----------------------------------------------------
!||====================================================================
      SUBROUTINE fillcne_pxfem(
     1   IEL , INOD, IXC,CEP, ADDCNE, CNE,CEL   )
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      "com01_c.inc"
#include      "com04_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IXC(NIXC,*),ADDCNE(0:*), CNE(*),IEL(*),INOD(*),
     .        CEP(*),CEL(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J, K, N, ADSKY(0:NPLYXFE+1),N1,PROC,NG,
     .        INDX, INDEX(NUMELC),P,II,NIN,INDEXC(2*NUMELC),
     .        WORK(70000),ITRI(NUMELC)
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
C   CALCUL DE CNE ADDCNE 
C-----------------------------------------------
      DO i = 0, nplyxfe + 1
        adsky(i) = addcne(i)
      ENDDO
C
C
C tri des elements locaux suivants num user
C
      DO i = 1, numelc
        itri(i) = ixc(7,i)
      ENDDO
      CALL my_orders(0,work,itri,indexc,numelc,1)
      
      indx = 0
      DO j=1,numelc
        i = indexc(j)
!!      DO   I = 1,NUMELC
       IF(iel(i) > 0 )THEN
         indx = indx + 1
         DO k=1,4
           n = ixc(k+1,i)
           n1 = inod(n)
           cne(adsky(n1)) = i
           adsky(n1) = adsky(n1) + 1
         ENDDO
         index(indx) = i
       ENDIF 
      ENDDO
C     
      DO proc = 1, nspmd
        nin = 0  
!!        DO I=1,INDX
          DO ii=1,numelc
!!           II = INDEX(I)
           ng = ii +  numels + numelq
           p = cep(ng) + 1 
           IF (p == proc) THEN
            j = iel(ii) 
            IF(j > 0) THEN 
              nin = nin + 1
              cel(j) = nin
            ENDIF  
           ENDIF
        ENDDO 
      ENDDO  
C      
      RETURN
      END