汽车转弯制动工况轮胎垂直载荷仿真计算
汽车转弯制动工况轮胎垂直载荷仿真计算
农机化研究
第1期
汽车转弯制动工况轮胎垂直载荷仿真计算
张利鹏 1,祁炳楠 2
(1.聊城大学 汽车与交通工程学院,山东 聊城 252059;2.燕山大学 信息工程学院,河北 秦皇岛 066004)
摘 要:在已有的轮胎力学模型的基础上,对汽车转弯制动工况车轮垂直载荷的变动情况作了深入的理论
- 104 -
2006 年 1 月
农机化研究
第1期
力 , F fy 为 前 悬 挂 质 量 横 向 惯 性 力 , Fry 为 后 非 悬 挂 质 量 横 向 惯 性 力 ,Fsx 为 悬 挂 质 量 的 纵 向 惯 性 力 ,Fsy 为 悬 挂 质 量 的 横 向 惯 性 力 ,前 悬 架 侧 倾 角 刚 度 Cεfx , 后 悬 架 侧 倾 角 刚 度 Cεrx , hrool 为 侧 倾 轴 离 地 高 度 , hs
横摆角加速度 侧向速度导致
2
1
1/S
横摆角速度
俯仰角加速度
5
1/S
1/S
侧仰角加速度 俯仰角
悬挂 前轴 变量 1 分量
悬挂质 C
量侧向
加速度
后轴
分量
f(u)
C
前轴 悬挂质量 2 侧彻角
铡度
前轴侧向 变动量
f(u)
1/S
1/S
侧仰角速度 侧仰角
3 侧仰加速度
C 前偏离分量
C 后轴侧倾角刚度 后偏
离分量 f(u)
纵向惯性力引起的左、右后轮载荷变动量为
偏离产生的附加力矩与向侧倾轴移动产生的惯性力 矩 的 合 力 矩 在 后 轴 的 分 量 ; Fryhr 为 后 非 悬 挂 质 量 的
∆z3x = ∆z4x = −∆z1x 。
十一自由度汽车转弯防抱制动系统仿真的研究
制 器,用 MABLAB / SIMULINK 提 供 的 Noniinear ConIroi
Design( NCD)来为 PID 控制器寻优。寻优得到 PID 控制器的
三 个 参 数 值:
Kp,Ki,Kd。
本文仿真
初始条件均设
为制动半径
R0 = 50m,汽 车初 速 V0 = l0m / s, 惯 性
·
误差的变化率 Es 的论域为[ - 0. 6,0. 6],语言变量分为“ 负
大”,“ 负中”,“ 负小”,“ 零”,“ 正小”,“ 正中”,“ 正大”,它们
的隶属度函数取三角形分布,分布情况如图 8 所示。取车轮
分泵制动液油压的改变量 ! 为输出语言变量,它的论域和子
·
项及子项的隶属度函数在论域上的分布与 "# 相同。
仿真结果表明制动过程与实际制动过程相符证明所采用的数学模型的可信以上仿真过程假设工况条件理想而事实上工况条件极其复杂除因系统本身参数变化导致pid参数不准确之外路面凹凸路面突变及车身跳动等均可对系统引入噪音
第 23 卷 第 08 期 文章编号:l006 - 9348(2006)08 - 0222 - 04
将产生不同的制动效果,表现为各车轮的滑移率曲线有较大
差别。两后轮的制动过程较平稳,超调量较小,所有车轮的调
— 224 —
图 6 白噪音工况下 PID 控制的车轮滑移率
由图 6 所示的曲线可知,当系统中引入噪音后,PID 控制 器的控制效果相当差,车轮滑移率远远大于给定值 0. 2。两后 轮在整个制动过程中基本处于危险区域,严重影响车体的方 向稳定性。两前轮从 0. l - 0. 5 秒钟处于危险区域,影响车体 的操作稳定性。针对 PID 控制器的不足,本文设计了鲁棒性 更强的模糊控制器。
各轮胎垂向载荷公式
各轮胎垂向载荷公式
各轮胎垂向载荷公式是指用数学公式计算轮胎在垂向方向上承
受的载荷大小,它与轮胎的结构、材料、气压等因素有关。
一般而言,轮胎垂向载荷公式可以表示为:
Fz = (k1 * P - k2 * Fz0) * (1 - k3 * γ) - k4 * P * Δh 式中,Fz表示轮胎垂向载荷大小,P表示轮胎内部气压大小,Fz0表示轮胎静载荷,γ表示轮胎侧向偏角,Δh表示轮胎变形量,k1、k2、k3、k4为系数。
轮胎垂向载荷公式的正确应用可以帮助车辆提高行驶稳定性、减少轮胎磨损和延长使用寿命。
因此,各车辆制造商都会对轮胎垂向载荷进行专门的测试和计算,以确保车辆在不同路况和负载情况下都能保持良好的性能表现。
- 1 -。
某重型车辆车桥轮毂结构优化后的计算机模拟
某重型车辆车桥轮毂结构优化后的计算机模拟摘要:近年来,对重型车辆的轻量化设计要求越来越高,在满足安全性能要求的前提下,质量稳定的轻量化产品将在竞争中更有影响力。
轮毂是重卡汽车车桥上的十分重要零部件,传统结构设计使用的材料较多,质量较大,结构十分笨重。
本文用有限元软件ansys对某重型车辆轮毂的优化结构进行了有限元分析,对轮毂轻量化设计的研究提供了参考资料。
引言近年来,随着能源和原材料的日趋紧张,对重型车辆的轻量化设计要求越来越高,在满足安全性能要求的前提下,质量稳定的轻量化产品将在竞争中更有影响力。
轮毂是重卡汽车车桥上的十分重要零部件,传统设计中为了保证安全性能,其设计安全系数都很高,所以结构中使用的材料较多,质量较大,结构十分笨重。
本文用有限元软件ansys对某重型车辆轮毂的优化结构进行了有限元分析,对轮毂轻量化设计的研究提供了参考资料。
1.有限元模型的建立某重型车辆轮毂原来质量33Kg,优化轮毂质量为24Kg,三维模型见图1。
a优化前b优化后图1 轮毂的三维模型轮毂、制动鼓、轮辋网格采用单元Solid92,划分网格后的有限元模型如图2所示。
轮辋、轮毂、制动毂的工况见表1。
轮毂内表面与参考点1通过RBE2耦合;制动毂内表面与参考点2通过RBE2耦合;轮毂安装半轴端盖处于参考点3通过RBE2耦合;轮辋与轮胎接触面与轮胎接地点通过RBE2耦合,假设轮胎半径500mm。
计算1分为4个工况,参考点1约束1-5自由度,参考点2约束6自由度,载荷施加在轮胎接地点处,具体如下表所示:表1 轮毂的工况图2划分网格后的有限元模型2.结果分析计算结果见图3、图4、图5。
不同工况下,原轮毂和优化后轮毂计算结果的应力对比见表2。
表2 不同工况原轮毂和优化后轮毂的应力a原轮毂 b优化后的轮毂图3垂直工况a原轮毂 b优化后的轮毂图4侧向工况a原轮毂 b优化后的轮毂图5制动工况通过对比可以看出,优化后的结构在不同工况下的屈服极限大大提高,因此优化后的轮毂性能更好。
三轴车辆行驶中的车轮垂直载荷转移计算研究
(7) 所研究车辆属于不等长的双横臂独立悬 架,考虑到《汽车构造》一书第 227 页的表述“不等 长的双横臂独立悬架,如两臂长度选择适当,可以 使车轮和主销角度以及轮距变化都不大”,忽略悬 挂 系 统 引 起 的 车 轮 轮 距 变 化 [4]。
r
rrr(ms2+mu2)g·dm / 2+Fsy2hs+mu2gayhu=(Fzml′+Fzmr′)dm (8)
r r
rr(ms3+mu3)g·dr / 2+Fsy3hs+mu3gayhu=(Fzrl′+Fzrr′)dr
式 中 :hs— ——簧 载 质 量 的 质 心 高 度 ;hu— ——非 簧 载
根据机械振动理论可以得到由各悬架变形量
对应车轮的垂直载荷为
Fzij =Kij Δij +Cij Δ觶 ij +mui g / 2
(3)
式中:Kij — ——对应车轮的悬架刚度;Cij — ——对应车
轮的悬架减震器阻尼系数;mui — ——各轴所承受的
非簧载质量;g— ——重力加速度。
如图 1 所示,假设车辆行驶方向如图所示,以
本文基于对某三轴轮式装甲车辆转向制动工 况行驶稳定性的研究,建立一个三轴轮式车辆的动 力学模型[3],动力学模型中涉及到垂直载荷的计算, 用于计算出各车轮所受地面的纵向力和侧向力。
1 模型假设
过弯的时候地面要给轮胎侧向力, 保证汽车
收稿日期:2015-04-14 修回日期:2015-04-24
的圆周运动所需的向心力, 相应地车身会产生一 个水平向外的侧向加速度, 这个惯性加速度作用 在悬架上,导致车身侧倾。 悬架系统为抵消车身的 侧倾需要抗侧倾力矩,外侧悬架压缩,内侧悬架拉 伸, 悬架弹簧等弹性元件变形产生的力相应地作 用 在 车 轮 胎 上 ,产 生 了 轮 胎 载 荷 转 移 [4]。
基于ansys workbench的赛车车轮疲劳仿真分析
10.16638/ki.1671-7988.2020.01.035基于ANSYS Workbench的赛车车轮疲劳仿真分析夏衍,陈双贺,田朔源,杨帆,修浩然(辽宁工业大学汽车与交通工程学院,辽宁锦州121001)摘要:赛车车轮是车辆承载的重要安全部件,行驶过程中,赛车车轮承受来自路面不同幅值、不同频率的激励除受垂直力外,还受因车辆起动、制动时扭矩的作用,转弯、冲击等来自多方向的不规则受力。
高速旋转的车轮直接影响车辆的平稳性和操纵性。
文章以Wonder7号铝合金车轮为研究对象,在CA TIA中建立赛车车轮的三维模型,并导入到ANSYS Workbench软件中生成轮辋和轮辐的几何模型。
根据计算极限工况下,对wonder7号车轮进行受力分析,并对车轮的受力载荷进行确定。
建立车轮的有限元模型并进行有限元分析。
为预测车轮的疲劳寿命,用Ansys中的Fatigue模块对车轮进行疲劳寿命分析,预测车轮疲劳破坏位置和使用寿命,对设计人员起了指导意义。
关键词:车轮设计;ANSYS分析;寿命预测中图分类号:U463.34 文献标识码:A 文章编号:1671-7988(2020)01-117-03Fatigue Simulation Analysis of Racing Car Wheels Based on ANSYS WorkbenchXia Yan, Chen Shuanghe, Tian Shuoyuan, Yang Fan, Xiu Haoran( School of automobile and traffic engineering, Liaoning University of Technology, Liaoning Jinzhou 121001 )Abstract: The car wheel is an important safety part of the vehicle load, driving process, the car wheel under the road from different amplitude, different frequency of the excitation in addition to the vertical force, but also by the vehicle starting, braking torque, turning, impact from multiple directions of irregular stress. The high speed rotating wheel directly affects the vehicle's stability and maneuverability. This paper takes Wonder7 aluminum alloy wheel as the research object, establishes the 3d model of racing wheel in CATIA, and imports it into ANSYS Workbench software to generate the geometric model of rim and spoke. According to the calculation of the limit condition, the stress analysis of wonder7 wheel was carried out, and the stress load of the wheel was determined.The finite element model of wheel is established and analyzed. In order to predict the Fatigue life of the wheel, the Fatigue life of the wheel is analyzed with the Fatigue module of Ansys, and the position and service life of the wheel are predicted.Keywords: Wheel design; ANSYS analysis; Life predictionCLC NO.: U463.34 Document Code: A Article ID: 1671-7988(2020)01-117-03引言车轮作为赛车上的重要部件,其可靠性与赛车的行驶安全性紧密相关。
魔术公式——精选推荐
Pacejka 的“魔术公式”是近期在汽车操纵动力学研究中比较流行的公式,它是用特殊正弦函数建立的轮胎的纵向力、横向力和自回正力矩模型。
因只用一套公式就完整地表达了纯工况下轮胎的力特性,故成为“魔术公式”。
“魔术公式”表达如下:sin(arctan((arctan())))vhY y Sy D C Bx E Bx Bx x X S =+=--=+其中,Y 表示侧向力、纵向力或回正力矩,X 表示侧偏角a 或滑移率s 。
现以侧向力为例说明公式中各系数的意义。
式中D --峰值因子,表示曲线的最大值B C D ××--侧偏角趋于零时轮胎的侧偏刚度,曲线中表示原点的斜率E --曲线形态因子,决定曲线最大值附近的形状C --决定曲线的形状,即曲线是像侧向力、纵向力还是回正力矩h S --水平方向漂移v S --垂直方向漂移除C 外,“魔术公式”中的参数都是垂直载荷和侧倾角g 的函数。
以下为轮胎纵向动力学特性在MA TLAB中的具体实现212x z zD a F a F =+1.65x C=2678x zz E a F a F a =++5234()///za F x z z x x B a F a F e C D =+52342122678sin(arctan((arctan())))()///1.65zx x x x x x x a F x z z x xx xz zxzz F D C B E B B B a F a F e C D C D a F a F E a F a F a k k k =--=+==+=++sin(arctan((arctan())))(arctan((arctan())))cos(arctan((arctan())))x xx x x x x z zx x x x xxxxxxxzdF dD C B E B B dF dF d B E B B D C B E B B CdF k k k k k k k k k =--+----2(arctan((arctan())))(arctan())1((arctan()))1((arctan()))x x x x z x x x x x x xx x x x z z z d B E B B dF dB dE d B B B B E B E B B dFdF dF k k k k k k k k k k k --=----+--222(arctan())()11()1()xx xxxx zz x z x z d B B dB dB B dB dF dF B dF B dFk k k k k k k k -=-=++ 555522343422234345(()///)()///(()///)(()///)zzzza F a F xzzx x z z xxz zza F a F x z z x x z z x x zdBd a F a F eC D d a F a F eC D dF dF dF dD a F a F eC D a F a F eC D a dF ++==-+-+sin(arctan((arctan())))x x x x x x x F D C B E B B k k k =--22cos(arctan((arctan())))(())1()1((arctan())x xxxxxxxxxxxxx x x x B D C BE BB C B E B dF B d B E B B kkkk kk k k ----+=+--车辆防抱死制动系统的控制技术研究车辆防抱死制动系统的控制技术研究 东大,侯光钰,张为公东大,侯光钰,张为公2.3.1 H.B.Pacejka 的魔术公式的魔术公式Pacejka 的“魔术公式”是汽车操纵动力学研究中应用比较广泛的轮胎力学模型,它用特殊的正弦函数建立轮胎的纵向力、横向力和回正力矩的函数表达式。
轮胎动态力学响应的计算与仿真
轮胎动态力学响应的计算与仿真随着工业的发展,越来越多的设备、汽车等需要使用轮胎。
而轮胎对于车辆的行驶、安全起着重要的作用。
为了更好地了解轮胎的性能,需要进行轮胎力学的计算与仿真。
本文将介绍轮胎动态力学响应的计算与仿真方法。
1. 轮胎基本构造轮胎是由胶料、纤维材料和金属材料组成的弹性土体结构。
一般分为胎面、肩部和斜肩带、直肩带、侧带及内部骨架几个部分。
其中胎面是与路面接触的部分,受到来自路面的压力和摩擦力,需要承担车辆的荷载和转向、制动等方向控制力。
肩部和斜肩带、直肩带起到转向力的传递作用,侧带主要起到承受侧向荷载和稳定轮胎的作用。
内部骨架则是轮胎的支撑结构,承担着承载力和刚度的作用。
2. 轮胎动态响应轮胎的动态响应包括从外部世界给予的荷载变化中,轮胎内部的变形和应力分布响应。
这些响应影响了轮胎的性能,如附着力、侧偏刚度、阻尼、刚度等。
因此,了解轮胎的动态响应对于轮胎的设计和应用至关重要。
3. 轮胎动态力学计算轮胎动态力学计算主要是解决轮胎与路面接触时的轮胎行驶性能问题。
主要涉及轮胎附着力、摩擦力、刚度和阻尼等参数的计算。
轮胎与路面接触时的附着力主要受到轮胎压力、侧向力和摩擦力的影响。
为了计算这些参数,需要建立轮胎力学模型和路面摩擦模型。
然后通过有限元法或其他数值计算方法求解,得到轮胎的动态响应。
4. 轮胎动态仿真轮胎动态仿真是建立轮胎力学模型,将其放置于车辆整车模型中,进行整车仿真。
仿真结果可以得到轮胎在车辆整体运动过程中的动态响应情况,如附着力、侧偏刚度、刚度、阻尼等参数。
通过仿真结果,可以评估轮胎性能,并进行轮胎设计和优化,提高轮胎的安全性能、保证车辆行驶的稳定性和舒适性。
5. 轮胎动态力学计算与仿真软件目前比较流行的轮胎动态力学计算和仿真软件包括ANSYS、ABAQUS、LS-DYNA等。
其中ANSYS的轮胎模块能够建立轮胎三维模型,进行轮胎刚度、阻尼、附着力、破坏等方面的仿真计算。
而ABAQUS则可以进行轮胎接触模型的建立和刚度、阻尼、摩擦等参数的计算。
用于车辆紧急制动仿真的动态轮胎模型
2007年12月农业机械学报第38卷第12期用于车辆紧急制动仿真的动态轮胎模型邓 李开军 何 乐 夏群生 【摘要】 提出一种可用于车辆紧急制动仿真的动态轮胎模型——L uGre 轮胎模型。
该动态模型不仅具有与经典稳态模型(如魔术公式)相似的稳态特性,可以方便地通过试验数据进行参数拟合和在线整定,而且能够精确捕捉汽车紧急制动过程中的瞬态特性,又由于采用微分方程形式来描述,更有利于开发高性能的电子制动控制系统来提高汽车的行驶安全性。
得到了L uGre 轮胎模型的基本特性,并分析了其在车辆紧急制动仿真中的稳态特性和动态特性。
关键词:动态轮胎模型 紧急制动 仿真中图分类号:U 46315;U 4631341文献标识码:AD ynam ic T ire M odel for Veh icle Em ergency Brak i ng Si m ula tionD eng Kun L i Kaijun H e L e X ia Q un sheng(T sing hua U n iversity )AbstractA dynam ic tire m odel 2L uGre tire m odel w h ich cou ld be u sed in the veh icle em ergency b rak ing si m u lati on w as in troduced .T h is m odel has si m ilar steady 2state characteristics com pared w ith“m agic fo r m u la ”.Its param eters cou ld be easily decided by the experi m en tal tire data and tuned on 2line .L uGre tire m odel cou ld no t on ly cap tu re tran sien t characteristics in the em ergency b rak ing p rocess ,bu t also take advan tage to the developm en t of h igh 2p erfo r m ance electron ic b rak ing con tro l system s to i m p rove the veh icle safety becau se of u sing the fo r m of o rdinary differen tial equati on .T he basic characters of the L uGre tire m odel w ere described and its steady 2state and tran sien t characteristics in the si m u lati on of the em ergency b rak ing p rocess have been analyzed .Key words D ynam ic tire m odel ,Em ergency b rak ing ,Si m u lati on收稿日期:2006-09-11邓 清华大学汽车安全与节能国家重点实验室 硕士生,100084 北京市李开军 清华大学汽车安全与节能国家重点实验室 硕士生何 乐 清华大学汽车安全与节能国家重点实验室 工程师夏群生 清华大学汽车安全与节能国家重点实验室 教授 引言在汽车动力学仿真中,一个合适的轮胎模型可以提供更精确的仿真结果。
轮毂承载能力仿真计算程序设计
轮毂承载能力仿真计算程序设计摘要:本文针对轮毂承载能力的仿真计算问题,设计并实施了一种轮毂承载能力仿真计算程序。
该程序基于有限元分析方法,通过建立轮毂模型、加载模型和边界条件,计算轮毂在实际工况下的最大承载能力。
实验结果表明,该程序能较准确地估计轮毂的承载能力,可为轮毂设计和选择提供参考依据。
1. 引言轮毂作为汽车重要的组件之一,承担着支撑车辆整体荷载的重要任务。
轮毂的承载能力直接影响着车辆的行驶安全和性能。
为了确保轮毂在运行过程中的稳定性和可靠性,需要进行承载能力的仿真计算。
本文旨在设计一种轮毂承载能力仿真计算程序,以辅助轮毂设计和选择。
2. 轮毂模型建立轮毂的承载能力计算首先需要建立轮毂的有限元模型。
通过CAD 软件将轮毂的几何形状进行精确建模,并划分网格单元。
然后,根据轮毂所使用的材料的力学性能参数,将这些参数分配给网格单元。
最后,根据轮毂的实际使用情况和工况,建立加载模型和边界条件。
3. 有限元分析方法本文采用有限元分析方法进行轮毂承载能力的仿真计算。
有限元分析方法是一种数值计算方法,通过将轮毂划分为有限个网格单元,对每个网格单元进行力学分析和计算,最终得到整个轮毂的受力和变形情况。
有限元分析方法具有精度高、灵活性好和计算结果可靠的优点。
4. 轮毂加载模型和边界条件轮毂的加载模型和边界条件是进行承载能力仿真计算的重要输入参数。
通过分析轮毂在实际工况下所受到的载荷,包括垂直载荷、侧向载荷和刹车力等,建立相应的数学模型。
然后,将这些载荷施加在轮毂模型上,并定义边界条件,如约束和支撑条件。
通过这些加载模型和边界条件,可以模拟轮毂在实际工况下的受力情况。
5. 轮毂承载能力的计算在完成轮毂模型的建立和加载模型的定义后,即可进行轮毂承载能力的计算。
通过有限元分析程序,对轮毂模型进行受力分析,得到轮毂在实际工况下的最大承载能力。
承载能力的计算结果可以直观地反映轮毂在受到特定负载时的变形和破坏情况,为设计和选择轮毂提供了重要的参考依据。
汽车转弯制动工况轮胎垂直载荷仿真计算
汽车转弯制动工况轮胎垂直载荷仿真计算
张利鹏;祁炳楠
【期刊名称】《农机化研究》
【年(卷),期】2006(000)001
【摘要】在已有的轮胎力学模型的基础上,对汽车转弯制动工况车轮垂直载荷的变动情况作了深入的理论分析,并进一步建立了仿真模型.这为汽车在该工况以及更复杂工况下的动力学仿真分析奠定了理论基础,方便了以后的计算.
【总页数】3页(P104-106)
【作者】张利鹏;祁炳楠
【作者单位】聊城大学,汽车与交通工程学院,山东,聊城,252059;燕山大学,信息工程学院,河北,秦皇岛,066004
【正文语种】中文
【中图分类】U461.3
【相关文献】
1.适用于汽车转弯制动仿真的轮胎力学模型 [J], 张利鹏
2.汽车转弯制动过程中的制动力动态分配控制 [J], 于彦;谢里阳;何辉;张立军
3.基于模糊PID的汽车高速转弯制动工况控制策略 [J], 王凡;张蕾
4.基于模糊PID的汽车高速转弯制动工况控制策略 [J], 王凡; 张蕾
5.转弯制动工况下某拖挂式房车行驶稳定性分析 [J], 张剑罡;徐希翼;许梁;张允明;苗立东
因版权原因,仅展示原文概要,查看原文内容请购买。
轿车制动性能仿真分析与计算
进行轿车制动性能的仿真分析和计算通常涉及以下步骤和考虑因素:
步骤一:建立轿车制动系统模型
1. 车辆参数收集:收集轿车相关参数,如质量、车轮半径、制动器类型等。
2. 制动系统建模:建立轿车制动系统的数学模型,包括制动盘、刹车片、制动液、制动缸等组成部分。
步骤二:制动力计算
1. 制动力分析:根据制动器和车速等参数计算制动力的大小。
2. 摩擦系数考虑:考虑轮胎与地面的摩擦系数,影响制动力的传递效果。
步骤三:制动距离仿真
1. 刹车距离计算:利用制动力和车辆动力学方程计算制动过程中的制动距离。
2. 不同情况考虑:考虑干燥、潮湿、结冰路面等不同路况对制动距离的影响。
步骤四:热力学分析
1. 制动系统热平衡:考虑制动过程中制动系统的热平衡问题,防止制动器过热失效。
2. 材料特性影响:考虑制动盘和刹车片材料的热特性对制动性能的影响。
步骤五:模拟验证与优化
1. 仿真验证:运用仿真软件(如ADAMS、Simulink等)进行制动性能的仿真验证。
2. 性能优化:根据仿真结果对制动系统进行优化设计,提高制动性能。
以上步骤涉及到多个领域的知识,包括车辆动力学、制动系统设计、热力学等。
在实际工程中,通常需要借助专业的仿真软件和工程计算工具来进行轿车制动性能的分析与计算。
变载荷工况下的轮胎动力学性能试验仿真
变载荷工况下的轮胎动力学性能试验仿真史宗莉;宋万清【摘要】针对变载荷下汽车轮胎动力学特性的复杂性,运用魔术公式轮胎模型对纵向力、侧向力以及回正力矩与纵向滑移率、侧偏角、垂直载荷的关系等轮胎特性进行了分析.用Matlab软件对魔术公式进行动力学仿真,研究了不同垂直载荷下的纵向力、侧向力、回正力矩和相同垂直载荷不同外倾角下的侧向力、回正力矩;分析了联合工况下的制动力和侧向力纵向滑移率的关系.试验结果表明:魔术公式轮胎动力学模型能够很好地模拟轮胎的动力学特性;一定载荷的增加,影响了纵向滑移率与纵向力、侧偏角与侧向力、侧偏角与回正力矩关系的因子数,提高了车辆的稳定性和操纵性.【期刊名称】《上海工程技术大学学报》【年(卷),期】2016(030)001【总页数】5页(P6-10)【关键词】轮胎动力学性能;魔术公式;侧偏特性;纵向滑移率【作者】史宗莉;宋万清【作者单位】上海工程技术大学电子电气工程学院,上海201620;上海工程技术大学电子电气工程学院,上海201620【正文语种】中文【中图分类】U461轮胎作为汽车与路面唯一接触的部件,随着对其力学特性的深入研究,车辆的操作性和稳定性将得到不断地改善和提高[1-2].轮胎力学特性的研究一般依托于轮胎力学模型的建立[3].目前,力学模型主要有理论模型(Gim模型和Fiala模型)、经验模型和半经验模型等3类.然而,精确的理论建模形式复杂,计算困难,简单的理论建模形式精确度较差.作为车辆模型的子模型,无法将其直接进行车辆操作稳定性仿真计算,只能作为定性分析的工具和半经验模型的基础;经验模型是直接从轮胎数据拟合而来,能够很好地吻合试验结果,但纯粹建立在轮胎数据之上,缺乏理论基础,外推性差[4-6];半经验模型的建模是基于理论研究和试验分析两个方面,同时具备了这两种模型的特点,因此,半经验模型不仅充分满足边界条件,而且精度高、外推性好,易于验证[7].UniTire统一轮胎模型和魔术公式轮胎模型是目前轮胎半经验模型的主体[8]. UniTire统一轮胎模型是建立在理论分析和试验研究基础上,以指数形式来拟合轮胎模型,用于预测轮胎的稳态特性[9].魔术公式[10-11]轮胎模型是以三角函数组合的形式来拟合轮胎试验数据,其突出特点为:模型中的系数具有明确的物理意义,数值获取方便,形式简洁,拟合参数较少且能够很好地拟合试验测得的侧向力、纵向力和回正力矩,统一性好、外推性强,适应各种工况下的轮胎受力情况.本文以魔术公式建立的稳态轮胎模型为基础,分析在变载荷工况下,轮胎的纵向滑移特性、侧偏特性和纵滑侧偏混合特性.图1为国际标准化组织(ISO)轮胎坐标系标准.图中,x轴为车轮平面与地平面的交线,车轮前进方向为正方向;z轴为地平面的垂线,垂直向上为正方向;y轴在地平面上,其正方向保证整个坐标系符合右手定则.将图中轮胎运动方向和轮胎滚动方向之间的夹角定义为侧偏角α;车轮的旋转角速度为Ω;轮胎的外倾角为γ;同时,坐标系上给出了轮胎印迹内的六分力,分别为纵向力Fx、侧向力Fy、垂直载荷Fz、翻转力矩Mx、滚动阻力矩My以及回正力矩Mz.魔术公式一般形式为式中:y(x)可以为侧向力Fy,也可以为纵向力Fx,或回正力矩Mz;根据y(x)定义的不同变量,x为轮胎的侧偏角α,或纵向滑移率σ;B为刚度系数;C为形状系数;D为峰值系数;E 为曲率系数;Sh为曲线水平位移;Sv为曲线垂直位移.改变式中系数B、C、D、E时,模拟出不同轮胎在不同路面上的纵向力、侧向力和回正力矩,如图2所示.2.1 纵向力(Fx)的特性稳态条件下轮胎纵向力和纵向滑移率的关系为纵向滑移特性,主要用于车辆在驱动和制动工况下纵向力的预测.由于轮胎的黏弹性,通常造成轮胎在行驶过程中边滚边滑.加速和制动所需要的摩擦力都来自于轮胎滚动速度和行驶速度之间的差值,而该差值正是造成接触区内轮胎滑动的原因.将此差值无量纲化,定义为向前行驶速度的百分比形式,称滑移率[12-13]σ(%),公式如下式中:reff为轮胎的有效半径;ωw为车轮角速度;vx为向前速度.由魔术公式得到单制动情况下纵向力Fx与纵向滑移率σ之间的关系为Fx=D sin[C arctan(BΦ)]Φ=(1-E)σ+arctan(Bσ)E/B式中:C =1.65;B=BCD/(CD);a1~ a8为魔术公式轮胎特定参数.在干燥路面,当外倾角γ=0°,将2、4、6和8 kN等4种变垂直载荷加入轮胎模型,纵向滑移率σ与纵向力Fx的关系,如图3所示.在同一载荷下,初期纵向滑移率和纵向力几乎呈线性关系一同增大.当纵向滑移率达到10%~15%这个区间范围内,纵向力会达到一个最大值,这时纵向驱动性能或制动力性能达到最大.随着纵向滑移率的进一步增加,导致轮胎的不平稳状态,纵向力下降至轮胎纯滑转.不同载荷下纵向力所能达到的最大值随着载荷的增加而增加,因此,适当的增加载荷,车辆的稳定性将提高.2.2 侧向力(Fy)的特性由魔术公式的一般形式,得到在纯转弯工况下侧向力Fy与侧偏角α之间的关系为Fy=D sin(C arctan(BΦ))+SvΦ=(1-E)(a+Sh)+arctan(B(a+Sh))E/B式中:C=1.30;BCD=a3sin(a4arctan(a5Fz))(1-a12|γ|);B=BCD/(CD);Sh=a9γ;γ.当外倾角γ=0°时,垂直载荷Fz为2、4、6 和8 kN时,侧偏角α和侧向力Fy的关系如图4所示.侧偏角α(绝对值)在较小的范围时,轮胎侧向力基本与侧偏角呈线性关系,此时的斜率为侧偏刚度.侧偏刚度的增加,能使得车辆平稳性提高,当侧偏角达到一定值时,轮胎侧向力不再随着侧偏角的增加而变化,轮胎侧向力达到附着极限.在同一侧偏角下,轮胎侧向力随着垂直载荷的增大而增大,侧偏刚度增加,车辆平稳性提高.在垂直载荷Fz为4 kN的情况下,分别在γ=0°、γ=±10°、γ=±15°时,侧向力分析如图5所示.随着车轮外倾角γ的增大,侧偏刚度略微减小,侧向力Fy增大,且具有方向性.2.3 回正力矩(Mz)的特性轮胎发生侧偏时,会产生作用于轮胎绕Ox轴的力矩,该力矩称为回正力矩.轮胎回正力矩的计算公式为Mz=Dsin(C arctan(BΦ))+SvΦ=(1-E)(a+Sh)+arctan(B(a+Sh))E/B式中:C=2.40;(1-a12|γ|);B=BCD/(CD);Sh=a9γ;γ.在不考虑外倾角时,将 2、4、6和8 kN等4种变垂直载荷下的侧向力进行对比,如图6所示.在侧偏角|α|< 4°时,回正力矩与侧偏角几乎呈线性关系.随着侧偏力加大,更多胎面达到附着极限,直到在|α|=4°时回正力矩达到最大值.大侧偏角时回正力矩会发生变号[14-15].当载荷增加一倍,增加的区域都集中在对回正力矩影响最大的区域,结果导致回正力矩随着垂直载荷的增加而加速增加.在载荷Fz=4 kN的情况下,分别在γ=0°、γ=±10°、γ =±15°时轮胎回正力矩如图7所示.回正力矩Mz随着车轮外倾角γ的增大而增大,γ对Mz影响具有方向性.侧偏和纵滑联合工况下,轮胎纵向力Fx与侧向力Fy分别与轮胎侧偏角α、纵向滑移率λ之间的关系为式中:Fx0为单工况下的纵向力;Fy 0为单工况下的侧向力.又有垂直载荷Fz=4 kN时,在α=2°、α=5°、α=8° 不同侧偏角下,侧向力Fy与纵向力Fx关于纵向滑移率σ的关系如图8所示.当车轮同时具有纵向和侧偏滑动工况时,则纵向和侧向两个力分别都要比单工况下的值小,最大侧向力和最大纵向力无法同时取到.在给定侧偏角的工况下,纵向滑移率的增加通常会使侧向力减小;反之,在给定制动工况下,侧偏角通常会使纵向力减小. 采用魔术公式进行轮胎动力学的建模、仿真和分析,结果表明:1) 魔术公式建立的轮胎模型,与实际轮胎受力相比,对侧向力、纵向力和回正力矩都可以取得较高的拟合精度,能够充分地展现汽车轮胎在变载荷下的动力学特性. 2) 车轮的纵向滑移率需控制在较小范围内,使得轮胎纵向力和侧向力都取得较大值,获得车辆侧向稳定性.3) 当侧偏角在较小范围时,侧向力与侧偏角几乎呈线性比例关系,当侧偏角达到一定值时,侧向力不再随着侧偏角变换,处于稳定状态.4) 回正力矩数值较小,对于轮胎动力特性的影响较小,一般进行忽略.5) 随着载荷的增加,轮胎的侧向力、纵向力和回正力矩都会相应增加,影响了纵向滑移率与纵向力、侧偏角与侧向力、侧偏角与回正力矩的关系因子数,一定程度上提高了车辆的稳定性和操纵性.【相关文献】[1] WANG J Y.Theory of Ground Vehicle[M].3rd ed.New York:John Wiley & Sows Inc.,2001.[2] REZA N,JAZAR D.Vehicle Dynamics:Theory and Applications [M].NewYork:Springer,2008.[3] THOMAS D.Fundamentals of Vehicle Dynamics[M].New York:Society of Automotive Engineering Inc.,2003.[4] RAJESH R.Vehicle Dynamics and Control[M].New York:Springer,2005.[5] PACEJKA H B.Tyre and Vehicle Dynamic[M].2nd ed.Burlington:Butterworth-Heinemann Ltd.,2006.[6] 任雷.复杂工况下的轮胎非稳态侧偏特性研究[D].长春:吉林工业大学,2000.[7] BESSELINK I J M,SCHMEITZ A J C,PACEJKA H B.An improved magic formula/swift tyre model that can handle inflation pressure changes[J].Vehicle SystemDynamics,2010,48(s1):337-352.[8] 卢荡,郭孔辉.汽车动力学仿真的轮胎翻倾力矩特性动力学模型[J].科学技术与工程,2004,4(7):601-603.[9] 袁忠诚,卢荡,郭孔辉.UniTire与Magic Formula稳态模型的对比研究[J].汽车技术,2006(2):7-11.[10] KIEBRE R,ANSTETT-COLLIN F,BASSET ing the magic formula model for induced longitudinal force at pure lateral slip[J].International Journal of VehicleDesign,2012(59):305-314.[11] BAKKER G E,PACEJKA H B,LIDNER L.A new tire model with an application in vehicle dynamics studies[C]∥Proceedings of Autotechnologies Conference and Exposition.Monte Carlo:SAE,1989:101-113.DOI:10.42/890087.[12] MACLAURIN B.A skid steering model using the magic formula[J].Journal of Terramechanics,2011,48(4):247-263.[13] KAGEYAMA I,KUWAHARA S.A study on tire modeling for camber thrust and camber torque[J].JSAE Review,2002,23(3):325-331.[14] 郭孔辉,隋军.轮胎垂直载荷分布与侧向力模型及纵向力之间的关系[J].汽车工程学报,1995,32(1):22-25.[15] VELENIS E,TSIOTRAS P,DECANUDASWIT C,et al.Dynamic tire friction models for combined longitudinal and lateral vehicle motion[J].Vehicle SystemDynamics,2005,43(1):3-29.。
桑塔纳2000轿车制动性能与仿真计算
桑塔纳2000轿车制动性能与仿真计算摘要本文首先对轿车制动系统进行了系统的理论分析,汽车制动性能是车辆安全行驶的重要保障,而汽车制动主要由三个方面来评价:制动效能、制动效能恒定性、制动时汽车方向的稳定性。
同过对汽车在平直路面上的制动进行力学分析,建立制动系统数学模型。
对桑塔纳2000轿车的制动力性能做研究,评价整车的制动力安全性,运用Visual C++对汽车的制动性能进行计算机仿真程序设计以及运用该程序对桑塔纳2000轿车的制动力分配,制动距离进行仿真计算。
程序着重于计算机仿真建模和图形结果的输出。
输出结果能够对汽车的制动性能进行有效评价,并对其进行优化设计,对汽车制动性能的计算机模拟仿真具有一定的参考作用,从而优化车辆的制动性能。
关键词:桑塔纳2000 制动性 Visual C++ 仿真Santana 2000 vehicle braking performance and SimulationABSTRACTThe first on the brake system of the system of theoretical analysis, The brake performance is the important guarantee of the traffic safety, and the brake of the automobile consists mainly in three ways to evaluate:brake efficiency, homeostasis of the brake efficiency , the stability of the braking car,s direction . By the analysis of the brake sexual mechanics we build mathematical model。
We research santana 2000,s performance of the power so that evaluate the system power security . Make use of visual C to design on the computer simulation program any distribute santana 2000 automobile, the performance of the brake The program concentrate on modeling and the results output。
后桥桥壳强度计算(垂直、牵引、制动、侧滑工况)
后桥桥壳强度计算(垂直、牵引、制动、侧滑工况)后桥桥壳强度计算(垂直、牵引、制动、侧滑工况)根据《汽车车桥设计》的方法进行计算简算。
数据仅供参考。
这种计算只适合设计初期的粗略计算判断。
实际设计时,需结合有限元分析软件,以处理桥壳的细节尺寸,使应力分布在更合理的状态。
本例子计算的桥壳结构如下截图所示(悬架按普通板簧悬架,车辆承载的力作用于桥壳方截面上面的的板簧盖板上-下图中未画出上盖板):说明:折算到桥壳板簧座处的受力时,应该要用10T减去桥总成的重量的。
不过本例子中是按作用在板簧座处的力为额定载荷,到轮胎上时,是额定载荷+桥重量。
这样算更保险点。
相当于叠加了一点安全系数。
注2:上述表格中的计算都是基于水平路面进行的计算。
所以计算结果仅能作为基础参考数据,起到数据统计对比的价值。
如果车桥使用的路况很恶劣,需要额外考虑坡度、凹坑、凸起等其他因素。
附件-计算表格:桥壳强度计算.xlsx项目代号单位数值两板簧座之间的弯矩M M Nm 28665桥壳截面宽B mm 135 桥壳板簧截面处的静弯曲应力σwj Mpa 98.7因是垂直静弯曲强度计算,所以按2.5倍计算。
地面对后驱动桥左轮的垂向反作用力Z2R N 88519.9地面对左右驱动轮的最大切向反作用力Pmax N 142135.4重力加速度g m/s^2 9.8驱动桥承受的侧向力P2 N板簧对桥壳的垂向作用力-左侧T2L N 1565.8原则上讲a+b 的值越大越好。
但是受空间和质量的限制,又不能做的太大,所以一般情况按a+b ≈rr/4。
离)。
因为此值一般都比较小,所以就省略了。
A-A 截面的垂向弯矩 M Nm 41743.0总结:上述计算的汇总信息如下:例如:公路用桥时,1为*.*g,2为*.*g,3为*.*g,4为*.*g。
(*.*为某一统计经验数值(即几点几个g。
g代表额定载荷)。
当有多个方向受力时,需要叠加各个实际载荷方向的数值)恶劣路面用桥时可在上述条件下叠加一些载荷,或是要求更高的安全系数,即在相同载荷条件下,允许的最大应力值必须更小一些。
弯道制动工况下车辆参数估计
弯道制动工况下车辆参数估计谭子胡;龚贤武;刘峰云;高齐;赵轩【摘要】[目的]为了准确估计车辆在弯道制动工况下的状态参数,观察松弛长度对车轮侧向力的影响,提出了一种基于无迹卡尔曼滤波器(UKF)的状态参数估计模型.[方法]该模型考虑了单个车轮驱/制动力对车辆稳定性的影响,将松弛长度作估计量来调节侧向力的估算精度.选择典型工况,应用Carsim与Matlab/Simulink联合仿真对模型进行验证.[结果]仿真表明,该模型估计的质心纵向速度、横摆率、车轮侧向力与参考数据100%拟合,质心侧向速度误差较大.[结论]该模型能有效地估计弯道制动工况下的车辆状态参数.【期刊名称】《甘肃农业大学学报》【年(卷),期】2019(054)003【总页数】8页(P178-185)【关键词】参数估计;车辆;弯道制动工况;无迹卡尔曼滤波器;松弛长度【作者】谭子胡;龚贤武;刘峰云;高齐;赵轩【作者单位】长安大学汽车学院,陕西西安710064;长城汽车股份有限公司,河北保定071000;长安大学电控学院,陕西西安710064;长安大学汽车学院,陕西西安710064;长安大学汽车学院,陕西西安710064;长安大学汽车学院,陕西西安710064【正文语种】中文【中图分类】V467.1+1随着电动汽车的发展,驱动电机的冗余控制使汽车在驱动力矩的分配方案上具有可控性.根据单个驱动电机的力矩反馈信号能估算出驱动轮的驱/制动力[1],相对于传统车参数估计模型而言,用准确的输入量代替未知参数,进一步提高了参数估计模型的估计精度.随着高精度GPS的应用,信息融合的方法得到广泛运用,通过载体姿态测量的方法直接获得车辆的横摆角或质心侧偏角[2].纵向车速、横摆角速度、质心侧偏角、纵向加速度以及侧向加速度可通过不同的传感器配置方案来测量与修正[3].这些能准确测量的参数可以作为参数估计模型的参考量,从而简化参数估计模型和提高估计精度.轮胎作为车身与路面的载体,准确估计车轮作用力,能提前预测车辆的运动状态和控制车辆.其中,轮胎侧向力是影响车辆横向稳定性的关键因素,对车辆的操纵稳定性具有重要意义[4].如何提出简单又相对准确的模型来估计车轮作用力和其它参数,成为研究的焦点.国内外学者提出了不同的车辆参数估计模型对侧向力进行估计,但未考虑制动下从动轮的制动力作用.Katriniok等[5]将经验轮胎公式引进模型的状态方程中,使状态方程离算化变得困难,计算量大.Doumiati和Baffet等[6-7]忽略从动轮纵向力,将驱动力之和作为估计量,利用扩展开卡尔曼滤波算法估计侧向力,曲线拟合程度较高,但在峰值处误差较大.Doumiati等[8]利用扩展卡尔曼滤波算法对侧倾时垂直载荷进行估计,然后再对侧向力进行估计,但未考虑车辆制动状态下从动轮的制动力作用.在国内,大多文献以四轮驱动车辆为研究对象,未考虑后轴驱动车辆从动轮的制动力对车辆稳定性的影响[9-11].解少博等[12]采用UKF滤波器对23个参数进行估计,将前轮纵向力也作为估计量,但计算量大,噪声参数不易调节.在Doumiati和Baffet[6-8],都采用松弛轮胎模型,将松弛长度作为常数,然而松弛长度不易测量.在Vantsevich等[13]中指出,松弛长度随垂直载荷与速度的变化而变化,Maurice[14]指出松弛长度随着侧偏角的大小而变化.因此,如何取松弛长度初值以及将松弛长度作为常数对车辆参数模型有什么影响,需要进一步观测与验证.上述文献中,忽略了从动轮的纵向力对车辆稳定性的影响;进行车轮驱动力比例分配时存在较大误差,同时没有考虑弯道制动工况下参数的估计效果;将松弛长度作为常数,未验证参考值的准确性.然而,在制动过程中从动轮的制动力对车辆影响比较大,在制动工况下忽略从动轮制动力会带来较大误差.针对以上问题,本文提出了基于UKF的前轮转向后轮驱动的电动车车辆参数估计模型,将驱动轮的驱/制动力作为已知量,减少驱动轮纵向力分配误差;考虑从动轮的制动力对车辆稳定性的影响;将松弛长度作为状态量,分析松弛长度对侧偏力的影响,在此基础上调节松弛长度初值.最后,通过弯道制动双移线工况和高速匀速双移线工况来验证观测器的合理性.1 动力学模型1.1 车辆模型七自由度的整车模型是研究四轮汽车动力学常用的模型[15-16],如图1所示.车辆模型为前轮转向、后轮驱动的双轨道车辆模型,忽略车轮的滚动阻力、空气阻力的影响,包含纵向运动、侧向运动、横摆运动和车轮的滚动运动.图1 后轮驱动的七自由度整车模型Figure1 The seven freedom degrees of rear drive vehicle model根据模型建立x、y、z 3个方向上的方程组,车辆质心纵向加速度、侧向加速度、横摆力矩方程、运动学方程、车轮侧偏角及车轮垂直载荷方程分别为:ax=(Fxflcosδfl+Fxfrcosδfr-Fyflsinδfl-Fyfrsinδfr+Fxrl+Fxrr)/m(1)ay=(Fxflsinδfl+Fxfrsinδfr+Fyflcosδfl+Fyfrcosδfr+Fyrl+Fyrr)/m(2)Mz=ls(Fxflcosδfl+Fxfrcosδfr-Fyflsinδfl-Fyfrsinδfr+Fxrr-Fxrl)+as(Fyflccosδfl+Fyfrcosδfr+Fxflsinδfr+Fxfrsinfr)-bs(Fyrt+Fyrr)(3)(4)(5)(6)式中:ax为纵向加速度,ay为侧向加速度,Mz为横摆力矩,Fxfl、Fxfr、Fxrl、Fxrr分别为左前轮、右前轮、左后轮、右后轮纵向力,Fyfl、Fyfr、Fyrl、Fyrr分别为左前轮、右前轮、左后轮、右后轮侧向力,δfl、δfr为左前轮、右前轮转角.(7)(8)式中:αfl、αfr、αrl、αrr分别为左前轮、右前轮、左后轮、右后轮侧偏角,vx、vy、γ分别为纵向速度、侧向速度和绕z轴的横摆角速度.FZfl、FZfr、FZrl、FZrr 为车轮的垂直载荷,ax和ay分别为纵向加速度和侧向加速度,hg为质心高度.1.2 从动轮制动力估计车辆为前置后驱电动车,从动轮纵向力可以根据dugoff轮胎模型计算出,驱动轮纵向力可以通过电机反馈信号及液压制动力矩信号来获取,数值更加准确.Dugoff 轮胎纵向力表达式如下:(9)式中,λ表示纵向滑移率,Cy表示轮胎侧偏刚度,Cx表示轮胎纵向刚度,α表示轮胎侧偏角,μ为路面附着系数,ε为速度影响系数,是一个与轮胎结构和材料有关的参数,它可以修正车轮滑移速度对轮胎力的影响[17].2 状态观测器设计2.1 松弛轮胎模型为了提高侧向力的估计精度,观测器单独考虑每个轮胎的侧向力,引入松弛轮胎模型[18],其定义如下:(10)式中,Cyij是侧偏刚度系数,αij为车轮侧偏角,σij为松弛长度,μ为摩擦系数,Fzij为轮胎垂直载荷,为车轮侧向力近似值.上述轮胎模型忽略纵向滑移率下的纯侧滑工况的简化模型,统一的轮胎胎压分布以及固定的轮胎/路面附着系数.在dugoff原始轮胎模型中,侧偏刚度为不受垂直载荷转移影响的常数.根据文献[8]表明,轮胎垂直载荷转移对侧偏刚度系数有较大的影响,变化的垂直载荷可以表达成关于垂直载荷的多项式.将侧偏刚度拟合成关于垂直载荷的多项式,最终的线性拟合公式为:Cyij(x)=0.001 214x2-15.34x-2 550(11)2.2 基于UKF的车辆参数模型设计由纵向速度、侧向速度、横摆率、侧向力和松弛长度构成状态向量如下:x(k)==由转向角、纵向力、侧向力近似值构成控制向量如下:=由纵向加速度、侧向加速度、横摆率和后轴平均速度(后轴各轮的速度可以通过速度传感器测出)构成观测向量如下:y(k)==[y1,y2,y3,y4]T为了估计前后轴松弛长度,引入状态模型如下:(12)将状态向量和控制向量带入到公式(4)、(5)、(6)、(10)、(12)中,构成状态方程如下:f1=[(u2+u3)cosu1-(x4+x5)sinu1+u4+u5]/m+x2x3f2=[(u2+u3)sinu1-(x4+x5)cosu1+x6+x7]/m-x1x3f3={l5[(-u2+u3)cosu1+(x4-x5)×sinu1+x4-x5]+as((x4+x5)×cosu1+(u2+u3)sinu1)-bs(x6+x7)}/Jzf6=0f7=0(13)将后轴平均速度作为观测量,后轴速度可以通过车轮上速度传感器获得,并与公式(1)、(2)、(6)构成观测方程组:(14)3 UKF原理无迹卡尔曼滤波器采用不敏变换,可以将估计参数的均值和协方差向真值逼近到三阶,精度与二阶的泰勒级数一致,比扩展卡尔曼滤波器(extended kalman filter,EKF)线性化误差更小一些[19].UKF观测器设计如下:(1) 采样点权值:(15)(2) 状态更新状态:(16)(3) 向前一步预测方差:(17)(4) 观测预测方程:(18)(5) 滤波器增益:(19)(6)状态估计方差:P(k+1|k+1)=P(k|k+1)-K(k)Pz×KT(k)(20)在2次试验中,将α=0.001,β=2,n=9,路面附着系数μ=0.85,纵向初始速度值与试验一致,其它变量设置较小的任意值.2次试验变量的初值为x0=[21 0 0 0 0 0 0 0.4 0.4]T,模型过程噪声和测量噪声对估计结果有重要影响,取Qm=diag[10^(-7),10^(-7),10^(-7),10^-6,10^-6,10^-6,10^-6,10^(-7),10^(-7)]′;Qs=diag[0.9,0.9,10^(-3),1]′,状态初始协方差为P0=I9*9的单位矩阵.4 试验验证与分析Carsim软件是专门针对车辆动力学的仿真软件,被广泛地应用于现代汽车控制系统的开发[20].本文采用在Simulink环境中搭建参数估计模型,在Carsim仿真软件中模拟驾驶环境,2个软件之间进行参数输入与输出,如图2所示.由上述可知,在该试验中纵向驱/制动力为已知量.从图3可以看出,弯道制动双移线试验中从动轮制动力(前轴)在制动过程中比较大,甚至出现峰值,因此从动轮的制动力对车辆的稳定性有很大的影响.图2 基于Matlab/Simulink与Carsim联合仿真平台Figure 2 The simlink platform based on the co-simulation Matlab/Simulink and Carsim图3 车轮驱/制动力Figure 3 The driving /braking force of tire4.1 弯道制动双移线试验为了验证观测器在弯道制动工况下的稳定性与准确性,选取典型双移线工况,初始速度为80 km/h,地面摩擦系数为0.85,以匀减速进行试验.图4为车辆质心位置的3个参数估计曲线,从图中可以看出,质心纵向速度和横摆率的估计值与参考值接近100%重合,侧向速度存在偏差.侧向速度的估计值产生偏差的原因:纵向速度远大于侧向速度,接近100倍,而侧向速度估计方程是关于纵向速度的函数,因此,过小的侧向速度,在计算过程中产生累计误差;另一方面,观测方程中纵向加速度为近似值,忽略横摆率的影响,进一步增加了误差;待估计的纵向力与侧向速度的绝对值大小而言,侧向速度值更小,在综合UKF增益的作用下,进一步增加了侧向速度误差.图4 弯道制动工况下车辆质心参数估计结果Figure 4 Estimated parameters of gravity center under streering and braking condition从图5可知,侧向力的估计值与参考值完全重合,说明设计的观测器能准确的估计弯道制动工况下的侧向力,松弛长度能更好地体现轮胎的侧向动力学特性.从图6可知,EKF观测器中松弛长度变化不大,可以当作常数,初值大小对观测器的影响不大;而在UKF观测器中,先出现峰值再向初值靠近,是一个动态变化的过程.出现峰值的原因是在峰值点处,车辆正处于转弯过程中,侧偏力突然增大,存在动态调节的过程.观测器动态调节松弛长度的过程,更能体现轮胎侧向特性变化.4.2 高速弯道双移线试验车辆在高速运动时对参数估计模型的稳定性较高,为了更好的验证观测器的稳定性与准确性,选取典型高速双移线工况,初始速度为80 km/h,地面摩擦系数为0.85,以恒定速度进行试验,行驶过程不进行减速.从图7~9可知,纵向速度、横摆率、车轮侧向力的估计值与参考值接近100%重合;当车辆在转弯时,侧向加速度接近0.4 g时,侧向力估计值仍能与参考值吻合;侧向速度与松弛长度的变化规律与弯道制动工况下的变化规律一样,试验结果表明该观测器在高速时也具有较好的观测效果.图7可知,当车辆高速行驶时,侧向速度变化比弯道制动工况下估计值波动大,这是因为车辆在弯道高速行驶时,车辆横摆力矩变化大,导致侧向速度波动大.图5 弯道制动工况下侧向力估计结果Figure 5 Estimation of lateral forces under streering and braking condition图6 弯道制动工况下松弛长度估计结果Figure 6 Estimation of relaxation length under streering and braking condition图7 高速弯道工况下车辆质心参数估计结果Figure 7 Estimated parameters of gravity center under streering and high speed condition图8 高速弯道工况下侧向力估计结果Figure 8 Estimation of lateral forces under streering and high speed condition图9 高速弯道工况下松弛长度估计结果Figure 9 Estimation of relaxation length under streering and high speed condition5 结论在设计前轮转向后轮驱动的电动汽车参数估计模型时,将驱动轮的驱/制动力作为已知状态量,能有效避免驱动力分配误差;考虑从动轮的纵向力,更能体现制动过程车辆的受力状况.在EKF观测器中,松弛长度变化不大,可以当作常数;在UKF观测器中,松弛长度在初始调节过程中出现波峰,但很快接近初值.将松弛长度作为状态量,能根据估计效果调整初值,同时也能体现车轮的侧向动力学特性,提高侧向力的估计精度. 通过弯道制动双移线实验和弯道高速恒速双移线试验,试验结果表明,设计的观测器能有效地观测车辆在制动和高速工况下的状态参数.参考文献【相关文献】[1] 褚文博,罗禹贡,陈龙,等.分布式电驱动车辆的无味粒子滤波状态参数联合观测[J].机械工程学报,2013,49(24):117-127.[2] 谢伯元,王建强,秦晓辉,等.基于车路协同的车辆状态估计方法[J].汽车工程,2014,36(8):968-973.[3] 余卓平,高晓杰.车辆行驶过程中的状态估计问题综述[J].机械工程学报,2009,45(5):20-33.[4] Davoodabadi I,Ramezani A A,Mahmoodi-k M,et al.Identification of tire forces using dual unscented kalman filter algorithm[J].Nonlinear Dynamics,2014,78(3):1907-1919. [5] Katriniok A,Abel D.Adaptive EKF-based vehicle state estimation with online assessment of local observability[J].IEEE Transactions on Control Systems Technology,2016,24(4):1368-1381.[6] Doumiati M,Ghandour R,Charara A,et al.Design,experimental validation,and comparison of two model-based EKF observers for lateral vehicle dynamics estimation;proceedings of the 2013[C]∥ 21st mediterranean conference on control and automation.IEEE Computer Society,2013.[7] Baffet G,Charara A,Lechner D,et al.Experimental evaluation of observers for tire-road forces,sideslip angle and wheel cornering stiffness[J].Vehicle System Dynamics,2008,46(6):501-520.[8] Doumiati M,Victorino A,Lechner D,et al.Observers for vehicle tyre/road forces estimation:Experimental validation[C]∥Selected Papers from the 9th International Symposium on Advanced Vehicle Control.Taylor and Francis Ltd.,2010.[9] 宋传学,肖峰,刘思含,等.基于无迹卡尔曼滤波的轮毂电机驱动车辆状态观测[J].吉林大学学报(工学版),2016,46(2):333-339.[10] 付翔,孙威,黄斌,等.四轮独立驱动车辆状态参数的UKF估计[J].武汉理工大学学报(交通科学与工程版),2017,41(3):379-384.[11] 张志勇,张淑芝,黄彩霞,等.基于UKF的分布式驱动电动汽车轮胎侧向力估计[J].西华大学学报(自然科学版),1-6.[12] 解少博,林程.基于无迹卡尔曼滤波的车辆状态与参数估计[J].农业机械学报,2011,42(12):6-12.[13] Vantsevich V V,Gray J P.Relaxation length review and time constant analysis for agile tire dynamics control[C]∥International design engineering technical conferences and computers and information in engineering conference.American Society of Mechanical Engineers (ASME),2015.[14] Maurice J P,Pacejka H B.Relaxation length behaviour of tyres[J].Vehicle System Dynamics,1997,27(sup001):339-342.[15] Li L,Jia G,Ran X,et al.A variable structure extended Kalman filter for vehicle sideslip angle estimation on a low friction road[J].Vehicle System Dynamics,2014,52(2):280-308.[16] 张一西,马建,唐自强,等.分布式驱动电动汽车直接横摆力矩控制策略与仿真验证[J].甘肃农业大学学报,2018,53(4):159-167,176.[17] Dugoff H,Fancher P S,Segel L.[C]∥An analysis of tire traction properties and their influence on vehicle dynamic performance.SAE International,1920.[18] Rajamani R.Vehicle dynamics and control[M].NewYork:Springer,2005.[19] 金学波.Kalman 滤波器理论与应用:基于MATALB实现[M].北京:科学出版社,2016.[20] 熊璐,陈晨,冯源.基于Carsim/Simulink联合仿真的分布式驱动电动汽车建模[J].系统仿真学报,2014,26(5):1143-1148,1155.。
极限工况下汽车操稳性下的轮胎模型
FY 2 bFY
2
考虑到δ角较小 cos 1
FY1 k11 FY 2 k22
FY M
k11 k22 Z ak11 bk2
2
线性二自由度汽车模型对前ห้องสมุดไป่ตู้角输入的响应
质心侧偏角 v
v ar u
u
ar
• 在高车速下,突然大角度急打方向盘,前、后轮侧 偏角突然迅速增大,达到轮胎的饱和状态,轮胎进入 非线性状态,导致各状态变量发生急剧的变化,横摆 角速度急剧增大,质心侧偏角急剧增大,与方向盘转 角成非线性关系。汽车失去稳定性,发生明显的侧 滑,不能正常行驶,容易发生严重的交通事故。
• 轮胎的侧偏力是汽车赖以实现方向变换的 原动力,也是保证行驶方向稳定性的重要条 件;由于轮胎的侧偏特性直接影响着汽车曲 线行驶时的稳态响应特性,通过改善轮胎的 侧偏特性可以达到改善汽车操纵稳定性的 目的。
• 极限工况:
大转向 高速行驶(80码以上) 湿滑路面。
• 汽车的操纵稳定性
包含相互联系的两个部分,一是操纵性,二是稳定性。 操纵性是指汽车能够确切地响应驾驶员转向指令的能力; 稳定性是指汽车受到外界路面挠动或突然阵风挠动后恢复原来 运动状态的能力;
• 轮胎模型:
由于轮胎在行驶过程中的受力非常复杂,根据路面性质、车速、 垂直载荷、摩擦产生的温度以及轮胎的形式等因素的变化而建 立了轮胎模型。
轮胎力建模
• 前轮转角
• 车速
轮
• 轮速
• 横摆角速度
胎
• 外倾角
模
• 轮胎垂向载荷
• 路面附着系数
型
侧向力 纵向力 回正力矩
轮胎模型的分类及各个典型模型
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
∆z3y = − (Fsyrhroll + Cεrxε + Fryhr ) Br
(5)
其 中 , Fsyrhroll 为 悬 挂 质 量 惯 性 力 矩 在 后 轴 的 分
量 , 有 Fsyf hroll = −M sasyahroll L ; Cεrxε 为 悬 挂 质 量 质 心
∆z2x = ∆z1x
本文深入地分析了汽车在转弯制动工况下,各 个车轮垂直载荷的变化情况,并对可能发生翻车事 故的状态进行了探讨,使车轮载荷的变化更加符合 现实情况。在此基础上,建立了车轮载荷变动的仿 真力学模型,此模型比较全面地反映了汽车在转弯 制动甚至更复杂的运动状态下,车轮载荷的变化,
图 5 整车载荷变动分析
参考文献:
分 析 ,并 进 一 步 建 立 了 仿 真 模 型 。这 为 汽 车 在 该 工 况 以 及 更 复 杂 工 况 下 的 动 力 学 仿 真 分 析 奠 定 了 理 论 基 础 ,
方便了以后的计算。
关键词:交通运输工程;垂直载荷;动力学仿真;转弯制动
中 图 分 类 号 :U461.3
文 献 标 识 码 :A
侧向速度 俯仰角速度俯仰角 5
f(u) 后非悬挂力矩
图 3 纵向载荷变动量仿真模型
(u) 纵向变动
1 纵向变动量
3.2 由侧向惯性力引起的轮胎载荷变动量
侧向惯性力引起的左前轮载荷变动量为
∆z1y = − (Fsyf hroll + Cεfxε + Ffyhf ) B f
(4)
式中
Fsyf hroll — 悬 挂 质 量 惯 性 力 矩 在 前 轴 的 分 量 ; Fsyf hroll = −Msasybhroll L ;
小 于 右 侧 轮 胎 ,所 以 单 考 虑 左 侧 轮 胎 垂 直 载 荷 情 况 。
当左侧两个轮胎载荷都为零后,车辆继续行驶下去
必然导致翻车,造成破坏性后果。因此当左侧两轮
胎载荷之和等于零即停止运行仿真。最终在限定的
初始条件下,可以根据受力情况建立各车轮的垂直
载荷仿真模型,如图 5 所示。
4 结论
∆z2 y = −∆z1y
Ff xhf 、 Frxhr — 非 悬 挂 质 量 的 惯 性 力 矩 , 有 Ff xhf = −mf hf (v&x − rvy − r2a) ; Frxhr = −mrhr (v&x − rvy + r2b) 。 纵向惯性力引起的右前轮载荷变动量为
侧向惯性力引起的左后轮载荷变动量为
[1] 张利鹏.汽车转弯制动性能分析与防抱控制仿真研 究[D].秦皇岛:燕山大学,2004.9-20.
[2] 程 军.用 MATLAB/SIMULINK 进行车辆控制系统的 设计[J].汽车研究与开发,1997,(3):7-10.
Tires Vertical Load Computation of Vehicle in Cornering with Braking
(1)
Fz3 = Fz4 = Wc ⋅ a 2L
(2)
a
b
as
s c
Fz (1,2)
hf
Wc Fz(3,4)
hr
L
图 1 静态垂直载荷分析
3 由惯性力引起的轮胎载荷变动量
汽车在转弯制动过程中,由于受到惯性力和质 心偏移的影响,使得各轮垂直载荷发生了很大的变 动,可分为由纵向惯性力引起的轮胎载荷变动量和 由侧向惯性力引起的轮胎载荷变动量两种,如图 2 所示。当悬挂质量处于静平衡位置时,以通过悬挂 质 量 质 心 S 垂 直 向 下 的 铅 锤 轴 为 z1 轴 ,向 上 为 正 ;z1 轴侧倾轴交点为 K ,以 K 为坐标中心,以通过 K 的 对 称 面 内 的 水 平 线 , 即 侧 倾 轴 作 为 x1 轴 , 由 右 手 定 则 确 定 y1 轴 , 此 坐 标 系 与 非 悬 挂 质 量 动 坐 标 系 xoyz 平 行 。 另 取 一 套 动 坐 标 系 x′Ky′z′ , 此 坐 标 系 固 结 在 悬挂质量上,随悬挂质量运动而运动。悬挂质量系 统相对于非悬挂质量系统的位置由绕侧倾轴线的转 角 ε 和绕俯仰轴线的转角 λ 确定,侧倾角速度为 p , 侧 倾 角 加 速 度 为 p& , 俯 仰 角 速 度 为 q , 俯 仰 角 加 速 度 为 q& , 各 转 角 、 角 速 度 和 角 加 速 度 均 以 右 手 定 则 确 定 正 方 向 。车 身 分 别 绕 x1 轴 进 行 侧 倾 运 动 ,绕 y1 轴 进 行 俯 仰 运 动 ,绕 z 轴 进 行 横 摆 运 动 。 Ffx 为 前 非 悬 挂 质 量 纵 向 惯 性 力 , Frx 为 后 非 悬 挂 质 量 纵 向 惯 性
2006 年 1 月
农机化研究
第1期
汽车转弯制动工况轮胎垂直载荷仿真计算
张利鹏 1,祁炳楠 2
(1.聊城大学 汽车与交通工程学院,山东 聊城 252059;2.燕山大学 信息工程学院,河北 秦皇岛 066004)
摘 要:在已有的轮胎力学模型的基础上,对汽车转弯制动工况车轮垂直载荷的变动情况作了深入的理论
λ
∆Z(1,2) x
Z′ Z1
Z
ay
∆Z(3,4) x
Fsyε y( y1 )
y′
ε
Fsy
M sy
S K (C)
Z′ ∆Z(2,4) y Z (Z1)
F1 y ∆Z(1,3) y
图 2 载荷变动分析
3.1 由纵向惯性力引起的轮胎载荷变动量
asy — 侧 向 加 速 度 ;
纵向惯性力引起的左前轮载荷变动量为
(1) 垂直载荷作用在车轮中心平面内,轮距和 轴距均为定值;
(2) 悬 挂 质 量 质 心 相 对 于 侧 倾 轴 线 的 高 度 作 为静态量来考虑;
(3) 非悬挂质量相对于地面无侧倾,仅考虑悬 挂质量的侧倾运动。
汽车动力学仿真是随着计算机软、硬件技术的 发展而发展起来的。作为一种图形化界面的仿真建 模软件,基于 MATLAB 语言环境的 SIMULINK 软件, 用户界面友好,各种功能模块化,可以直接用鼠标 拖放模块,建立信号连接,进行建模,是目前工程 界常用的仿真工具[2]。本文将用其建立汽车转弯制 动工况车轮垂直载荷的仿真模型。
∆z1x = (Fsxhroll + M sx + Ffxhf + Frxhr ) 2L
(3)
式中
Fsxhroll — 悬 挂 质 量 的 惯 性 力 矩 , Fsx = −M sasx ; asx — 纵 向 加 速 度 ;
Ffyhf — 前 非 悬 挂 质 量 的 惯 性 力 矩 , 有 Ffyhf = −m f hf (ay + ar&) ;
为进行比较复杂工况的汽车系统动力学仿真研究提 供理论基础,以避免重复的计算。
整车质量
C
纵向速 度导数
1 俯仰角 加速度
5 侧向 速度
7 侧向速 度导数
2
3 横摆角 加速度
4 侧仰角 加速度
5
纵向 速度
f(u)
In Dut1 if{ }
左前
1
u1>0) u1
else
f(u)
轮载荷 1
3
假设
5 纵变 4 7 纵向变动量
为 静 止 状 态 悬 挂 质 量 质 心 离 地 高 度 ,∆zix(i=1~4) 为 由 纵 向 惯 性 力 引 起 的 4 个 车 轮 载 荷 变 动 量 ,∆ziy(i=1~4) 为 由 横向惯性力引起的 4 个车轮载荷变动量。
ax
λ
Ffx
X(X1) X
L
a
b
M sx Fsx
S K C Frx
文 章 编 号 :1003─ 188X(2006)01─ 0104─ 03
1 引言
汽车在转弯制动联合工况下,每个轮胎的载荷 分配都要受到转向和制动所产生的加速度以及车身 俯仰和侧倾的影响,会发生很大的变动,从而导致 各 轮 垂 直 载 荷 的 不 一 致 。从 轮 胎 附 着 特 性 分 析 可 知 , 轮胎上的垂直载荷分配的不同会造成各个车轮所受 到的地面附着力的差异,对汽车的操纵稳定性产生 很 大 的 影 响 [1]。 在 研 究 整 车 系 统 运 动 之 前 , 必 须 详 细分析各个车轮在运动过程中的载荷分配。影响轮 胎的垂直载荷的因素非常多,全部考虑起来是极其 复杂的,为使问题简化,作如下假设:
⎪⎪Fz2 = (Wcb+ Fsxhs +Msx + Ffxhf + Frxhr) 2L+
⎪⎪ (Fsyfhroll +Cεfxε + Ffyhf ) Bf ⎪⎨Fz3 =[Wca−(Fsxhs +Msx + Ffxhf + Frxhr)] 2L−
⎪ ⎪
(Fsyrhroll +Cεrxε + Fryhr) Br
⎪⎪Fz4 =[Wca−(Fsxhs +Msx + Ffxhf + Frxhr)] 2L+
⎪⎩ (Fsyrhroll +Cεrxε + Fryhr) Br
(6)
车轮上的载荷不能小于零,所以当车轮载荷小
于 等 于 零 时 用 零 代 替 。由 于 所 模 拟 的 车 辆 为 左 转 弯 ,
受到惯性离心力的作用,必然是左侧轮胎垂直载荷
纵向惯性力引起的左、右后轮载荷变动量为
偏离产生的附加力矩与向侧倾轴移动产生的惯性力 矩 的 合 力 矩 在 后 轴 的 分 量 ; Fryhr 为 后 非 悬 挂 质 量 的
∆z3x = ∆z4x = −∆z1x 。
惯 性 力 矩 , Fryhr = −mrhr (ay − br&) 。