4569846498

486 Multibody Systems Approach to Vechicle Dynamics
FYSCFACT*FY
ENDIF
ENDIF
C Aligning Torque based on final FY including camber force
C Not quite right as Camber force has no pneumatic trail.
C Pay attention if using this for motorcycle work;
C consider reworking so that TZ does not include camber force.
C Aligning torque form is a bit poor anyway.
C -- Contact Patch Length
CP_LEN(R1**2 - (R1FZ/CZ)**2)**0.5*2.0
IF(ABS(ALPHA) .GT. 1.D-10) THEN
IF(ABS(ALPHA) .LE. ALPHA_C) THEN
TZ-FY*CP_LEN*(1-ABS(ALPHA)/ALPHA_C)/4
ELSE
TZ0.0
ENDIF
ENDIF
C
C
Other Calculations as FIALA Tyre:
Rolling resistance moment (TY)
CALL STEP(ABS(WSPIN), .25D0*WSPNMX, ZERO, WSPNMX,
* ONE, IORD, SCFACT, ERFLG)
IF (WSPIN .GE. 0.0) THEN
TY-SCFACT*CR*FZ
ELSE
TYSCFACT*CR*FZ
ENDIF
C Compute righting moment due to lateral Contact Patch Shift (TX)
C Use CG as “shape factor” to add to or subtract righting moment
C from MSC.ADAMS’ Toroidal assumption. CG 1 fatter than toroid
C CG 1 more like blade
C
C TXFZ*0.5*SIN(2*CAMBER)*R2/2*(CG - 1)
C Copy the calculated values for FX, FY, FZ, TY & TZ to FSAE
C and TSAE arrays and return. Scale them to start with to
C let the simulation wake up gently. But not FZ, otherwise
C we fall through the floor.
1000 CALL STEP(TIME, ZERO, ZERO, ONE, ONE, IORD, TSCALE,
ERFLG)
FSAE(1) FX*TSCALE
FSAE(2) FY*TSCALE
FSAE(3) FZ
TSAE(1) TX*TSCALE
TSAE(2) TY*TSCALE
TSAE(3) TZ*TSCALE
FPROP(1) 0.0
FPROP(2) 0.0
C
RETURN
END