cm27in3.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/.
!||====================================================================
!||    cm27in3                ../starter/source/materials/mat/mat027/cm27in3.F
!||--- called by ------------------------------------------------------
!||    c3init3                ../starter/source/elements/sh3n/coque3n/c3init3.F
!||    cinit3                 ../starter/source/elements/shell/coque/cinit3.F
!||    cmaini3                ../starter/source/elements/sh3n/coquedk/cmaini3.F
!||--- calls      -----------------------------------------------------
!||    ancmsg                 ../starter/source/output/message/message.F
!||    fretitl2               ../starter/source/starter/freform.F
!||--- uses       -----------------------------------------------------
!||    message_mod            ../starter/share/message_module/message_mod.F
!||====================================================================
      SUBROUTINE cm27in3(ELBUF_STR,
     .                   GEO ,IGEO,PM ,IPM ,IX ,NIX,
     .                   NLAY,IR  ,IS ,IMAT )
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE elbufdef_mod            
      USE message_mod
      USE names_and_titles_mod , ONLY : nchartitle
C-----------------------------------------------
C INITIALISE LES DIRECTIONS DE FISSURES
C INITIALISE LES EPAISSEURS ET LES MATERIAUX DES COUCHES
C
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      "vect01_c.inc"
#include      "param_c.inc"
#include      "scr17_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IGEO(NPROPGI,*),IX(NIX,*),NIX,IPM(NPROPMI,*),NLAY,IR,IS,IMAT
      my_real GEO(NPROPG,*),PM(NPROPM,*)
      TYPE(elbuf_struct_), TARGET :: ELBUF_STR
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IPTHK,IPMAT,IPPOS,MID,MLN,IMID,IPID,PID
      INTEGER I,N,I1,I2,I3,J,IGTYP,II,L_DMG,ILAYER,IT,NPTT
      my_real thkl,pos0,dp
      CHARACTER(LEN=NCHARTITLE) :: TITR,TITR1
      my_real, DIMENSION(:), POINTER :: dir_dmg
C
      TYPE(buf_lay_) ,POINTER :: BUFLY
      TYPE(L_BUFEL_) ,POINTER :: LBUF
C======================================================================|
      IF(npt==0)THEN
        imid=ix(1,1)
        ipid = ix(nix-1,1)
        pid = igeo(1,ipid)
        mid = ipm(1,imid)
        mln = nint(pm(19,imid))
        CALL fretitl2(titr,igeo(npropgi-ltitr+1,pid),ltitr)
        CALL fretitl2(titr1,ipm(npropmi-ltitr+1,imid),ltitr)
        CALL ancmsg(msgid=23, anmode=aninfo, msgtype=msgerror, i1=pid, c1=titr, i2=mid, c2=titr1, i3=27)
      ENDIF
C
      igtyp=nint(geo(12,imat))
      IF (igtyp /= 51 .AND. igtyp /= 52) THEN
        DO n=1,npt
          ilayer = n
          IF (nlay > 1) THEN
            lbuf => elbuf_str%BUFLY(ilayer)%LBUF(ir,is,1)
          ELSE
            lbuf => elbuf_str%BUFLY(1)%LBUF(ir,is,ilayer)
          ENDIF
C
          l_dmg = elbuf_str%BUFLY(1)%L_DMG
          dir_dmg => lbuf%DMG(1:l_dmg*llt)
C
          DO i=lft,llt
            dir_dmg(i) = one
            dir_dmg(i+llt) = zero
          ENDDO
        ENDDO
      ELSEIF (igtyp == 51) THEN
        DO ilayer=1,nlay
          nptt = elbuf_str%BUFLY(ilayer)%NPTT
          DO it=1,nptt
            lbuf => elbuf_str%BUFLY(ilayer)%LBUF(ir,is,it)
C
            l_dmg = elbuf_str%BUFLY(ilayer)%L_DMG
            dir_dmg => lbuf%DMG(1:l_dmg*llt)
C
            DO i=lft,llt
              dir_dmg(i) = one
              dir_dmg(i+llt) = zero
            ENDDO
          ENDDO
        ENDDO
      ENDIF
C
      ipthk = 300
      ippos = 400
      ipmat = 100
      IF(igtyp/=11 .AND. igtyp/=17 . and. igtyp/=51) THEN
        thkl = one / npt
        pos0 =-half*(one + thkl)
        DO n=1,npt                                                
          i1=ippos+n                                              
          i2=ipthk+n                                              
          i3=ipmat+n                                              
          dp =  n*thkl
          DO i=lft,llt                                            
            geo(i1,imat) = pos0 + dp 
            geo(i2,imat) = thkl                              
          ENDDO                                                   
        ENDDO                                                     
      ENDIF
C-----------
      RETURN
      END