七自由度车辆数学模型
七自由度车辆图和数学模型
图一.七自由度车辆动力学模型纵向力平衡方程:121234()()cos ()sin x y x x y y x x m V r V F F F F F F δδ∙-⋅=+-+++侧向力平衡方程:121234()()sin ()cos y x x x y y y y m V r V F F F F F F δδ∙+⋅=+++++绕Z 轴力矩平衡方程:11212211224334[()sin ()cos ][()cos ()sin ]2()()2w z x x y y x x y y w x x y y t I r F F F F a F F F F t F F F F bδδδδ∙⋅=++++-+-+--+四个车轮的力矩平衡方程:tw w xi bi di I R F T T iw ∙⋅=-⋅-+上述方程中:δ为前轮转角;Vx ,Vy 分别为纵向、横向车速;β为质心侧偏角;γ为横摆角速度;Fxi 、Fyi 、Fzi 分别为轮胎纵向力、侧向力、垂向力;i=1、2、3、4,为分别对应的车轮;m 为整车质量;ms 为悬挂质量;a 、b 为前后轴到质心的距离;l=a+b 为前后轴距;tw1为前轴轮距;tw2为后轴轮距;d= 122tw tw +为平均轮距;Iz 为整车绕Z 轴的转动惯量;h 为质心到地面的距离;纵向加速度为x x y a V r V ∙=-⋅ 侧向加速度为 y y x a V r V ∙=+⋅各轮胎垂向载荷公式:123422222222z s x s y z s x s y z s x s y z s x s y b h h bF mg m a m a l l d l b h h bF mg m a m a l l d l a h h aF mg m a m a l l d l a h h aF mg m a m a l l d l=--⋅=-+⋅=+-⋅=++⋅各轮胎侧偏角公式:112212223242arctan()arctan()arctan()arctan()w w w w t t t t Vy arVx r Vy arVx r Vy brVx r Vy brVx rαδαδαα+=--+=-+-=---=-+各车轮轮心在车轮坐标系下的纵向速度:11122324()cos ()sin 2()cos ()sin 222w t w t w t w t t V Vx r Vy ar t V Vx r Vy ar t V Vx rt V Vx rδδδδ=-++=+++=-=+各车轮滑移率的计算:1111222233334444w t t w t t w t t w t t w R V V w R V V w R V V w R V V λλλλ-=-=-=-=轮胎模型:Dugoff 轮胎模型纵向力:()1x ixiiC F f s λλ=-侧向力:tan ()1y i yiiC F f s αλ=-(1)i S λ=-11()(2)1S f S S S S >⎧=⎨-<⎩。
图一.七自由度车辆动力学模型
图一.七自由度车辆动力学模型纵向力平衡方程:()()cos ()sin x y xfl xfr yfl yfr xrl xrr m V r V F F F F F F δδ•-⋅=+-+++ ①侧向力平衡方程:()()sin ()cos y x xfl xfr yfl yfr yrl yrr m V r V F F F F F F δδ•+⋅=+++++ ②绕Z 轴力矩平衡方程:12[()sin ()cos ][()cos ()sin ]2()()2w z xfl xfr yfl yfr xfr xfl yfl yfr w xrr xrl yrl yrr t I r F F F F a F F F F t F F F F b δδδδ•⋅=++++-+-+--+ ③四个车轮的力矩平衡方程:tw w xi bi di I R F T T i w •⋅=-⋅-+ ④上述方程中:δ为前轮转角;Vx ,Vy 分别为纵向、横向车速;β为质心侧偏角;γ为横摆角速度;Fxi 、Fyi 、Fzi 分别为轮胎纵向力、侧向力、垂向力;i=左前轮fl 、右前轮fr 、左后轮rl 、右后轮rr ,为分别对应的车轮;m 为整车质量;ms 为悬挂质量;a 、b 为前后轴到质心的距离;l=a+b 为前后轴距;tw1为前轴轮距;tw2为后轴轮距; Iz 为整车绕Z 轴的转动惯量;hg 为质心到地面的距离;各轮胎垂向载荷公式:_1_1_1_122222222gg z fl x y w gg z fr x y w gg z rl x y w gg z rr x y w h h b b F mg mV mV l l t l h h b b F mg mV mV l l t lh h aa F mg mV mV l l t lh h a a F mg mV mV l l t l••••••••=--⋅=-+⋅=+-⋅=++⋅⑤各轮胎侧偏角公式:11222222arctan()arctan()arctan()arctan()wwww y fl t x y fr t x y rl t x y rr t x V arV rV arV rV brV rV br V rαδαδαα+=--+=-+-=---=-+ ⑥各车轮轮心在车轮坐标系下的纵向速度:1_1_2_2_()cos ()sin 2()cos ()sin 222w t fl x y w t fr x y w t rl x w t rr x t V V r V ar tV V r V ar t V V rtV V rδδδδ=-++=+++=-=+ ⑦其中_t fl V 、_t fr V 、_t rl V 、_t rr V 为轮胎坐标系下的轮胎纵向速度。
微型轿车七自由度模型的仿真分析与试验研究
关键词: 汽车; 平顺性; 仿真; 试验
中图分类号: U 462 2+ 4
文献标识码: A
Study on Sim ulation and Test for a Seven D O F M odel of a M ini Car
L IANG X in cheng, ZHANG Jun, DA I X in, ZHOU F eng jun ( N ational Eng ineering L aboratory for E lectric Vehicle, Beijing Inst itute o f T echno logy, Be ijing 100081, Ch ina)
K ey w ord s: Car; R ide com for;t S im ulat ion; T est
在汽车的各种性能中, 平顺性是一个比较重要 的评价指标. 平顺性较差会缩短元件的寿命, 还会 造成乘员工作效率低, 影响人的身心健康. 长期处 于不舒适的振动环境中容易引发各种心脏疾病. 汽 车的振源主要有路面不平度激励、发动机激励、传 动轴不平衡激励、轮胎激励、侧向风激励等. 路面 不平度对汽车产生的激励具有随机性, 其统计特性 主要采用路面功率谱密度来描述. 通常汽车模型参 数越多, 与实际越吻合. 但参数过多, 模型求解异
1 2
kfA
(Z1
-
q1 ) 2
+
1 2
kfB
(Z2
-
q2 ) 2
+
1 2
kfC
(
Z3
-
q3 ) 2
+
1 2
kfD
(Z4
-
q4 ) 2.
( 5)
取广义坐标 { u }T = [ Z1、 Z 2、 Z3、 Z4、 Zb、 b、
基于7自由度模型的整车振动分析
基于7自由度模型的整车振动分析整车振动分析是一种重要的工程技术手段,用于研究整车在正常行驶过程中的振动性能。
在整车振动分析中,7自由度模型是一种经典的方法。
7自由度模型是指整车被简化为一个由7个自由度构成的系统,各个自由度分别表示车体的平动、横摆、侧倾、纵向加速度、纵向速度、纵向位移和纵向倾斜。
这种简化模型可以提供足够的信息,以便评估整车在各种工况下的振动特性。
整车在振动分析中的振动特性主要包括固有频率、振型和振幅。
固有频率是指整车固有振动的频率,是整车振动特性的一个重要指标。
振型是指在特定频率下,整车各个自由度的振动形态。
振幅是指振动的幅度或位移,反映了整车振动的强度。
整车振动分析的过程一般包括以下几个步骤:第一步是建立整车的7自由度模型。
这一步需要将整车分解为各个自由度,并确定各自由度之间的耦合关系。
第二步是确定整车的质量分布和各种刚度参数。
这些参数可以通过实际测试和数值模拟得到。
第三步是求解整车的特征值问题,得到整车的固有频率和振型。
求解特征值问题可以采用有限元方法或其他数值求解方法。
第四步是分析整车在各种工况下的振动响应。
根据不同的激励条件(如路面不平度、加速度等),可以计算出整车各个自由度的振动响应。
第五步是评估整车的振动性能。
通过对振动响应进行分析和比较,可以评估整车在各种工况下的振动性能,并对设计进行优化。
整车振动分析的结果对车辆的安全性、舒适性和稳定性有着重要的影响。
一方面,合理的振动控制可以提高车辆的抗振性能,减少碰撞风险。
另一方面,舒适的振动特性可以提高驾驶者和乘客的乘坐体验,提高车辆的市场竞争力。
总的来说,基于7自由度模型的整车振动分析是一种重要的工程技术手段,可以帮助工程师评估和优化整车的振动性能,提高车辆的安全性、舒适性和稳定性。
通过适当的振动控制和优化设计,可以为用户提供更好的驾驶体验,并提高整车的市场竞争力。
基于七自由度整车模型的汽车振动特性分析
d o p t e d t o s i mu l a t e t h e d y n a mi c c h a r a c t e r i s t i c s o f he t s y s t e m a n d he t i n l f u e n c e r u l e s o f s o me p a r a me t e r s o f t h e V C — h i e l e ,s u c h a s d a mp a n d s t i f f n e s s ,o n he t d y n a mi c c h a r a c t e is r t i c s o f v e h i c l e a r e s t u d i e d .T h e r e s u l t s r e v e a l ha t t
王
( 安徽科 技学 院 摘
娟, 李 同杰 , 孟令启
凤阳 2 3 3 1 0 0 )
机 电与车 辆工程 学 院 , 安徽
要: 本文建 立 了独立 悬 架汽车整 车七 自由度动 力 学模 型 , 推 导 了振 动 微 分 方程 , 运 用 数值 仿 真 的 手段
研究了汽车悬架阻尼 系 数、 刚度等设计参数对汽车振动特性的影响规律 。结果表明 , 在共振 区域附近, 相
WA N G J u a n , L I T o n g — j i e , ME N G L i n g— q i
( C o l l e g e o f M e c h a r t r o n i c s a n d V e h i c l e E n g i n e e i r n g , A n h u i S c i e n c e nd a T e c h n o l o g y U n i v e r s i t y , F e n g y a n g 2 3 3 1 0 0 , C h i n a )
7自由度
1.1 车辆动力学模型
建立七自由度车辆模型如图2所示,假定车身为一刚性体,车轮简化为质点,悬架类型可简化为四个独立悬架,七个自由度分别是车身的垂向位移、侧倾角、俯仰角和四个车轮的垂向位移。
建立纵向力平衡方程:
()()cos ()sin x y xfl xfr yfl yfr xrl xrr m V r V F F F F F F δδ∙
-⋅=+-+++(1)
建立侧向力平衡方程:
()()sin ()cos y x xfl xfr yfl yfr yrl yrr m V r V F F F F F F δδ∙
+⋅=+++++(2)
建立绕Z 轴力矩平衡方程:
12[()sin ()cos ][()cos ()sin ]
()()22
z xfl xfr yfl yfr xfr xfl w w yfl yfr xrr xrl yrl yrr I r F F F F a F F l l
F F F F F F b δδδδ∙
⋅=++++-+-+--+(3) 建立四个车轮的运动微分方程:
tw i w xi bi di I w R F T T ∙
⋅=-⋅-+
式中:δ为前轮转角;β为质心侧偏角;γ为横摆角速度;V x ,V y 分别为纵向、横向车速;F x 、F y 、F z 分别为轮胎纵向力、侧向力、垂向力;fl 、fr 、rl 、rr 分别对应左前轮、右前轮、左后轮、右后轮;m 为整车质量;a 、b 为前后轴到质心的距离;l=a+b 为前后轴距;l w1为前轴轮距;l w2为后轴轮距;I Z 为绕Z 轴的转动惯量。
Ms 为悬挂质量,tw1为前轴轮矩,tw2为后轴轮矩。
车辆系统动力学结构模型
mb2
Ib2
kp
mw4
zw4
cp
z b2 b2
mw3
zw3
I hp
m 2lc
oc
hp
zc c
cs
kp
mb1
Ib1
cp
zb1 b1
mw2
zw2
mw1
zw1
29
车辆部件受力分析
Fs(2)
Mc. g
Fs (1)
Fp(4) Mb. g Fp(3)
Fp(2) Mb. g Fp(1)
Fw(4)
Fw(3)
Fw(2)
Mb Kp Cp zb
Mw
zw
19
三自由度系统振动方程
(7.3)
M c zc cs (zs zb ) k s (zs zb ) 0
M b zb cs (zs zb ) k s (zs zb ) c p (zb zw ) k p (zb zw ) 0
M w zw c p (zb zw ) k p (zb zw ) 0
Fs
Fg
mx Fg Fs
Fg mg
mx cx kx mg
Fc kx cx
x (Fg Fs ) / m
xn1
xn1
xn xn
xnt (1/ 2 )xnt2 (1 )xnt xn1t
xn1t 2
(2)
17
m 1000kg, k 108 N/m, c 104 N s/m, t 10-4, 0.5
26 22 30 34
17 18
33
29 21 25
11,12 3,4
9,10 1,2
1-8 轮轨力 17-20 中央悬挂力 25-28 抗蛇行减振器阻尼力 33-36 牵引拉杆力
汽车2自由度和7自由度动力学建模仿真
汽车2自由度和7自由度动力学建模仿真案场各岗位服务流程销售大厅服务岗:1、销售大厅服务岗岗位职责:1)为来访客户提供全程的休息区域及饮品;2)保持销售区域台面整洁;3)及时补足销售大厅物资,如糖果或杂志等;4)收集客户意见、建议及现场问题点;2、销售大厅服务岗工作及服务流程阶段工作及服务流程班前阶段1)自检仪容仪表以饱满的精神面貌进入工作区域2)检查使用工具及销售大厅物资情况,异常情况及时登记并报告上级。
班中工作程序服务流程行为规范迎接指引递阅资料上饮品(糕点)添加茶水工作要求1)眼神关注客人,当客人距3米距离时,应主动跨出自己的位置迎宾,然后侯客迎询问客户送客户注意事项15度鞠躬微笑问候:“您好!欢迎光临!”2)在客人前方1-2米距离领位,指引请客人向休息区,在客人入座后问客人对座位是否满意:“您好!请问坐这儿可以吗?”得到同意后为客人拉椅入座“好的,请入座!”3)若客人无置业顾问陪同,可询问:请问您有专属的置业顾问吗?,为客人取阅项目资料,并礼貌的告知请客人稍等,置业顾问会很快过来介绍,同时请置业顾问关注该客人;4)问候的起始语应为“先生-小姐-女士早上好,这里是XX销售中心,这边请”5)问候时间段为8:30-11:30 早上好11:30-14:30 中午好 14:30-18:00下午好6)关注客人物品,如物品较多,则主动询问是否需要帮助(如拾到物品须两名人员在场方能打开,提示客人注意贵重物品);7)在满座位的情况下,须先向客人致歉,在请其到沙盘区进行观摩稍作等待;阶段工作及服务流程班中工作程序工作要求注意事项饮料(糕点服务)1)在所有饮料(糕点)服务中必须使用托盘;2)所有饮料服务均已“对不起,打扰一下,请问您需要什么饮品”为起始;3)服务方向:从客人的右面服务;4)当客人的饮料杯中只剩三分之一时,必须询问客人是否需要再添一杯,在二次服务中特别注意瓶口绝对不可以与客人使用的杯子接触;5)在客人再次需要饮料时必须更换杯子;下班程序1)检查使用的工具及销售案场物资情况,异常情况及时记录并报告上级领导;2)填写物资领用申请表并整理客户意见;3)参加班后总结会;4)积极配合销售人员的接待工作,如果下班时间已经到,必须待客人离开后下班;1.3.3.3吧台服务岗1.3.3.3.1吧台服务岗岗位职责1)为来访的客人提供全程的休息及饮品服务;2)保持吧台区域的整洁;3)饮品使用的器皿必须消毒;4)及时补充吧台物资;5)收集客户意见、建议及问题点;1.3.3.3.2吧台服务岗工作及流程阶段工作及服务流程班前阶段1)自检仪容仪表以饱满的精神面貌进入工作区域2)检查使用工具及销售大厅物资情况,异常情况及时登记并报告上级。
汽车理论第七章
单质量系统受力分析
简化得二阶常系数齐次方程
Mz cz kz cq kq
对应的齐次方程为 Mz cz kz 0
令2n c , m
02
K m
f0
0 2
1
2
K m
n c 0 2 Mk
n称为系统的阻尼系数,ω0称为系统的固有频率
z
A e nt
1.悬架结构
• 影响因素 • 1.弹簧刚度K • 2 .悬架弹性特性 • 3.减振器阻尼力系数 C。
2.轮胎
• 提高轮胎缓冲性能的方法有:
• (1)增大轮胎断面、轮轮宽度和空气容量,并相 应降低轮胎气压。
• ( 2)改变轮胎结构型式,如采用子午线轮胎,它因 胎体的径向弹性大,可以缓和不平路面的冲 击并吸 收大部分冲击能量使平顺性得到改善。
• 所以
E[q 2 (t)]
Sq ( f )df
二、路面不平度的自功率谱密度函数——路面谱
1.路面不平度的功率谱通常把相对基准平面的高度q , 沿着道路走向长度l的变化q(l) 称为道路纵(断面)曲线 或不平度函数。
根据测量的路面不平度随机数据,在计算机上处理
得到路面不平度功率谱 Gq (n)
动挠度 f d 对路面输入q的频率响应函数为
fd
2M
2 / 02
q
M2 jck
1 2 02 j 2 0
动挠度 f d 对路面不平度q的幅频特性为
1
fd q
(12
4 )2 4
2
2
2
第四节 “人体—座椅”系 统参数对振动的影响
• ISO2631-1:1997 (E)《人体承受全身振动 评价—第一部份:一般要 求 》
基于旋量理论的7自由度机械臂运动学建模与分析
基金项目:山西省自然科学基金项目(201801D121183)ꎻ山西省重点研发计划项目(201803D421028ꎻ201903D421051)第一作者简介:李冠琦(1996 )ꎬ男ꎬ山西吕梁人ꎬ硕士研究生ꎬ研究方向为康复机器人机构ꎮDOI:10.19344/j.cnki.issn1671-5276.2022.01.027基于旋量理论的7自由度机械臂运动学建模与分析李冠琦ꎬ武建德ꎬ李瑞琴(中北大学机械工程学院ꎬ山西太原030051)摘㊀要:机械臂模仿人手臂的7自由度会拥有冗余自由度ꎮ基于旋量理论计算7自由度机械臂的正向运动学解ꎬ从数值上验证矩阵指数先分块展开比直接泰勒展开准确ꎮ用Newton-Raphson数值迭代法求逆解ꎮ通过编写Matlab程序对正逆解互相验证ꎮ研究发现逆解的求解有适用范围ꎬ较之传统的D-H法ꎬ使用0螺距的螺旋轴ꎬ会使建模更加简洁ꎮ关键词:旋量ꎻ机械臂ꎻNewton-Raphson数值迭代法中图分类号:TP241㊀㊀文献标志码:A㊀㊀文章编号:1671 ̄5276(2022)01 ̄0105 ̄03KinematicsModelingandAnalysisofSeven-degree-freedomRoboticArmBasedonSpinorTheoryLIGuanqiꎬWUJiandeꎬLIRuiqin(InstituteofMechanicalEngineeringꎬNorthUniversityofChinaꎬTaiyuan030051ꎬChina)Abstract:Theroboticarmmimickingthe7DOFofthehumanarmwillhaveredundantdegreesoffreedom.Theforwardkinematicsolutionoftheseven-DOFroboticarmiscalculatedbasedonthespinortheoryꎬanditisverifiednumericallythatthematrixexponentialfirstblockexpansionismoreaccuratethanthedirectTaylorexpansion.TheNewton-Raphsonnumericaliterationmethodisusedtofindtheinversesolution.Theforwardandinversesolutionsareverifiedagainsteachotherbywritingamatlabprogram.Itisfoundthatthereisarangeofapplicabilityfortheinversesolutionꎬandtheuseofa0-pitchspiralaxisresultsinamoreconcisemodelingthanthetraditionalD-Hmethod.Keywords:spinortheoryꎻroboticarmꎻNewton-Raphsonnumericaliterativemethod0㊀引言机械臂是机器人技术领域中应用最为广泛的自动化装置ꎬ在工业制造㊁医疗康复等领域都有其应用实例ꎬ拟人化机械臂技术也逐渐成熟[1-4]ꎮ匹兹堡大学生物医学团队在BCI机械臂假肢里引入实时触觉反馈ꎬ使完成任务的时间减少一半ꎮ波士顿动力开发Stretch移动式7R机械臂已应用于仓库运输ꎮ20世纪80年代学者们开始讨论将旋转轴从数学中引入机器人研究的可行性ꎬ旋量理论[5-9]日益成熟ꎬ然而对其的应用介绍却鲜见报道ꎮ本文建模7自由度串联机械臂ꎬ着重对螺旋轴这一单位矢量在运动学的应用加以详细描述ꎮ6自由度逆解有通用公式ꎬ7R逆解包含32个实根ꎬ出现虚空间或自运动歧ꎮ7R逆解可结合Newton-Raphson数值迭代法求解[10]ꎮ使用Matlab进行正解中的化简与直接展开等运算过程的比较ꎬ并且同时使用开源Python对结果进行比较ꎮ使用基于空间(space简称s)坐标系的空间雅可比矩阵ꎬ找出数值迭代法求逆解的适用范围ꎬ基于正解的前提下ꎬ验证逆解的准确性ꎮ1㊀旋量理论螺旋(旋量 旋转向量)理论在正运动学中的应用ꎮ1.1㊀M矩阵(0位置㊁起始位置)M矩阵为当所有的关节转角都为0的时候ꎬ操作空间坐标系(body简称b)在s坐标系中的位置和姿态矩阵ꎮ1.2㊀旋转向量Sң与指数积公式将每一个转动关节视为0螺距的轴ꎬ假设除了最后一个关节转动θʎ外ꎬ其他关节都是固定的ꎮ在s坐标系里ꎬ关节7上螺旋轴的向量形式如式(1)所示ꎬ在高维度上ꎬRn是n维的欧几里得空间ꎬS7ңɪR6㊁ω7ңɪR3ꎬV7ңɪR3ꎮ螺旋轴的矩阵形式S7如式(2)所示ꎮS7ң=ωң7Vң7éëêêùûúú(1)S7=ω7Vң700éëêêùûúú(2)T07=eS7θ7M(3)501 博看网 . All Rights Reserved.特殊正交群是所有有效的3ˑ3旋转矩阵的集合群SO(3):包括R㊁ω㊁eωθꎮ特殊的欧几里得集合群或刚体运动群或R3中的同质变换矩阵Se(3)表示位姿:包括M㊁eSθ㊁T07ꎮeSθ=I+Sθ+S2θ22!+S3θ33!+ =eωθf(θ)ν01éëêêùûúú(4)式(4)中ꎬ可利用特性ω3=-ω来化简ꎬ且eωθ有Rodrigues公式:f(θ)=Iθ+(1-cosθ)ω+(θ-sinθ)ω2(5)eωθ=I+sinθω+(1-cosθ)ω2(6)依次解锁一个角度ꎬ往前代值ꎬ得到T07=eS1θ1 eS7θ7M(7)2㊀数值法逆解使用非线性寻根的Newton-Raphson方法ꎬ有寻根㊁不存在根时寻找近似解㊁存在多个解时寻找最优解的优势ꎮ给定一个初始值ꎬ然后代入迭代式求解直到出现误差范围内的解ꎮ设正向运动学函数为f(θd)ꎬ末端执行器的位置向量为νꎬ非线性寻根Newton-Raphson是找到目标函数的解ꎮ几何视角如图1所示ꎻ式(8)是解析视角ꎮ式(9)-式(11)是计算雅可比矩阵ꎮ图1㊀迭代法几何过程Δθ=J-1(θ0)[vd-f(θi)](8)JS(θ)=JS1(θ1)ңJS2(θ2)ң JSn(θn)ң[](9)JS1(θ1)ң=S1ң(10)JSi(θi)ң=[eS1θ1 eSi-1θi-1]∗Siң(11)式中∗为其伴随矩阵ꎮ对这种算法的进一步改进:1)末端执行器的位置描述f(θd)变更为正运动学计算出的矩阵T07ꎮ2)误差调整ꎮ用螺旋轴SMTң两分量的模代替末端执行器每次迭代的位置变化ꎮ3)引进伪逆矩阵J†避免求解奇异时无解的情况ꎮ在Matlab里编程为pinv(J)ꎮ当前基于b坐标系ꎬ变换为基于s坐标系:SMTbң=log(TMT07(θi))(12)SMTsң=(Tsb)∗SMTbң(13)式中:Tsb是坐标转移矩阵ꎻ∗为求其伴随矩阵ꎮ改进算法的流程图如图2所示ꎮ图2㊀改进的数值迭代法流程图3㊀验证正逆解3.1㊀正运动学方程㊀绘制三维模型图ꎬ并建立7R示意图(图3)进行验证ꎮ如图3(b)所示的坐标平面ꎬy方向定义为a并依次标号ꎬz方向定义为b也依次序标号ꎮ标示7个螺距为0的右手螺旋轴S1ң-S7ңꎮM=1000010240000100001éëêêêêùûúúúú(14)S7ң=[0㊀0㊀1㊀-a7㊀0㊀0]T=[0㊀0㊀1㊀-300㊀0㊀0]T(15)U a U 7Rb 7R .图3㊀三维模型图及7R示意图根据式(3)ꎬ使用矩阵分块后化简的运算过程ꎬ编程并计算结果ꎮT(θ)=c7-s70-2700s7s7c702700c7-30000100000éëêêêêùûúúúú(16)根据式(3)ꎬ发现将矩阵指数直接泰勒展开(采用了3种计算方法:Pada法㊁特征值法㊁6次的泰勒展开)得出的601结果虽然一致ꎬ但是它是虚数形式ꎬ增加了计算量ꎮ结果:TM_pada(θ)=c7-s70-2700s7s7c701350e-θ7i+1350eθ7i-30000100001éëêêêêùûúúúú(17)使用先进行矩阵分块然后利用特性化简的运算方法ꎬ该结果更准确㊁后期的运算量更小ꎮ代入式(7)ꎬ使用该方法分别在Matlab和Python中运算ꎬ得出的结果一致:T=c5-6c1234c7-s1234c7-s1234c7-c5-6c1234s7-s5-6c1234r1c1234s7+c5-6s1234c7c1234c7-c5-6s1234s7-s5-6s1234r2s5-6c7-s5-6s7c5-6r30001éëêêêêêùûúúúúúr1=150s1234-5-1350s12347-300s123-675s123457-6-675s123467-5-150s12345-1350s1234-7+600s1234+300s12+675s12345-6-7+675s12346-5-7-300s1+750s12345-6-750s12346-5r2=150c1234-5+1350c12347+300c123+675c123457-6+675c123467-5-150c1234-5+1350c1234-7-600c1234-300c12-675c12345-6-7-675c12346-5-7+300c1-750c12345-6+750c12346-5r3=2700s7c5s6-s5c6-1500s5s6-c5(1500c6-1500)-1200c5+12003.2㊀逆运动学方程1)第一次验证设初始位置为正解的0位置ꎬ转动角度为θlistꎮθlist=πꎻπ2ꎻπ3ꎻπ4ꎻπ5ꎻπ6ꎻπ7[]=[3.142ꎻ1.571ꎻ1.047ꎻ0.785ꎻ0.628ꎻ0.524ꎻ0.449]给逆解的初始值[3ꎬ1.5ꎬ1ꎬ0.6ꎬ0.5ꎬ0.4ꎬ0.3]ꎮ解得[3.161ꎬ1.563ꎬ1.059ꎬ0.762ꎬ0.628ꎬ0.523ꎬ0.449]ꎮ2)第二次验证因发现误差较大:1)改用角度制ꎻ2)迭代次数增加到1000次ꎻ3)精度调整为eω<0.0001ꎬev<0.0001ꎮ经过多组数据实验发现ꎬ相差3ʎ以内ꎬ位置完全重现ꎻ相差10ʎ以内ꎬ第一角度有0.1ʎ的偏差ꎻ相差15ʎ以上的逆运算ꎬ第一个角度有1ʎ以上的偏差ꎮ所以ꎬ该方法求逆解有完全重现的适用范围ꎬ需要把初始解猜测在真实解的附近15ʎ以内(表1)ꎮ表1㊀正逆解互相验证单位:(ʎ)㊀验证项目θ1θ2θ3θ4θ5θ6θ7Δθ初始角度10203040506070 猜测角度7172737475767计算的逆解102030405060700猜测角度6162636465666计算的逆解10.00120.00129.99939.9995060700.001猜测角度1112131415161计算的逆解10.0620.03229.94939.9595060700.060猜测角度0102030405060计算的逆解10.120.05329.91639.9325060700.100猜测角度-551525354555计算的逆解10.7420.40329.37339.4845060700.740续表1验证项目θ1θ2θ3θ4θ5θ6θ7Δθ猜测角度-641424344454计算的逆解11.0320.56629.12839.2765060701.030猜测角度-13-3717273747计算的逆解26.72824.46724.37534.4295060706.728猜测角度-15-5515253545计算的逆解24.46633.10318.25724.17350607014.466猜测角度-16-7414243444计算的逆解63.432-26.264-22.72585.55750607053.4324㊀结语本文重点介绍了如何使用形如螺旋楼梯的转向量计算正逆解ꎮ结合旋量与Newton-Raphson数值法求逆解ꎮ为涵盖转动机构可能出现的问题ꎬ采用目前串联机构中最复杂的7R机构并且使机构尽可能复杂ꎬ但是设计中没有涵盖helical螺旋和cylindrical圆筒等机构ꎮ在比较了不同的运算方法后ꎬ计算正解ꎬ得出先将矩阵分块㊁再利用特性化简的方法更好的结论ꎮ在计算逆解的过程中ꎬ发现当初始猜测的第一个角度超过真实解20ʎꎬNewton-Raphson数值法不会重现正解ꎮ验证结果表明旋量形式美观ꎬ建模快捷ꎬ编程明了ꎬ适用于机械臂中的串联机构㊁并联机构㊁转动关节㊁平动关节㊁螺旋关节等的运动学㊁动力学特性建模和分析ꎮ参考文献:[1]FLESHERSNꎬDOWNEYJEꎬWEISSJMꎬetal.Abrain-computerinterfacethatevokestactilesensationsimprovesroboticarmcontrol[J].Scienceꎬ2021ꎬ372(6544):831 ̄836.[2]KEVINMLꎬFRANKCP.Modernrobotics:mechanicsplanningandcontrol[M].Illinois:CambridgeUniversityPressꎬ2017.[3]刘世平ꎬ曹俊峰ꎬ孙涛ꎬ等.基于BP神经网络的冗余机械臂逆运动学分析[J].中国机械工程ꎬ2019ꎬ30(24):2974 ̄2977ꎬ2985.[4]赵京ꎬ王鑫ꎬ张自强ꎬ等.基于肘部自运动的主从异构7自由度机械臂运动映射及其几何逆解[J].机械工程学报ꎬ2020ꎬ56(15):181 ̄190.[5]YOUWSꎬLEEYHꎬOHHSꎬetal.Designofa3D-printableꎬrobustanthropomorphicrobothandincludingintermetacarpaljoints[J].IntelligentServiceRoboticsꎬ2019ꎬ12(1):1 ̄16.[6]BINDURAꎬNELOYAAꎬALAMSꎬetal.Sigma-3:Integrationandanalysisofa6DOFroboticarmconfigurationinarescuerobot[C]//20194thInternationalConferenceonRoboticsandAutomationEngineering(ICRAE).Singapore:IEEEꎬ2019:6 ̄11.[7]常健ꎬ王亚珍ꎬ李斌.基于力/位混合算法的7自由度机械臂精细操控方法[J].机器人ꎬ2016ꎬ38(5):531 ̄539.[8]张昌ꎬ武玉强.基于P-Rob六自由度机械臂运动学建模与仿真[J].包装工程ꎬ2020ꎬ41(11):166 ̄173.[9]WIEDMEYERWꎬALTOÉPꎬAUBERLEJꎬetal.Areal-time-capableclosed-formmulti-objectiveredundancyresolutionschemeforseven-DoFserialmanipulators[J].IEEERoboticsandAutomationLettersꎬ2021ꎬ6(2):431 ̄438.收稿日期:20210401701 博看网 . All Rights Reserved.。
基于7自由度模型的整车振动分析
基于7自由度模型的整车振动分析
7自由度模型是指将整车简化为具有七个独立振动自由度的系统,包括车体的纵向、横向和垂向振动、车体的横摇、纵摇和俯仰振动以及车轮的转动。
这个模型的基本假设是车体和车轮之间的振动是相互独立的。
整车振动分析的目的是评估整车在不同运动条件下的振动响应,找出潜在的振动问题并提供解决方案。
以下是整车振动分析的主要内容:
1.路面激振分析:通过设置不同的路面输入谱,模拟车辆在真实道路上行驶的振动情况。
通过分析车体和车轮的响应,可以评估整车的舒适性和稳定性。
2.悬架系统分析:悬架系统是减震和减振的关键部件,它可以减少车辆振动对乘坐舒适性和行驶稳定性的影响。
整车振动分析可以评估悬架系统的性能,并优化悬架参数以提高车辆的悬挂性能。
3.车体结构分析:车体结构的刚度和材料特性对整车的振动性能有着重要影响。
通过建立车体的有限元模型,可以评估不同结构参数和材料的振动响应,并进行优化设计。
4.动力系统分析:动力系统的工作也会对整车的振动产生影响。
通过分析发动机和传动系统的振动响应,可以评估动力系统的性能,并优化系统的设计以降低振动。
整车振动分析可以采用各种计算方法和仿真工具进行,包括有限元方法、多体动力学模拟、多物理场耦合仿真等。
这些分析结果可以为整车设计提供指导,并减少开发过程中出现的振动问题。
总之,基于7自由度模型的整车振动分析是一种有效的评估整车振动性能的方法,可以帮助优化整车的舒适性、稳定性和振动抑制能力。
汽车 七自由度动力学模型 作用
汽车七自由度动力学模型作用
汽车的七自由度动力学模型是一种用于描述汽车运动的数学模型,它考虑了汽车在空间中的三个平移自由度和三个旋转自由度,以及一个额外的自由度,通常用于描述车辆的垂直运动(车身的垂直加速度)。
以下是汽车七自由度动力学模型的主要作用:
仿真与建模:汽车七自由度动力学模型可以用于仿真和建模车辆在不同路况、操控条件下的运动行为。
通过模拟车辆的运动特性,可以评估车辆的稳定性、操控性和舒适性等方面的性能。
车辆控制系统设计:汽车七自由度动力学模型对于设计车辆的控制系统具有重要意义。
通过对车辆动力学特性的深入理解,可以优化车辆的悬架系统、制动系统、转向系统等,提高车辆的安全性和稳定性。
底盘调校与优化:基于七自由度动力学模型的仿真分析可以帮助工程师进行底盘调校和优化。
调校车辆的悬架几何、阻尼特性、操纵杆参数等,以实现更好的悬架调谐和操控性能。
碰撞安全评估:七自由度动力学模型可以用于评估车辆在碰撞事件中的动力学响应,包括车身的变形、乘员的受力情况等。
这有助于设计更安全的车身结构和安全气囊系统。
驾驶辅助系统设计:对车辆动力学行为的深入理解可以指导驾驶辅助系统的设计和优化,如防抱死制动系统(ABS)、电子稳定控制系统(ESC)等,提高车辆在各种路况下的稳定性和操控性。
总的来说,汽车七自由度动力学模型是汽车工程中的重要工具,对于设计、优化和评估车辆性能具有重要的作用,可以帮助提高车辆的安全性、操控性和舒适性。
1。
引入橡胶衬套的车辆七自由度振动模型分析
1.2 橡胶衬套粘弹性模型
橡胶衬套的粘弹性是通过动刚度和损耗角进行 频率相关性描述。目前,通常描述橡胶衬套粘弹性 的模型有:Kelvin-Voigt模型、三参数Maxwell模型和 Frequency-Dependent模型。三种粘弹性模型的具体内容 如表1所示。
表1 粘弹性模型图
模
Kelvin-Voigt
综上所述,国内外学者研究整车振动模型对车辆平
顺性的影响时,建立的整车振动模型很少考虑到橡胶衬 套力学特性模型,研究将橡胶衬套力学特性结合到整车 模型中更加匮乏,这对平顺性分析有一定的影响。车辆 是一个复杂的振动系统,在研究车辆平顺性时,整车振 动模型通常采用七自由度模型。因此,本文通过构建橡 胶衬套力学特性模型,建立引入橡胶衬套静态和动态力 学特性模型的整车振动模型,通过MATLAB建立其数 学模型并进行仿真,将经典七自由度振动数学模型与引 入橡胶衬套静态、动态力学特性模型的整车振动数学模 型对平顺性评价指标的影响进行对比分析,研究橡胶衬 套对于车辆平顺性的影响。
型
模型
三参数 Maxwell
模型
FrequencyDependent
模型
模 型 图
通过橡胶衬套的静态力学试验和动态力学试验,可 得到橡胶衬套的线刚度、扭转刚度、动刚度以及损耗角 等参数。
基于橡胶衬套的试验数据和粘弹性的三参数 Maxwell模型、Kelvin-Voigt模型和Frequency-Dependent 模型,分别对橡胶衬套的动刚度及损耗角进行拟合对 比,其结果如图1~图3所示。图1~图3中,虚线为试验 数据;实线为模型拟合曲线。
1)车辆行驶状态为匀速行驶,静平衡下的车辆位 移为线性位移;
2)模型忽略轮胎阻尼,用等效刚度代替; 3)车辆簧载与非簧载质量视为刚体,且弹簧、减 震器、轮胎位移可以通过刚度线性表示; 4)路面激励为垂向方向,车轮、减震器等只考虑 垂向的运动、刚度等; 5)车辆左右对称。 图中:q1,q2,q3,q4为分别为路面激励;Z1,Z2, Z3,Z4分别为轮心位移; m1,m2,m3,m4为分别为簧下质量;Kt1,Kt2, Kt3,Kt4 分别为轮胎刚度; Kd1,Kd2,Kd3,Kd4分别为悬架刚度;Cd1,Cd2,
(完整word)七自由度整车模型及参数
七自由度整车模型及参数七自由度线性整车模型如图1.1所示.图中各符号意义如下:s M 、θI 、φI ——悬挂质量、悬挂质量的侧倾转动惯量和俯仰转动惯量;1t m 、2t m 、3t m 、4t m ——非悬挂质量(分别为前左、前右、后左、后右,下同); 1s k 、2s k 、3s k 、4s k ——悬架刚度系数; 1t k 、2t k 、3t k 、4t k ——轮胎刚度; 1s c 、2s c 、3s c 、4s c -—阻尼器阻尼系数; 1u 、2u 、3u 、4u —-作用于悬架的控制力;1r x 、2r x 、3r x 、4r x —-地面扰动输入;1t x 、2t x 、3t x 、4t x —-非簧载质量位移;1s x 、2s x 、3s x 、4s x ——悬挂质量与悬架连接处的位移;c x 、θ、φ-—悬挂质量的垂直位移、侧倾角、俯仰角;xf l 、xr l -—悬挂质量质心至前后车轴的距离;ylf l 、ylr l --前后悬挂质量质心至左轮的距离.图1.1 七自由度整车模型令地面扰动输入向量T r r r r x x x x w ][4321=、车轮位置向量T t t t t t x x x x x ][4321=、悬挂质量运动向量T c C x X ][φθ=、悬架控制力向量T u u u u u ][4321=、悬挂质量与悬架的四个连接点处的位置向量T s s s s t x x x x x ][4321=、悬架动挠度向量T st st st st st x x x x x ][4321=(1st x 、2st x 、3st x 、4st x 分别表示前左、前右、后左、后右悬架动挠度),易知,t s st x x x -=。
根据悬架的特点和几何关系可以得出:C s HX x = (1)式中 ⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡----=xr ylrxr ylr xf ylfxf ylf l l l l l ll l H 1111 根据牛顿第二定律,系统总的运动方程可表示为:)()()()()(w x K u x X H C x HX K xM u H x X H C H x HX K H X M tttCStCstt T tC S T t C s T C b ----+-=+----= (2)式中),,(φθI I M diag M s b =,),,,(4321s s s s s k k k k diag K =,),,,(4321s s s s s c c c c diag C =,),,,(4321t t t t t m m m m diag M =,),,,(4321t t t t t k k k k diag K = 令[]T t C x X x=ˆ则上述运动微分方程可以表示为: u E w E x K x D xM 21ˆˆˆ+=++ (3) 式中⎥⎦⎤⎢⎣⎡=t bM M M ,⎥⎦⎤⎢⎣⎡+--=t s s s T s T K K H K K H H K H K ,⎥⎦⎤⎢⎣⎡--=s s s T s T C HC C H H C HD ,⎥⎦⎤⎢⎣⎡=t K E 01,⎥⎦⎤⎢⎣⎡-=I H E T 2 选取状态向量[]T x x x ˆˆ~= ,则式(3.5)的状态空间形式为: u B w B x A x 21~~~~~++= (4)式中 ⎥⎦⎤⎢⎣⎡--=--D M K M I A 110~,⎥⎦⎤⎢⎣⎡=-1110~E M B ,⎥⎦⎤⎢⎣⎡=-2120~E M B。
车辆系统动力学-复习提纲
1. 简要给出完整约束与非完整约束的概念2-23,24,25,1)、约束与约束方程一般的力学系统在运动时都会受到某些几何或运动学特性的限制,这些构成限制条件的具体物体称为约束,用数学方程所表示的约束关系称为约束方程。
2)、完整约束与非完整约束如果约束方程只是系统位形及时间的解析方程,则这种约束称为完整约束。
完整约束方程的一般形式为:式中,qi为描述系统位形的广义坐标(i=1,2,…,n);n为广义坐标个数;m为完整约束方程个数;t为时间。
如果约束方程是不可积分的微分方程,这种约束就称为非完整约束。
一阶非完整约束方程的一般形式为:式中,qi为描述系统位形的广义坐(i = 1, 2, …,n);为广义坐标对时间的一阶与数;n为广义坐标个数;m为系统中非完整约束方程个数;t为时间。
2. 解释滑动率的概念3-7,81.滑动率S车轮滑动率表示车轮相对于纯滚动(或纯滑动)状态的偏离程度,是影响轮胎产生纵向力的一个重要因素。
为了使其总为正值,可将驱动和被驱动两种情况分开考虑。
驱动工况时称为滑转率;被驱动(包括制动,常以下标b以示区别)时称为滑移率,二者统称为车轮的滑动率。
参照图3-2,若车轮的滚动半径为rd,轮心前进速度(等于车辆行驶速度)为uw,车轮角速度为ω,则车轮滑动率s定义如下:车轮的滑动率数值在0~1之间变化。
当车轮作纯滚动时,即uw=rd ω,此时s=0;当被驱动轮处于纯滑动状态时,s=1。
3. 轮胎模型中表达的输入量和输出量有哪些?3-22,23轮胎模型描述了轮胎六分力与车轮运动参数之间的数学关系,即轮胎在特定工作条件下的输入和输出之间的关系,如图3-7所示。
根据车辆动力学研究内容的不同,轮胎模型可分为:(1)轮胎纵滑模型主要用于预测车辆在驱动和制动工况时的纵向力。
(2)轮胎侧偏模型和侧倾模型主要用于预测轮胎的侧向力和回正力矩,评价转向工况下低频转角输入响应。
(3)轮胎垂向振动模型主要用于高频垂向振动的评价,并考虑轮胎的包容特性(包含刚性滤波和弹性滤波特性)。
利用七自由度车辆模型估计汽车状态参数_管欣
16 期
管
欣, 等: 利用七自由度车辆模型刚度和阻尼矩阵分别表示为
1
1. 1
七自由度车辆模型
模型概述 本文建立的七自由度车辆模型如图 1 所示。假
定车身为一刚性体, 四个车轮为质量点。 悬架类型 可简化为四个独立悬架,B1 — B4 四个点分别表示 四个独立悬架和车体的链接点。 七个自由度分别 侧倾角 φ 、 俯仰角 θ 和四个 是车身的垂向位移 Z c 、 Z u2 , Z u3 , Z u4 , 车轮的垂向位移 Z u1 , 用列向量 Z 表 示: Z = [ Zc φ θ Z u1 Z u2 Z u3 Z u4
A∈Rn ˑ n , B∈Rn ˑ m , C∈Rp ˑ n , D∈Rn ˑ q , p≥q。 式( 9) 中, C, D 满秩, 假设矩阵 B, 干扰 ξ 未知但有界, 即 d > 0。 ξ< d, 假设存在线形坐标变换 T0 , 使得动态系统可变 成式( 10 ) 。 1 ( t) = A11 x1 ( t) + A12 y( t) + B1 u( t) x
[ 1 —2 ]
。
本文针 对 整 车 主 动 悬 架 控 制 过 程 中 较 难 直 提出基于滑模观测 接测得的车辆关键 状 态 参 数 , 器理论 , 建 立 了 包 含 车 身 的 垂 向 振 动 加 速 度、 侧 倾角 、 俯仰角和四个 车 轮 的 垂 向 运 动 在 内 的 汽 车 七自由 度 动 力 学 模 型 和 状 态 估 计 模 型 。 文 中 对 该状态估计模型进行了仿真分析和场地实验验 采用该方法可以有效地 证分析 。 分析结果表明 , 估计汽车的侧倾角 速 率 、 俯仰角速率和垂向加速 度 。 该方法 还 为 汽 车 先 进 控 制 系 统 中 的 状 态 参 数估计提供了一种准确且低成本的实时软测量 技术 。
最新汽车2自由度和7自由度动力学建模仿真
1 路面模型的建立在分析主动悬架控制过程时,路面输入是一个不可忽略的重要因素,本文利用白噪声信号为路面输入激励,)(2)(2)(twUGtxftxggππ+-=•其中,f为下截止频率,Hz;G0为路面不平度系数,m3/cycle;U0为前进车速,m/sec;w为均值为零的随机输入单位白噪声。
上式表明,路面位移可以表示为一随机滤波白噪声信号。
这种表示方式来源于试验所测得的路面不平度功率谱密度(PSD)曲线的形状。
我们可以将路面输入以状态方程的形式加到模型中:⎪⎩⎪⎨⎧=+=•XCYWFXAXroadroadroadroadroad1,2,2,==-==roadroadroadgroadCUGBfAxXππ;D=0;考虑路面为普通路面,路面不平系数G=5e-6m3/cycle;车速U=20m/s;建模中,路面随机白噪声可以用随机数产生(Random Number)或者有限带宽白噪声(Band-Limited White Noise)来生成。
本文运用带宽白噪声生成,运用MATLAB/simulink建立仿真模型如下:图1 路面模型2 汽车2自由度系统建模图2 汽车2自由度系统模型根据图2所示,汽车2自由度系统模型,首先建立运动微分方程:()()()()()b b s b w s b w w w t w g s b w s b w m x K x x C x x m x K x x K x x C x x =----⎧⎪⎨=--+-+-⎪⎩整理得:⎪⎪⎩⎪⎪⎨⎧+--+-+-+-=-+-+-+-=gw t b w t s b w s b w s b w s w b b s b b s w b s b s bx m K x m K K x m K x m C x m C x x m K x m K x m C xb m C x式中:s C 为悬架阻尼,s K 为悬架刚度,t K 为轮胎刚度,b m 为车身质量,w m 为车轮质量,b b b x xx 、、分别为车身位移、速度、加速度,w w w x x x 、、分别为车轮位移、速度、加速度,g x 为路面输入。
汽车7自由度振动模型建模及在汽车平顺性能仿仿研究中的应用
汽车7自由度振动模型建模及在汽车平顺性能仿仿研究中的
应用
张庆才;雷雨成
【期刊名称】《上海汽车》
【年(卷),期】1999(000)003
【摘要】本文叙述了采用多刚体系统动力学建立汽车7自由度振动模型的方法。
并利用该模型对某汽车的平顺性进行了研究。
结果证明了该建模方法的正确性和数值计算的稳定性,同时也证明了设计结果的可行性。
【总页数】3页(P16-18)
【作者】张庆才;雷雨成
【作者单位】同济大学;同济大学
【正文语种】中文
【中图分类】U461.56
【相关文献】
1.坐姿人体四自由度动力学模型研究r——集中参数模型及其在汽车乘坐舒适性研究中的应用 [J], 白先旭;程伟;徐时旭;钱立军
2.汽车零件仿金涂装及应用研究 [J], 吴涛;王纳新;牛丽媛
3.数学建模在高职汽车类专业中的应用研究
——以汽车装配顺序优化为例 [J], 王雪
4.数学建模在高职汽车类专业中的应用研究--以汽车装配顺序优化为例 [J], 王雪
5.论仿拟在汽车广告语翻译中的应用 [J], 李卫丽;朱虹
因版权原因,仅展示原文概要,查看原文内容请购买。
七自由度车辆数学模型
图一.七自由度车辆动力学模型纵向力平衡方程:()()cos ()sin x y xfl xfr yfl yfr xrl xrr m V r V F F F F F F δδ•-⋅=+-+++ ①侧向力平衡方程:()()sin ()cos y x xfl xfr yfl yfr yrl yrr m V r V F F F F F F δδ•+⋅=+++++ ②绕Z 轴力矩平衡方程:12[()sin ()cos ][()cos ()sin ]2()()2w z xfl xfr yfl yfr xfr xfl yfl yfr w xrr xrl yrl yrr t I r F F F F a F F F F t F F F F b δδδδ•⋅=++++-+-+--+ ③四个车轮的力矩平衡方程:tw w xi bi di I R F T T i w •⋅=-⋅-+ ④上述方程中:δ为前轮转角;Vx ,Vy 分别为纵向、横向车速;β为质心侧偏角;γ为横摆角速度;Fxi 、Fyi 、Fzi 分别为轮胎纵向力、侧向力、垂向力;i=左前轮fl 、右前轮fr 、左后轮rl 、右后轮rr ,为分别对应的车轮;m 为整车质量;ms 为悬挂质量;a 、b 为前后轴到质心的距离;l=a+b 为前后轴距;tw1为前轴轮距;tw2为后轴轮距; Iz 为整车绕Z 轴的转动惯量;hg 为质心到地面的距离;各轮胎垂向载荷公式:_1_1_1_122222222gg z fl x y w gg z fr x y w gg z rl x y w gg z rrx y w h h b b F mg mV mV l l t l h h b b F mg mV mV l l t lh h aa F mg mV mV l l t lh h a a F mg mV mV l l t l••••••••=--⋅=-+⋅=+-⋅=++⋅(此处是Ax Ay )⑤各轮胎侧偏角公式:11222222arctan()arctan()arctan()arctan()www w y fl t x y fr t x y rl t x y rr t x V arV rV arV rV br V r V brV rαδαδαα+=--+=-+-=---=-+⑥各车轮轮心在车轮坐标系下的纵向速度:1_1_2_2_()cos ()sin 2()cos ()sin 222w t fl x y w t fr x y w t rl x w t rr x t V V r V ar tV V r V ar t V V rtV V rδδδδ=-++=+++=-=+ ⑦ 其中_t fl V 、_t fr V 、_t rl V 、_t rr V 为轮胎坐标系下的轮胎纵向速度。
车辆多自由度振动系统的建立
车辆多自由度振动系统的建立【摘要】为了评价汽车的平顺性,建立了汽车七自由度振动模型。
通过拉格朗日方程推导出了数学模型。
依托MATLAB软件,结合编程语言与SIMULINK 平台仿真分析了模型的频域响应和时域响应。
结果表明,此种多自由度振动模型符合实车情况。
【关键词】车辆;多自由度;振动模型1.引言随着生活水平的提高和生活节奏的加快,人们对汽车乘坐舒适性的要求越来越高[1]。
汽车平顺性是影响汽车乘坐舒适性的重要原因[2]。
所谓汽车平顺性是指汽车在一般行驶速度范围内行驶时,能保证成员不会因车身振动而引起不舒服和疲劳的感觉,以及保证所运货物完整无损的性能[3]。
车辆的振源主要有路面不平度,发动机,传动轴不平衡,轮胎以及侧向风的激励[4]。
一般汽车参数模型越多就越接近实际情况;然而参数过多就会导致模型复杂,求解困难[5]。
因此,合理选择车辆振动系统的自由度数至关重要。
2.汽车多自由度振动系统模型在汽车振动模型中,前轮与后轮可以近似为弹簧;前悬架与后悬架可以近似为弹簧与阻尼的并联;驾驶员座椅可以近似为弹簧与阻尼的并联;动力总成前悬置与后悬置同样可以近似为弹簧与阻尼的并联,如图1所示。
此模型具有七个自由度:前轮跳动,后轮跳动,车身跳动,车身俯仰,座椅跳动,动力总成跳动,动力总成俯仰。
模型中参数:mb为车身质量,Ib为车身绕惯性轴的转动惯量;kf为前悬架刚度,cf为前悬架阻尼,lf为前悬架到质心的距离,kr为后悬架刚度,cr为后悬架阻尼,lr为后悬架到质心的距离;mtf为前轮胎质量,ktf为前轮胎刚度,mtr为后轮胎质量,ktr为后轮胎刚度;mp为动力总成质量,Ip为动力总成绕惯性轴的转动惯量,lp为动力总成质心至汽车车身质心的水平距离,kp1为动力总成前悬置刚度,cp1为动力总成前悬置阻尼,lp1为动力总成前悬置到汽车质心的水平距离,kp2为动力总成后悬置刚度,cp2为动力总成后悬置阻尼,lp2为动力总成后悬置到汽车质心的水平距离;ms为座椅系统等效质量,ks为座椅系统等效刚度,cs为座椅系统等效阻尼,ls为座椅系统距离质心位置的水平距离;qf为来自于地面的前轮激励刚度,qr为来自于地面的后轮激励刚度。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
图一.七自由度车辆动力学模型
纵向力平衡方程:
()()cos ()sin x y xfl xfr yfl yfr xrl xrr m V r V F F F F F F δδ•
-⋅=+-+++ ①
侧向力平衡方程:
()()sin ()cos y x xfl xfr yfl yfr yrl yrr m V r V F F F F F F δδ•
+⋅=+++++ ②
绕Z 轴力矩平衡方程:
1
2
[()sin ()cos ][()cos ()sin ]2()
()2
w z xfl xfr yfl yfr xfr xfl yfl yfr w xrr xrl yrl yrr t I r F F F F a F F F F t F F F F b δδδδ•
⋅=++++-+-+--+ ③
四个车轮的力矩平衡方程:
tw w xi bi di I R F T T i w •
⋅=-⋅-+ ④
上述方程中:δ为前轮转角;Vx ,Vy 分别为纵向、横向车速;β为质心侧偏角;γ为横摆角速度;Fxi 、Fyi 、Fzi 分别为轮胎纵向力、侧向力、垂向力;i=左前轮fl 、右前轮fr 、左后轮rl 、右后轮rr ,为分别对应的车轮;m 为整车质量;ms 为悬挂质量;a 、b 为前后轴到质心的距离;l=a+b 为前后轴距;tw1为前轴轮距;tw2为后轴轮距; Iz 为整车绕Z 轴的转动惯量;hg 为质心到地面的距离;
各轮胎垂向载荷公式:
_1_1_1_122222222g
g z fl x y w g
g z fr x y w g
g z rl x y w g
g z rr
x y w h h b b F mg mV mV l l t l h h b b F mg mV mV l l t l
h h a
a F mg mV mV l l t l
h h a a F mg mV mV l l t l
•••••
•
••=--⋅
=-+⋅
=+-⋅
=++⋅
(此处是Ax Ay )
⑤
各轮胎侧偏角公式:
1
1
22222
2
arctan(
)arctan()
arctan(
)arctan(
)
w
w
w w y fl t x y fr t x y rl t x y rr t x V ar
V r
V ar
V r
V br V r V br
V r
αδαδαα+=--+=-+-=---=-+
⑥
各车轮轮心在车轮坐标系下的纵向速度:
1
_1_2
_2_()cos ()sin 2()cos ()sin 2
22w t fl x y w t fr x y w t rl x w t rr x t V V r V ar t
V V r V ar t V V r
t
V V r
δδδδ
=-
++=+++=-=+ ⑦
其中_t fl V 、_t fr V 、_t rl V 、_t rr V 为轮胎坐标系下的轮胎纵向速度。
各车轮滑移率的计算:
________fl
t fl
fl t fl
fr t fr
fr t fr rl t rl
rl t rl rr t rr
rr t rr
w R V V w R V V w R V V w R V V λλλλ-=-=-=-=
⑧。