原理:魔术公式
c*sin(b*arctan(a*x-d*(a*x-arctan(a*x))))
BCD=tan θ=abc
下列操作就是求BCD
c*sin(b*atan(a*x-d*(a*x-atan(a*x))))
用atan
打开carsim,找到轮胎(tire)部分:
求出每个轮胎的垂向载荷(以参考垂直载荷为例圈1),打开2
数据中表格最上面一行代表了不同的垂向载荷,x轴是侧偏角的绝对值,其余部分代表的是侧向力的绝对值。点击Excel,
导出之后的界面是这样的:
第一列是侧偏角,第一行是垂直载荷。
下一步操作是复制选中下列的数据(除了第一行,数据任选)
打开Matlab,
点击 新建变量,将上述复制的数据粘贴过去,保存,然后Matlab工作区就会出现下图的unnamed的数据库
之后在新建脚本,将整列数据分别赋值给x、y,注意一定是列向量。
形如这种:
完成上述操作后,点击运行,将x,y数据导入到工作区
在命令行窗口输入cftool,回车,调用曲线拟合工具。导入x、y选择自定义公式拟合。
+回车
然后选择下列
这是拟合公式,其中BCDE用abcd代替
改完之后,发现曲线并没有拟合各个点,这是因为设置的startpoint不合适,即没有找到合适的系数来拟合曲线。
点击Fit opion,输入合适的参数 如下
然后看下图
曲线拟合有点趋势了arctan不能识别 用atan
Carsim轮胎数据保持基本拟合就好了,一般的d设置为0.95,酌情修改
将abc相乘就是BCD的值,此时的单位是N/deg,最后的单位是N/rad,用得到的值除以pi/180,即得到轮胎侧偏刚度。