add_mass_stat.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"
#include "units_c.inc"
#include "scr18_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	add_mass_stat (ms, in, stifn, stifr, itab, totmas)

Function/Subroutine Documentation

add_mass_stat()

subroutine add_mass_stat	(		ms,
			in,
			stifn,
			stifr,
		integer, dimension(*)	itab,
			totmas )

Definition at line 31 of file add_mass_stat.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE my_alloc_mod
      USE plot_curve_mod
      USE message_mod
C-----------------------------------------------
C     TRI DES DT NODAUX ET IMPRESSIONS
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER ITAB(*)
      my_real
     .   ms(numnod),in(numnod),stifn(numnod),stifr(numnod),totmas
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
#include      "units_c.inc"
#include      "scr18_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,N,COMPT,K,NVAL
      INTEGER :: IERROR
      INTEGER, DIMENSION(:), ALLOCATABLE :: PERM
      my_real, DIMENSION(:), ALLOCATABLE,TARGET ::
     .   dt
      my_real, DIMENSION(:),POINTER ::
     .   tmp
      my_real dtsca,target_dt_tab(20),per_adm_tab(20),target_dt_tab2(20),summ,sumk
C=======================================================================
C
C--------------------------------------------------------------------------------------
C     DUPLICATED FROM OUTRIN - STATS ON ADDED MASS + TARGET NODAL TIME STEP ESTIMATION
C--------------------------------------------------------------------------------------
C
      ALLOCATE(dt(2*numnod))
      CALL my_alloc(perm,numnod)
      mass0_start = totmas
      dtsca = zep9
C
      tmp => dt(numnod+1:numnod*2)
C 
      DO i=1,numnod
        IF((ms(i)/=zero).AND.(stifn(i)>em20))THEN
          dt(i)=ms(i)/stifn(i)
        ELSE
          dt(i)=ep30       ! -1- free nodes dt=1e30 instead of dt=1416=sqrt(1e6).
                           ! -2- nodal time step from VOID elem (stifn<em20) is dt=EP30 too.
        ENDIF
      ENDDO
C
C  --- Sorting
C
      summ = zero
      sumk = zero
 
      DO i=1,numnod
        tmp(i)=i
        perm(i) = i
        IF (dt(i) < ep30) THEN
          summ = summ + ms(i)
          sumk = sumk + stifn(i)
        ENDIF
      ENDDO
C
      CALL myqsort(numnod,dt,perm,ierror)
      tmp(1:numnod) = perm(1:numnod)
 
      DEALLOCATE( perm )
C
      nval = 20
      DO i=1,20
        per_adm_tab(i) = 0.005*i
      ENDDO
C
C----- Step 1 - determination target time step for targeted added mass (from 1% to 10%) - dtsca = 0.9
C
      dtsca = zep9
      CALL find_dt_target(ms,stifn,target_dt_tab,per_adm_tab,dt,tmp,dtsca,totmas,nval,numnod)
C
C----- Step 2 - determination target time step for targeted added mass (from 1% to 10%) - dtsca = 0.67
C
      dtsca = zep67
      CALL find_dt_target(ms,stifn,target_dt_tab2,per_adm_tab,dt,tmp,dtsca,totmas,nval,numnod)
C
C  -------------------------------------------------------------------------------------------------------------------------
C
      IF (dt(1) < ep20) THEN       
C
C----- Visual output of added mass estimation - only if nodal time step is significant
C
        WRITE(iout,2001)
        WRITE(iout,2002)
        CALL plot_curve(target_dt_tab, per_adm_tab, 20, input_size_x=60, input_size_y=24,
     .                  input_txt_x="TARGET TIME STEP",input_txt_y="DM/M - at first cycle - scaling factor = 0.9",input_symbol="*")
C
        WRITE(iout,3000) '            DM/M at first cycle        ','Target time step - scaling factor = 0.9',
     .                   'Target time step - scaling factor = 0.67'
        WRITE(iout,3000) '  -------------------------------------','---------------------------------------',
     .                   '---------------------------------------'
C
        DO i=1,5
          k = 2*i
          WRITE(iout,3001) per_adm_tab(k),target_dt_tab(k),target_dt_tab2(k)
        ENDDO
C
        WRITE(iout,4001)
        WRITE(iout,4002)
C
      ENDIF
 
C-----------
      DEALLOCATE( dt)
 
C-----------
 2001 FORMAT(//,'     INITIAL ADDED MASS ESTIMATION for /DT/NODA/CST')
 2002 FORMAT(   '     ----------------------------------------------',//)
C
 3000 FORMAT(2x,a,2x,'|',2x,a,2x,'|',2x,a)
 3001 FORMAT(17x,1pg9.3,17x,'|',17x,1pg9.3,17x,'|',17x,1pg9.3)
C
 4001 FORMAT(//,'    WARNING : These added mass estimations are valid at first cycle only.')
 4002 FORMAT('    Because of element deformation or contacts the added mass can increase during the run',//)
C----------
 
      RETURN