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)

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 )

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
      my_real
     .   tt, dt1
      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,K,ISENS 
      my_real
     .   ts,te
C-----------------------------------------------
C
!     Retrieve sysid of the sensors at cycle=0
      IF (ncycle == 0 .AND. ip == 1) THEN
       DO i=1,nel
        isens = int(bpreld(i,10))
        IF (isens > 0)THEN            
         DO k=1,nsensor
          IF(isens==sensor_tab(k)%SENS_ID) GOTO 100  ! do it in staerter !!!
         ENDDO
100      bpreld(i,10)=k
        ENDIF
       ENDDO
      ENDIF
      
       DO i=1,nel
        ts = bpreld(i,1)
        te = bpreld(i,2)
        IF(bpreld(i,3) == zero)THEN
          isens = int(bpreld(i,10))
        
          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
           sig(i,1) = bpreld(i,4) !SIG(I,1)+BPRELD(I,4)
           sig(i,2) = bpreld(i,5) !SIG(I,2)+BPRELD(I,5)
           sig(i,3) = bpreld(i,6) !SIG(I,3)+BPRELD(I,6)
           sig(i,4) = bpreld(i,7) !SIG(I,4)+BPRELD(I,7)
           sig(i,5) = bpreld(i,8) !SIG(I,5)+BPRELD(I,8)
           sig(i,6) = 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 (TT <= (TE+EM01*(TE-TS))) THEN       ! sb : seems to be an error 
         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
       ENDIF
 
       
      RETURN