轮胎模型计算的Bond图方法
- 格式:pdf
- 大小:239.83 KB
- 文档页数:3
1前言Bond图是一种功率流程图[1],亦称功率键合图,可通过描述功率在系统中的传递过程来建立系统的模型和计算过程。
它是一种模块化方法,利用基本模块可作系统模型的扩展,亦可作模型环节的细化,其过程清晰,形式规范。
因此,在建立系统Bond图时,无需事先列出完整的系统微分方程;而且在根据Bond图进行仿真计算时,每一步的计算式不超过一阶方程,并可以使用非线性算式。
对于线性或线性化处理的系统,Bond图又很容易转化为传递函数或状态方程,从而得到简约形的系统模型。
本文通过Bond图建立了车辆轮胎的计算模型,计算结果与转鼓试验台的试验结果吻合较好,说明Bond图的轮胎计算模型是有效的。
2轮胎模型2.1轮胎物理模型轮胎的物理模型如图1所示,它包括轮辋、轮胎胎壁、带束、胎冠及滑移模型。
在车轮高速运动中,带束可以假设为一个刚性带环,其半径为轮胎有效滚动半径r e。
带束与轮辋通过胎壁的弹性和阻尼相联系。
轮轴上的驱动转矩经过胎壁传递至带束,并通过胎冠的弹性变形和滑移产生地面对车轮的驱动力。
胎冠变形和滑移与路面作用力的关系,在图1表达为阻尼与弹簧相串联的模型,但这是高度非线性的阻尼器与弹簧,只有当滑移率在原点附近时,才表现出较好的线性关系。
图1轮胎物理模型1.轮辋2.带束3.路面4.滑移模型5.胎冠模型 6.胎壁径向刚度和阻尼7.胎壁扭转刚度和阻尼8.胎壁2.2轮胎模型的Bond图Bond图用汇点代表分流功率的环节,在汇点的分流路线上使用了两个变量表示功率(如力和速度)。
每个汇点可带有阻抗、刚度和惯量等参数,代表功率在该环节上的消耗。
Bond图的汇点有0型和1型两种。
在0型汇点中,各路线上的力相等,而速度之和为零;1型汇点则相反,各路线上的速度相等,但力之和为零,两者都反映了各环节中输入输出能量的守恒。
Bond图用半箭头代表功率的流向,箭头端部短线表示力参数与速度参数是前因后果的关系,即ω=f(M);箭头尾部短线则表示相反的因果关系,即M=f(ω)。
轮胎稳态模型的分析综述_张向文(2)4轮胎经验模型轮胎经验模型是直接根据试验测试数据拟合得到的模型,与试验结果较接近,而经验模型公式简单,便于计算和实际应用,但需要大量的试验数据。
由于试验条件限制和路面状况的多变性,难以得到所有路面状况和所有轮胎运动状态的试验数据。
因此,经验模型只是根据有限的试验数据得到,模型外推性不好,参数没有明确的物理意义。
4.1多项式模型多项式模型由S.Germann 等人提出,其利用简单的多项式函数近似描述轮胎与路面摩擦系数和滑移率之间的关系[42]:μ=a 0+a 1s +a 2s 2(50)式中,参数a 0、a 1和a 2需要通过试验数据进行辨识。
根据辨识的参数,利用多项式模型可以方便的进行摩擦系数求解和汽车控制系统设计,但是该模型仅在滑移率较小时误差较小,当滑移率逐渐增大时,误差会越来越大。
4.2Burckhardt 模型Burckhardt 模型是M.Burckhardt 提出的一种摩擦系数μ与滑移率s 的关系模型[43,44]:μ(s )={c 1[1-exp (-c 2s )]-c 3s }e-c 4v(51)式中,c i (i =1,…,4)随路面状况的变化而变化,可以通过试验测试数据拟合得到;e -c 4v反映速度变化引起的摩擦系数变化。
若忽略速度变化影响,Burckhardt 模型可以简化为[43~47]:μ(s )=c 1[1-exp (-c 2s )]-c 3s(52)根据简化模型,利用试验测试数据可以拟合得到不同路面状况下的参数如表1所列。
表1不同路面状况下Burckhardt 模型各参数的典型值为了分析Burckhardt 模型特性,利用式(52)和式(9)、式(10)进行仿真研究。
利用表1的参数仿真不同路面状况下纵向、侧向摩擦系数随纵向滑移率和侧偏角的变化如图27所示,仿真中α=8°,s x =0.2。
(a )随纵向滑移率的变化(b )随侧偏角的变化图27Burckhardt 模型不同路面状况下的纵向摩擦系数和侧向摩擦系数轮胎稳态模型的分析综述*张向文1王飞跃2高彦臣3(1.桂林电子科技大学;2.中国科学院自动化研究所复杂系统管理与控制国家重点实验室;3.软控股份有限公司)觹基金项目:国家自然科学基金项目(60804059);广西自然科学基金项目(2010GXNSFA013130);中国科学院复杂系统与智能科学重点实验室开放课题。
第21卷第4期2007年7月山东理工大学学报(自然科学版)JournalofShandongUniversityofTechnology(NaturalScienceEdition)V01.21No.4Jul.2007文章编号;16726197(2007)040100—03基于模型的轮胎参数估计方法孙春玲,韩加蓬,张锋(山东理S-大学变通与车辆S-程学院,山东淄博255049)摘要:利用汽车ABS轮速传感器信号,在建立轮胎滚动的纵向力模型的基础上,运用最小二乘法等算法对轮胎半径和纵向刚度进托估计,并通过线性仿真进行进一步验证,由此来监测轮胎压力的变化情况,为间接TPMS的研究提供了一种较好的方法.关键词:轮胎压力;监测系统;轮胎模型;滑移率;轮胎半径中图分类号:U463.341文献标识码:AEstimationoftireparametersbasedonamodelSUNChun—ling.HANJia-peng.ZHANGFeng(SchoolofTrafficandVehicleEngineering.ShandongUniversityofTechnology·Zibo250049,China)Abstract:ThispapermakesuseofsignalsoftheABSwheelspeedsensors,estimatesthetireradialandlongitudinalstiffnesstimelybyleast—squearsalgorithmcombiningtheestablishingofthe10ngitudinalforcemodelofrollingtire,andmakesfurthervalidationbylinearsimula—tionformoniteringthechangeoftirepressureandprovidingagoodmethodfortheresearchoftheindirectTPMS.Keywords:tirepressure;monitoringsystem;tiremodel;slip;tireradius近几年来,由于轮胎问题而引发的重大交通事故数目呈逐年上升趋势.据报道,在全国众多车辆事故中,因轮胎引发的事故约占20.8%,高速公路每百公里事故率为普通公路的4倍多,其中70%的是由于轮胎爆破而引起的,而在美国这一比例则达到80%.轮胎故障问题造成了巨大损失,怎样防止爆胎已成为安全驾驶的一个重要问题.因此,汽车在行驶中,实时在线、不断地监测轮胎气压是否正常对保证行驶安全至关重要.世界各国学者都在努力找到监测轮胎压力的好方法.目前,轮胎压力监测系统(tirepressuremonito收稿日期:20061015作者简介:孙春玲(1981),女,硬士研究生ringsystem,TPMS)主要分为直接型TPMS和间接型TPMS两种类型.直接型TPMS利用安装在轮胎里的压力传感器来直接测量轮胎的气压,但是需要增加新的硬件,成本较高.间接型TPMS利用汽车现有的传感器采集所需信号,并通过软件算法监测轮胎压力,不需要增加新的硬件,安装费用较低.本文(假设汽车是前轮驱动)利用ABS轮速传感器信号,在建立轮胎纵向力模型的基础上,运用最dx--乘法等算法对轮胎半径和纵向刚度进行实时估计,并通过线性仿真进一步验证并得到较好的结果,通过分析二者的变化来监测轮胎压力的变化 万方数据第4期孙春玲,等:基于模型的轮胎参数估计方法101情况,作为轮胎压力监控系统的理论依据,为间接TPMS的研究提供了一种较好的方法.1轮胎模型的建立对于自由滚动的轮胎,其角速度和速度之间满足如下关系o3V=R·训式中:V为轮胎中心的速度;R为轮胎的滚动半径;W为轮胎的角速度.然而,当车辆行驶的时候,由于轮胎和地面之间的作用力,轮胎始终存在着滑移.轮胎滑移率o]可以定义为S一一—V—--iR-.w(1)式中s为车轮滑移率.魔术公式01是根据车轮滑移率预测轮胎纵向力的典型轮胎模型,它来源于经验数据,在限定的载荷、路面、轮胎特性和其他因素条件下,将车轮滑移率与轮胎纵向力联系在一起.从不同摩擦系数路面上魔术公式轮胎模型的典型力滑移率曲线中(见图1)可以看出:当轮胎滑移率在3%以内时,力和滑移率两者之间是近似的线性关系.而汽车在正常行驶时,轮胎滑移率很少超过1%~2%.因此,在这一线性区域,两者之间的关系可以定义为F—C,S(2)式中:F为轮胎传递的纵向力;C;为轮胎的纵向刚度.裟4∞。
为方便具有轮胎非线性分析需求的用户熟悉Marc(Marc Mentat)中轮胎建模的方法和流程,针对某汽车轮胎的装配、充气、在路面承受后轮胎的变形和应力分析进行描述。
一、使用Marc Mentat建立轮胎二维轴对称模型模拟装配到轮辋和充气过程首先根据轮胎的结构和尺寸参数在Mentat中建立下图所示的有限元模型(195/65R15汽车轮胎),用户可以使用Marc Mentat直接创建轮胎截面的几何模型和有限元模型,也可以利用Mentat提供的接口,将其他CAD或CAE软件创建的模型导入到Mentat中进行后续材料参数、边界条件、分析参数等的定义。
Mentat提供了多种商用的CAD和CAE软件的接口,具体可以参考Marc用户手册的介绍。
本例出示的轮胎模型包括橡胶胎面、带束层、胎冠(tread、base、rubber)以及布帘等加强结构(bead、rebar1、rebar2)。
如下图所示:轮胎截面有限元模型(二维轴对称)橡胶材料部分可以采用Marc提供的Mooney模型定义,根据实际材料特性输入相应的材料参数即可。
Marc Mentat提供了多种模拟橡胶材料的本构模型和实验曲线拟合工具,用户可以根据供应商提供的或实测的该橡胶材料的实验曲线(应力-应变曲线)选择合适的超弹性材料模型进行拟合,并由Mentat自动计算和应用材料参数到模型中。
具体步骤可参考Marc用户手册或基础培训教程中的相关介绍。
详细内容可参考mar103教程中的介绍.轮胎各部分材料类型分布对于加强材料,这里包括了两类,一类是金属圈bead结构,直接采用各向同性材料本构模型,输入相应的结构材料参数,例如弹性模量、泊松比等即可。
另一类加强筋材料采用嵌入式模型(本例中加强筋单元嵌入到rubber基体材料中),用于模拟轮胎橡胶材料中嵌入的布帘和加强筋结构,这些结构可以指定为Marc中的rebar单元来模拟。
Marc支持一种基体材料中同时嵌入多层和多种加强筋材料的定义,这些加强筋结构可以分层分布在基体材料不同的厚度处、加强筋的铺设方向、截面积以及数量等均可以根据实际结构定义。
1前言Bond图是一种功率流程图[1],亦称功率键合图,可通过描述功率在系统中的传递过程来建立系统的模型和计算过程。
它是一种模块化方法,利用基本模块可作系统模型的扩展,亦可作模型环节的细化,其过程清晰,形式规范。
因此,在建立系统Bond图时,无需事先列出完整的系统微分方程;而且在根据Bond图进行仿真计算时,每一步的计算式不超过一阶方程,并可以使用非线性算式。
对于线性或线性化处理的系统,Bond图又很容易转化为传递函数或状态方程,从而得到简约形的系统模型。
本文通过Bond图建立了车辆轮胎的计算模型,计算结果与转鼓试验台的试验结果吻合较好,说明Bond图的轮胎计算模型是有效的。
2轮胎模型2.1轮胎物理模型轮胎的物理模型如图1所示,它包括轮辋、轮胎胎壁、带束、胎冠及滑移模型。
在车轮高速运动中,带束可以假设为一个刚性带环,其半径为轮胎有效滚动半径r e。
带束与轮辋通过胎壁的弹性和阻尼相联系。
轮轴上的驱动转矩经过胎壁传递至带束,并通过胎冠的弹性变形和滑移产生地面对车轮的驱动力。
胎冠变形和滑移与路面作用力的关系,在图1表达为阻尼与弹簧相串联的模型,但这是高度非线性的阻尼器与弹簧,只有当滑移率在原点附近时,才表现出较好的线性关系。
图1轮胎物理模型1.轮辋2.带束3.路面4.滑移模型5.胎冠模型 6.胎壁径向刚度和阻尼7.胎壁扭转刚度和阻尼8.胎壁2.2轮胎模型的Bond图Bond图用汇点代表分流功率的环节,在汇点的分流路线上使用了两个变量表示功率(如力和速度)。
每个汇点可带有阻抗、刚度和惯量等参数,代表功率在该环节上的消耗。
Bond图的汇点有0型和1型两种。
在0型汇点中,各路线上的力相等,而速度之和为零;1型汇点则相反,各路线上的速度相等,但力之和为零,两者都反映了各环节中输入输出能量的守恒。
Bond图用半箭头代表功率的流向,箭头端部短线表示力参数与速度参数是前因后果的关系,即ω=f(M);箭头尾部短线则表示相反的因果关系,即M=f(ω)。
此外,Bond图用TF表示动力参数的转换环节,如转矩和转速通过半径为r的滚筒转换为力和线速度。
从驱动功率输入车轮轮轴到地面切向推力驱动车轮行走的功率传递过程,可以得到图2所示的Bond图。
图中M为转矩,ω为转速,R为阻抗,K 为刚度,I为惯量,F z为法向力,F x为纵向力,V x 为速度,V bx为滑移速度,r e为有效滚动半径,R a为轮胎模型计算的Bo nd图方法同济大学徐志新周钅宏【Abstract】The Bond g ra p h is a kind of modularized p ower p rocedure chart,which establishes model and calculatin g p rocedure of a s y stem b y descri p tion of t he transmittin g p rocess of p ower in t he s y stem.The Bond g ra p h of motor vehicle t y re model is established and calculation p rocedure is g iven out accordin g to t he g ra p h.Result of calculation coincides better wit h test result on chassis d y namome2 ter and demonstrates t he effectiveness of t he Bond g ra p hic calculation met hod of t y re model.【摘要】Bond图是一种模块化的功率流程图,它通过描述功率在系统中的传递过程,建立系统模型及计算过程。
建立了车辆轮胎模型的Bond图,并根据该图给出了计算流程。
计算结果与转鼓试验结果吻合较好,说明轮胎模型的Bond图计算方法是有效的。
To p ic words:T y re,Model,Calculation method主题词:轮胎模型计算方法滚动阻力系数;下角标符号中r 代表轮辋,s 代表胎壁,b 代表带束,c 代表轮胎着地,x 代表轮心运动。
图2轮胎模型的Bond 图在图2中,因轮辋是刚体,各处的转速ωr 相同,故为1型汇点,输入的转矩M r 在克服轮辋转动惯性I r 及阻尼R r 后传入胎壁;胎壁允许弹性变形,输入与输出的速度不同,所以为0型汇点,它传递的转矩由胎壁角位移刚度K s 及阻尼R s 决定;带束设为刚环,可作为1型汇点,输入的转矩除克服轮胎转动惯性I b 外,还需克服滚动阻力矩M R a 及阻尼R b ;其后输出的转矩及转速ωb ,经轮胎有效滚动半径r e 转换为带束的周向力和线速度V b ;带束线速度V b 与轮心速度V x 的差V bx 就是滑移速度,根据滑移特性,这将产生路面切向推力F x ,所以轮胎与路面的接触是0型汇点;该切向推力作用于轮心,克服轮心运动惯性I x 、负载和运动阻尼R x ,驱动车轮行走,因此轮心运动为1型汇点。
根据图2所示的Bond 图模型可得到轮胎模型的计算框图(见图3)。
在计算中,角速度ω、弹性转矩M K 及阻尼转矩M R 可由式(1)~式(3)得到。
ω=1I ΘM Id t +ω(0)(1)M K=K Θωd t +M K(0)(2)M R =R ω(3)速度Vx 、阻尼力F R 的计算与此类似。
在本文中,滑移率模型采用线性Brush 公式(刷子模型)[2]:F X =2a 2C Pλ(4)式中,λ为轮胎滑移率;C P 为胎冠刚度;a 为轮胎半接地长度。
由图3可清楚地看出轮轴驱动力矩转换为车轮驱动力的计算过程。
在该过程中,各环节位置清楚,模块分解细致,计算关系明了,各参数对仿真结果的影响易于考察。
从车轮动力传递的实际过程来看,Bond 图的计算模型是合理的。
图3轮胎模型计算框图3Bond 图模型计算及试验结果3.1模型参数提取及主要试验参数平面内运动的轮胎模型的主要特征参数有径向刚度、转动刚度、纵向滑移刚度、有效滚动半径及轮胎/车轮惯量等。
轮胎径向刚度通过对轮胎垂直加载,压向一平面,由测得的轴荷及轮轴垂直位移得到;轮胎转动刚度主要由胎壁角位移刚度决定,将轮胎套在一与胎体内径相同的钢环上,在保证轮胎没有径向位移及胎体保持圆弧状的情况下,通过加载并测量轮辋的相对角位移而得到;纵向滑移刚度由实测的稳态滑移率特性曲线在原点处的斜率得到。
表1列出了主要试验参数及轮胎参数。
3.2仿真及试验结果为了考察模型的计算效果,在转鼓试验台上进行了轮胎制动试验(轮胎规格为185/70R14)。
首先表1主要试验参数及轮胎参数参数数值轮轴静载/N6000转鼓速度/km ・h-120最大制动力矩/N ・m1240制动时间/s8轮胎/轮轴转动惯量/k g ・m20.636/0.627轮胎有效滚动半径/m 0.3021轮胎接地长度/m0.137胎冠纵向刚度/N ・m -21.107轮胎转动刚度/N ・m ・rad -17.7×104轮胎转动阻尼/N ・m ・s ・rad-150轮胎径向刚度/N ・m -12.71×105轮胎径向阻尼/N ・s ・m-1320(下转第41页)71104117×107的差距呢?笔者认为:a 1原有计划经济的原因。
在旧的经营模式下,企业领导干好干坏一个样,而采用CAD/CAM 技术投资大,有较大风险,效益回报有一定的滞后期,导致企业领导不够重视,怕影响正常的生产。
b 1管理上的缺陷。
由于条块分割,重复引进,企业相互之间缺乏必要的交流和协作,影响了CAD/CAM 技术效益的发挥。
c 1人才培养的不足。
像福特、通用等汽车公司之所以CAD/CAM 技术应用得好,是因为几十年来一直大力开展CAD/CAM 应用而积淀下来的宝贵经验和培养出了一支高水平的技术队伍。
有钱可以买到先进的软、硬件,但一支高水平的技术队伍需要若干年的培养。
现在国内既懂计算机软、硬件,又有丰富专业知识的人才奇缺,而这恰是企业最为宝贵的财富。
随着科学技术的迅猛发展,CAD/CAM 技术正向集成化、智能化、网络化、标准化、可视化等方向发展。
新的技术正在不断出现,在这新旧交替之际,只要我们抓住机遇,就能跳跃前进,赶上世界科技发展的步伐;反之,如不能把握机会,则将更加落后。
可以相信,中国的汽车CAD/CAM 事业前途光明,但需要大家付出艰辛的努力。
(责任编辑郝旭辉)原稿收到日期为1998年10月15日。
(上接第11页)在轮轴上施加6000N 恒定垂直载荷,再由215m 转鼓以20km/h 的线速度带动车轮转动。
通过钳式制动器对轮胎进行制动力矩加载,制动器的制动压力由一液压伺服系统控制。
为了减少轮胎瞬态特性对制动过程的影响,采用了慢制动试验,制动压力逐步增加再逐步降低,全部制动时间为8s 。
将制动过程中轮胎与转鼓的转速、轮胎滑移速度和轮胎制动力等实测数据记录下来,然后将实测的制动力矩送入Bond 图计算模型,作为计算模型的输入,因此模型的仿真结果与实测结果是通过相同的输入得到的。
模型仿真及试验结果见图4(图中实线为仿真结果,虚线为试验结果)。
图4模型仿真与台架试验计算结果与试验结果的对比表明,两者在量值及变化趋势方面都能较好的吻合,说明轮胎模型的Bond 图计算方法是相当有效的。
计算与试验结果的差别主要来自于两个方面,一是计算时采用了线性Brush 模型,该模型描述了零点附近滑移率的稳态特性,而实际试验的结果总带有瞬态和非线性特性;二是轮胎的有关参数存在一定的测量误差。
若采用动态的轮胎模型,如根据轮胎松弛长度变化建立的一阶滑移率模型[3]等,计算结果会有进一步改善。
Bond 图模块化结构的优点正是在于,当进一步细化模型或改变滑移率模型以便用于不同行驶条件的仿真时,整个计算模块无需大的改动。
4结论a 1用Bond 图方法进行轮胎模型计算,具有计算简单、调整方便及过程详尽的特点。
计算结果与实测试验结果的比较表明,该计算方法是有效的。
b 1模块化结构的Bond 图计算方法,对于轮胎模型细化、扩展及变化均有较好的适应性。
参考文献1Karno pp D ,Rosenber g R 1S y stem D y namics 1Wile y :A U 2nified A pp roach ,197512Pace j ka H B 1Modellin g of The Pneumatic Tire and It s Im p act on Vehicle D y namic Behaviour 1Ober p faffenhofen :Carl Cranz Gesellschaft ,19851(责任编辑路人)修改稿收到日期为1998年12月9日。