4569846498
Appendix B: Fortran tyre model subroutines 479 C C C C C C C C Compute lateral force CFYA0 SHFYA8*FZPA9A10*GAMMA DFY(A1*FZPA2)*(1-A15*GAMMA**2)*FZP IF(ALPHASHFY.LT.0.0)THEN DUMFY-1.0 ELSE DUMFY1.0 ENDIF EFY(A6*FZPA7)*(1-(A16*GAMMAA17)*DUMFY) BFY((A3*SIN(2*ATAN(FZP/A4)))*(1-A5*ABS(GAMMA)))/ (CFYDFY) SVFYA11*FZPA12(A13*FZP**2A14*FZP)*GAMMA PHIFY(1-EFY)*(ALPHASHFY)(EFY/BFY)* ATAN (BFY*(ALPHASHFY)) FYPDFY*SIN(CFY*ATAN(BFY*PHIFY))SVFY Change sign FYFYP Compute self aligning moment CTZC0 SHTZC11*FZPC12C13*GAMMA DTZ(C1*FZP**2C2*FZP)*(1-C18*GAMMA**2) IF(ALPHASHTZ.LT.0.0)THEN DUMTZ-1.0 ELSE DUMTZ1.0 ENDIF ETZ(C7*FZP**2C8*FZPC9)*(1-(C19*GAMMAC20)* DUMTZ) ETZETZ/(1-C10*ABS(GAMMA)) BTZ((C3*FZP**2C4*FZP)*(1-C6*ABS(GAMMA))* EXP(-C5*FZP))/(CTZDTZ) SVTZC14*FZPC15(C16*FZP**2C17*FZP)* GAMMA PHITZ(1-ETZ)*(ALPHASHTZ)(ETZ/BTZ)*ATAN(BTZ* (ALPHASHTZ)) TZPDTZ*SIN(CTZ*ATAN(BTZ*PHITZ))SVTZ C C Convert to Nmm and change sign C TZTZP*1000.0 C C Copy the calculated values for FX, FY, FZ, TY & TZ to FSAE C and TSAE arrays C 1000 FSAE(1)FX FSAE(2)FY
480 Multibody Systems Approach to Vechicle Dynamics C FSAE(3)FZ TSAE(1)0.0 TSAE(2)0.0 TSAE(3)TZ FPROP(1)0.0 FPROP(2)0.0 RETURN END B.3 The Harty tyre model subroutine C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C MDI TIRSUB : Prodrive Concept Tyre Model A Quick & Simplified Tyre Model which plugs in as the FIALA model does, with a “TIRE” statement. Unlike FIALA, critical slip angle is broadly independent of load and initial cornering stiffness is strongly load dependent. The model does handle comprehensive slip. Lateral force generation is zero at peak longitudinal force slip ratio (typically about 20%) but returns to a value around one tenth of the peak lateral force as the wheel progresses beyond that limit. This may result in poor post-spin performance. The force generated with locked wheels is currently aligned with the wheel plane; this is incorrect. Longitudinal force generation is assumed to be symmetric for tractive and braking slip. This is not generally true beyond the critical slip ratio for real tyres but is reasonable up to that point. This tyre will over estimate longitudinal forces for tractive slip and slightly underestimate them for braking slip in the post-critical regions. Camber thrust is included as for the motorcycle tire model using “taut string” logic. Lateral migration of the contact patch is not included, unlike a motorcycle tyre model. Aligning Torque calculation includes the lateral force due to camber. This is not quite right as the camber force mechanism has no pneumatic trail associated with it. Pay attention if using this for motorcycle work; consider reworking it so that TZ does not include the camber force. The form of the aligning torque is a bit poor and would benefit from some more thought; pneumatic trail collapses linearly with lateral force.
- Page 452 and 453: Simulation output and interpretatio
- Page 454 and 455: Simulation output and interpretatio
- Page 456 and 457: Simulation output and interpretatio
- Page 458 and 459: Simulation output and interpretatio
- Page 460 and 461: Simulation output and interpretatio
- Page 462 and 463: Simulation output and interpretatio
- Page 464 and 465: 8 Active systems 8.1 Introduction M
- Page 466 and 467: Active systems 443 mechanical actua
- Page 468 and 469: Active systems 445 Table 8.1 MSC.AD
- Page 470 and 471: Active systems 447 Body acceleratio
- Page 472 and 473: Active systems 449 impressions of t
- Page 474 and 475: Active systems 451 For this reason,
- Page 476 and 477: Appendix A: Vehicle model system sc
- Page 478 and 479: Appendix A: Vehicle model system sc
- Page 480 and 481: Appendix A: Vehicle model system sc
- Page 482 and 483: Appendix A: Vehicle model system sc
- Page 484 and 485: Appendix A: Vehicle model system sc
- Page 486 and 487: Appendix A: Vehicle model system sc
- Page 488 and 489: Appendix A: Vehicle model system sc
- Page 490 and 491: Appendix A: Vehicle model system sc
- Page 492 and 493: Appendix A: Vehicle model system sc
- Page 494 and 495: Appendix A: Vehicle model system sc
- Page 496 and 497: Appendix B: Fortran tyre model subr
- Page 498 and 499: Appendix B: Fortran tyre model subr
- Page 500 and 501: Appendix B: Fortran tyre model subr
- Page 504 and 505: Appendix B: Fortran tyre model subr
- Page 506 and 507: Appendix B: Fortran tyre model subr
- Page 508 and 509: Appendix B: Fortran tyre model subr
- Page 510 and 511: Appendix C: Glossary of terms Agili
- Page 512 and 513: Appendix C: Glossary of terms 489 a
- Page 514 and 515: Appendix C: Glossary of terms 491 t
- Page 516 and 517: Appendix C: Glossary of terms 493 e
- Page 518 and 519: Appendix C: Glossary of terms 495 m
- Page 520 and 521: Appendix C: Glossary of terms 497 h
- Page 522 and 523: Appendix C: Glossary of terms 499 t
- Page 524 and 525: Appendix C: Glossary of terms 501 s
- Page 526 and 527: References 503 Blundell, M.V. The m
- Page 528 and 529: References 505 Hudi, J. AMIGO - A m
- Page 530 and 531: References 507 Palmer, D. Course no
- Page 532 and 533: References 509 European ADAMS News
- Page 534 and 535: Index Abuse loads, 148 3 g bump, 14
- Page 536 and 537: Index 513 Elk Test, 1 El-Nasher, 29
- Page 538 and 539: Index 515 generalised translational
- Page 540 and 541: Index 517 steer axis inclination, 1
Appendix B: Fortran tyre model subroutines 479<br />
C<br />
C<br />
C<br />
C<br />
C<br />
C<br />
C<br />
C<br />
Compute lateral force<br />
CFYA0<br />
SHFYA8*FZPA9A10*GAMMA<br />
DFY(A1*FZPA2)*(1-A15*GAMMA**2)*FZP<br />
IF(ALPHASHFY.LT.0.0)THEN<br />
DUMFY-1.0<br />
ELSE<br />
DUMFY1.0<br />
ENDIF<br />
EFY(A6*FZPA7)*(1-(A16*GAMMAA17)*DUMFY)<br />
BFY((A3*SIN(2*ATAN(FZP/A4)))*(1-A5*ABS(GAMMA)))/<br />
(CFYDFY)<br />
SVFYA11*FZPA12(A13*FZP**2A14*FZP)*GAMMA<br />
PHIFY(1-EFY)*(ALPHASHFY)(EFY/BFY)*<br />
ATAN (BFY*(ALPHASHFY))<br />
FYPDFY*SIN(CFY*ATAN(BFY*PHIFY))SVFY<br />
Change sign<br />
FYFYP<br />
Compute self aligning moment<br />
CTZC0<br />
SHTZC11*FZPC12C13*GAMMA<br />
DTZ(C1*FZP**2C2*FZP)*(1-C18*GAMMA**2)<br />
IF(ALPHASHTZ.LT.0.0)THEN<br />
DUMTZ-1.0<br />
ELSE<br />
DUMTZ1.0<br />
ENDIF<br />
ETZ(C7*FZP**2C8*FZPC9)*(1-(C19*GAMMAC20)*<br />
DUMTZ)<br />
ETZETZ/(1-C10*ABS(GAMMA))<br />
BTZ((C3*FZP**2C4*FZP)*(1-C6*ABS(GAMMA))*<br />
EXP(-C5*FZP))/(CTZDTZ)<br />
SVTZC14*FZPC15(C16*FZP**2C17*FZP)*<br />
GAMMA<br />
PHITZ(1-ETZ)*(ALPHASHTZ)(ETZ/BTZ)*ATAN(BTZ*<br />
(ALPHASHTZ))<br />
TZPDTZ*SIN(CTZ*ATAN(BTZ*PHITZ))SVTZ<br />
C<br />
C Convert to Nmm and change sign<br />
C<br />
TZTZP*1000.0<br />
C<br />
C Copy the calculated values for FX, FY, FZ, TY & TZ to FSAE<br />
C and TSAE arrays<br />
C<br />
1000 FSAE(1)FX<br />
FSAE(2)FY