IF ((IBASE.EQ.4) .AND. ( DSNOW .LE. (1.0/XCM2IN))) RRES1 = 0.0
C
C ..... Traction calculation
C
IF( NVEH(I) .EQ. WHEEL )THEN
AREA1 = GCA(I,JPSI)
XNWORT = NWHL(I)
XWID1 = SECTW(I)
ELSE
AREA1 = GCA(I,JPSI) / 2.0
XNWORT = 2.0
XWID1 = TRAKWD(I)
END IF
C
Weight on a single traction element
WGHT1 = WGHT(I) / XNWORT
PRES1 = WGHT1/AREA1
C
IF( IBASE .NE. 4 )THEN
C
Snow on road (IBASE=1)
C
Snow on frozen or strong (RCI>100) soil (IBASE=2)
C
Snow on unfrozen soil (RCI<=100) (IBASE=3)
IF (SIGMA0.GE.0.55) THEN
C
hard packed snow
TRCN1 = 0.321 * (PRES1*X2KPA)**0.97
ELSE
TRCN1 = 0.851 * (PRES1*X2KPA)**0.823
END IF
ELSE
C
Snow on ice base (IBASE=4)
IF ( DSNOW .LE. (1.0/XCM2IN)) THEN
TRCN1 = 0.1 * (PRES1*X2KPA)
ELSE
IF (SIGMA0.GE.0.55) THEN
C
Disturbed snow
TRCN1 = 0.321 * (PRES1*X2KPA)**0.97
ELSE
C
Undisturbed snow
TRCN1 = 0.127 * (PRES1*X2KPA)**1.06
END IF
END IF
END IF
C
C
TRCN1 is a stress at this point
TRCN1 = (TRCN1/X2KPA) * AREA1
C
TFOW(I) = TRCN1 * XNWORT / WGHT(I)
RTOW = RRES1 * XNWORT / WGHT(I) + RHARD(I)
C
RTOWPB(I)= RTOW
RTOWT(I) = RTOW
DOWPB(I) = TFOW(I) - RTOW
C
20 CONTINUE
C
46
Previous Page
