Simple combined slip and brush tire model

Objective: clear all close all clc Fz=4000;mu=0.85;lt=0.2;kt=3.5e6;ky=2.5e6; Cx=lt^2*kt/2; C_alpha=lt^2*ky/2; for i=3:2:11 slip_angle=(i-1)*pi/180; for j=1:201, Skid_ratio=(j-1)*0.005; term1=((Cx*Skid_ratio)^2+(C_alpha*tan(slip_angle))^2 )^0.5; lc=mu*Fz*(1-Skid_ratio)*lt / (2*term1); if lc>=lt Fx(j)=Cx*Skid_ratio/(1-Skid_ratio);…