boltst.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	boltst (ip, bpreld, sig, tt, nel, npt, nsensor, sensor_tab, fun_id, sens_id)

Function/Subroutine Documentation

boltst()

subroutine boltst	(	integer, intent(in)	ip,
			bpreld,
			sig,
			tt,
		integer, intent(in)	nel,
		integer, intent(in)	npt,
		integer, intent(in)	nsensor,
		type (sensor_str_), dimension(nsensor), intent(in)	sensor_tab,
		integer, intent(in)	fun_id,
		integer, intent(in)	sens_id )

Definition at line 36 of file boltst.F.

C-----------------------------------------------  
C   M o d u l e s
C-----------------------------------------------  
      USE sensor_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      "com01_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NPT,NSENSOR,NEL,IP,FUN_ID,SENS_ID
      my_real
     .   tt
      my_real
     .   bpreld(nel,*),sig(nel,6)
      TYPE (SENSOR_STR_) ,DIMENSION(NSENSOR) ,INTENT(IN) :: SENSOR_TAB
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,ISENS 
      my_real
     .   ts,te,sfac
C-----------------------------------------------
!!! see preload_solid_ini.F90
      
       DO i=1,nel
        ts = bpreld(i,1)
        te = bpreld(i,2)
        IF(bpreld(i,3) == zero)THEN
          isens = sens_id
        
          IF (isens > 0) THEN
           ts = ts+sensor_tab(isens)%TSTART
           te = te+sensor_tab(isens)%TSTART
          ENDIF
        ENDIF
 
        IF(bpreld(i,3) < one)THEN
          IF(tt >= ts )THEN
           sfac = bpreld(i,10)
           sig(i,1) = sfac*bpreld(i,4) !SIG(I,1)+BPRELD(I,4)
           sig(i,2) = sfac*bpreld(i,5) !SIG(I,2)+BPRELD(I,5)
           sig(i,3) = sfac*bpreld(i,6) !SIG(I,3)+BPRELD(I,6)
           sig(i,4) = sfac*bpreld(i,7) !SIG(I,4)+BPRELD(I,7)
           sig(i,5) = sfac*bpreld(i,8) !SIG(I,5)+BPRELD(I,8)
           sig(i,6) = sfac*bpreld(i,9) !SIG(I,6)+BPRELD(I,9)
           bpreld(i,3) = real(ip)/real(npt)
           
!          Added to shift all the preloading duration in case of sensor use        
           bpreld(i,1) = ts
           bpreld(i,2) = te
          ENDIF
        ENDIF
       ENDDO
      
! Added to compensate loss of stress when filling the gap
       IF(bpreld(1,3) == one)THEN
         ts = bpreld(1,1)
         te = bpreld(1,2)
         IF(fun_id == 0)THEN
           IF (tt <= (ts+em01*(te-ts))) THEN
            DO i=1,nel
              sig(i,1) = bpreld(i,4)
              sig(i,2) = bpreld(i,5)
              sig(i,3) = bpreld(i,6)
              sig(i,4) = bpreld(i,7)
              sig(i,5) = bpreld(i,8)
              sig(i,6) = bpreld(i,9)
            ENDDO
           ELSE IF (tt >= te) THEN
            DO i=1,nel
              bpreld(i,3) = two
            ENDDO
           ENDIF
         ELSE ! ramping inital stress need more time to reach the pick
           IF (tt < te) THEN
             sfac = bpreld(1,10)
             sig(1:nel,1:6) = sfac*bpreld(1:nel,4:9)
           ELSE 
             bpreld(1:nel,3) = two
           ENDIF
         ENDIF
       END IF !(BPRELD(1,3) == ONE)THEN
 
       
      RETURN