01.05.2017 Views

4569846498

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

478 Multibody Systems Approach to Vechicle Dynamics<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

C<br />

CS TPROP(4)<br />

CA TPROP(5)<br />

CR TPROP(7)<br />

DZ TPROP(8)<br />

U TPROP(11)<br />

Convert sign on alpha<br />

RALPHA CPROP(2)<br />

RGAMMA CPROP (3)<br />

CG TPROP (6)<br />

ALPHA -RALPHA<br />

GAMMA RGAMMA<br />

Initialize force values<br />

FX 0.D0<br />

FY 0.D0<br />

FZ 0.D0<br />

TY 0.D0<br />

TZ 0.D0<br />

IF(DEFL .LE. 0.D0) THEN<br />

GOTO 1000<br />

ENDIF<br />

Calculate normal loads due to stiffness (always .LE. zero)<br />

FZDEFL -DEFL*CZ<br />

Calculate normal loads due to damping<br />

FZDAMP -2.D0*SQRT(AMASS*CZ)*DZ*(DEFLD)<br />

Calculate total normal force (fz)<br />

FZ MIN (0.0D0, (FZDEFLFZDAMP) )<br />

Convert to kN and change sign<br />

FZP -FZ<br />

Compute longitudinal force<br />

IF(ABS(SLIP) .LE. ABS(SSTAR)) THEN<br />

FX -CS*SLIP<br />

ELSE<br />

FX1 U*ABS(FZ)<br />

FX2 (U*FZ)**2/(4.D0*ABS(SLIP)*CS)<br />

FX -(FX1-FX2)*SIGN(1.0D0,SLIP)<br />

ENDIF

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!