使用魔术公式的轮胎模型

使用魔术公式的轮胎模型使用魔术公式的轮胎模型主要有Pacejka’89、Pacejka’94、MF-Tyre、MF-Swift四种。Pacejka’89和’94轮胎模型Pacejka’89和’94轮胎模型是以魔术公式主要提出者H.B.Pacejka教授命名的，根据其发布的年限命名。目前有两种直接被ADAMS引用。魔术公式是用三角函数的组合公式拟合轮胎试验数据，用一套形式相同的公式就可以完整地表达轮胎的纵向力Fx、侧向力Fy、回正力矩Mz、翻转力矩Mx、阻力矩My以及纵向力、侧向力的联合作用工况，故称为“魔术公式”。魔术公式的一般表达式为：BxBxEBxCDxYarctanarctansin式中Y(x)可以是侧向力，也可以是回正力矩或者纵向力，自变量x可以在不同的情况下分别表示轮胎的侧偏角或纵向滑移率，式中的系数B、C、D依次由轮胎的垂直载荷和外倾角来确定。Pacejka’89轮胎模型认为轮胎在垂直、侧向方向上是线性的、阻尼为常量，这在侧向加速度常见范围≤0.4g，侧偏角≤5°的情景下对常规轮胎具有很高的拟合精度。此外，由于魔术公式基于试验数据，除在试验范围的高精度外，甚至在极限值以外一定程度仍可使用，可以对有限工况进行外推且具有较好的置信度。魔术公式正在成为工业标准，即轮胎制造商向整车厂提供魔术公式系数表示的轮胎数据，而不再是表格或图形。基于魔术公式的轮胎模型还有较好的健壮性，如果没有某一轮胎的试验数据，而使用同类轮胎数据替代仍可取得很好的效果。图基于魔术公式的轮胎模型的输入和输出变量Pacejka’89轮胎力与力矩的计算轮胎纵向力计算公式为：VXSBXBXEBXCDF111arctanarctansin其中X1为纵向力组合自变量：X1=(κ+Sh)，κ为纵向滑移率（负值出现在制动态，-100表示车轮抱死）C——曲线形状因子，纵向力计算时取B0值：C=B0D——巅因子，表示曲线的最大值：ZZFBFBD221BCD——纵向力零点处的纵向刚度：ZFBZZeFBFBBCD5423输入变量纵向滑移率κ侧偏角α[rad]外倾角γ[rad]车轮垂直载荷FZ[kN]输出变量纵向力FX[N]侧向力FY[N]翻转力矩MX[Nm]滚动阻力矩MY[Nm]回正力矩MZ[Nm]魔术公式B–刚度因子：B=BCD/(C×D)Sh——曲线的水平方向漂移：109BFBSZhSv——曲线的垂直方向漂移：Sv=0E——曲线曲率因子，表示曲线最大值附近的形状：8726BFBFBEZZ图轮胎属性文件中的纵向力计算系数数据块图Pacejka’89轮胎纵向力示例轮胎侧向力计算公式为：VYSBXBXEBXCDF111arctanarctansin此时的X1为侧向力计算组合自变量：X1=(α+Sh)，α为侧偏角C——曲线形状因子，侧向力计算时取A0值：C=A0D——巅因子，表示曲线的最大值：ZZFAFAD221BCD——侧向力零点处的侧向刚度：5431arctan2sinAAFABCDZB–刚度因子：B=BCD/(C×D)Sh——曲线的水平方向漂移：8109AAFASZh曲线形状因子巅因子计算系数BCD计算系数曲线水平漂移计算系数曲线曲率因子计算系数Sv——曲线的垂直方向漂移：131211AFAFASZZVE——曲线曲率因子，表示曲线最大值附近的形状：76AFAEZ图轮胎属性文件中的侧向力计算系数数据块图Pacejka’89轮胎纵向力示例轮胎回正力矩计算公式为：VZSBXBXEBXCDM111arctanarctansin此时的X1为回正力矩计算组合自变量：X1=(α+Sh)，α为侧偏角C——曲线形状因子，回正力矩计算时取C0值：C=C0D——巅因子，表示曲线的最大值：ZZFCFCD221BCD——回正力矩零点处的扭转刚度：ZFCZZeCFCFCBCD564231B–刚度因子：B=BCD/(C×D)Sh——曲线的水平方向漂移：131211CFCCSZh曲线形状因子巅因子计算系数BCD计算系数曲线水平漂移计算系数曲线曲率因子计算系数曲线垂直漂移计算系数Sv——曲线的垂直方向漂移：171615214CFCFCFCSZZZVE——曲线曲率因子，表示曲线最大值附近的形状：1098271CCFCFCEZZ图轮胎属性文件中的回正力矩计算系数数据块图Pacejka’89轮胎回正力矩示例侧偏刚度（LateralStiffness）侧偏刚度在Pacejka’89和’94轮胎模型中假定是一个常量，在轮胎属性文件的参数PARAMETER数据段中通过LATERAL_STIFFNESS语句设定。侧向形变De：De＝Fy/LATERAL_STIFFNESS；翻转力矩：Mx＝-Fz×De；纵向力和侧偏角联合作用的回正力矩Mz；MZ=MZ,MF+Fx×De，这里MZ,MF为魔术公式计算所得的回正力矩。滚动阻力（Rollingresistance）滚动阻力系数RR同样是在轮胎属性文件中规定的具体值，滚动阻力矩My：My=Fz×Re×RR曲线形状因子巅因子计算系数BCD计算系数曲线水平漂移计算系数曲线曲率因子计算系数曲线垂直漂移计算系数这里：Re为轮胎的滚动半径；RR为滚动阻力系数；Fz垂直载荷（kN）。平滑过渡（Smoothing）是否使用平滑过渡也在轮胎属性文件中规定：USE_MODE=1或2：关闭平滑过渡USE_MODE=3或4：使用平滑过渡轮胎属性文件TR_rear_pac89.tir全文（示例整车模型MDI_Demo_Vehicle.asy使用的）：$---------------------------------------------------------------------MDI_HEADER[MDI_HEADER]FILE_TYPE='tir'FILE_VERSION=2.0FILE_FORMAT='ASCII'(COMMENTS){comment_string}'Tire-XXXXXX''Pressure-XXXXXX''TestDate-XXXXXX''Testtire''NewFileFormatv2.1'$--------------------------------------------------------------------------UNITS[UNITS]LENGTH='mm'FORCE='newton'ANGLE='radians'MASS='kg'TIME='sec'$--------------------------------------------------------------------------MODEL[MODEL]!usemode1234!-------------------------------------------!smoothingXX!combinedXX!PROPERTY_FILE_FORMAT='PAC89'轮胎模型关键词FUNCTION_NAME='TYR900'解算器函数USE_MODE=4.0平滑过渡模式$----------------------------------------------------------------------DIMENSION[DIMENSION]UNLOADED_RADIUS=340.6轮胎自由半径WIDTH=255.0轮胎宽度ASPECT_RATIO=0.35高宽比$----------------------------------------------------------------------PARAMETER[PARAMETER]VERTICAL_STIFFNESS=310.0纵向刚度系数VERTICAL_DAMPING=3.1纵向阻尼系数LATERAL_STIFFNESS=190.0侧偏刚度ROLLING_RESISTANCE=0.0滚动阻力系数$-----------------------------------------------------------LATERAL_COEFFICIENTS[LATERAL_COEFFICIENTS]a0=1.65000a1=-34.0a2=1250.00a3=3036.00a4=12.80a5=0.00501a6=-0.02103a7=0.77394a8=0.0022890a9=0.013442a10=0.003709a11=19.1656a12=1.21356a13=6.26206$-------------------------------------------------------------------longitudinal[LONGITUDINAL_COEFFICIENTS]b0=2.37272b1=-9.46000b2=1490.00b3=130.000b4=276.000b5=0.08860b6=0.00402b7=-0.06150b8=1.20000b9=0.02990b10=-0.17600$----------------------------------------------------------------------aligning[ALIGNING_COEFFICIENTS]c0=2.34000c1=1.4950c2=6.416654c3=-3.57403c4=-0.087737c5=0.098410c6=0.0027699c7=-0.0001151c8=0.1000c9=-1.33329c10=0.025501c11=-0.02357c12=0.03027c13=-0.0647c14=0.0211329c15=0.89469c16=-0.099443c17=-3.336941注意：属性文件中的单位数据块[UNITS]不用于魔术公式的系数a,b,c。