低地板大型公共汽车翻滚碰撞分析印度博士巴巴萨海布安贝德卡M.docx
- 文档编号:23717555
- 上传时间:2023-05-20
- 格式:DOCX
- 页数:20
- 大小:2.78MB
低地板大型公共汽车翻滚碰撞分析印度博士巴巴萨海布安贝德卡M.docx
《低地板大型公共汽车翻滚碰撞分析印度博士巴巴萨海布安贝德卡M.docx》由会员分享,可在线阅读,更多相关《低地板大型公共汽车翻滚碰撞分析印度博士巴巴萨海布安贝德卡M.docx(20页珍藏版)》请在冰豆网上搜索。
低地板大型公共汽车翻滚碰撞分析印度博士巴巴萨海布安贝德卡M
低地板大型公共汽车翻滚碰撞分析
MarathwadaUniversity
2.2有限元模型
2.2.1划分网格
通过对车辆模型的分析计算获取车辆在碰撞过程中各组成部分和子系统变形及相互作用的信息。
对碰撞行为描述的精确性取决于CAD参数的设置和网格划分质量的好坏。
CAD几何模型的形状和尺寸应与实际车辆完全一致。
有限元模型网格划分密度要合理,使计算时间较为合理。
图2.1有限元建模流程图
图2.1所示为客车有限元建模方法。
使用HyperMesh划分网格,它是一款功能强大的有限元前后处理软件,能建立车辆有限元模型、显示结果及对进行数据分析。
首先利用Pro-E或CATIA建立CAD模型,并转化为IGES格式。
客车CAD模型由厂商提供,将CAD模型导入HyperMesh。
第一步从模型中提取中面,如图2.2所示。
第二步使用geomcleanup和defeature选项对中间层进行几何清理,为划分网格做好准备。
图2.2抽取中面
2.2.2网格质量评价标准
HyperMesh中有默认的网格质量评价标准,如下:
最小边长——一个网格单元最小的边长
最大边长——一个网格单元最大的边长
纵横比——一个网格单元最大边长与最小边长的比值
扭曲度——一个网格单元偏离平面的角度
最小(最大)四边形内角——四边形单元最小(最大)内角
最小(最大)三角形内角——三角形单元最小(最大)内角
三角形单元百分率——三角形单元占所有网格单元的比例
表2.1列出了网格质量评价标准,该标准适用于网格划分的全过程。
最小单元的尺寸不应小于5mm,以缩短网格划分的时间。
表2.1网格质量评价标准
2.2.3非结构部件
一些非结构部件如发动机、油箱、蓄电池隔间、车顶空调等与车身非常近,需要把它们的质量均分到各个节点上去。
由于这些部件在汽车翻滚中变形很小,所以定义它们为刚性体,只影响客车模型的质量和重心(CoG)。
车身内部的附件也要进行建模,以考察它们在客车翻滚试验中与车身有限元模型的相互作用。
为了使客车有限元模型的重心(CoG)与真实客车的重心(CoG)重合,还需要把一些质量分配到各个节点上去。
对发动机模型划分网格采用长方体单元,对轮胎模型划分网格采用封闭球体单元,这样便于定义内部压力。
图2.3客车非结构部件
2.3定义材料参数
车身结构使用的材料主要是钢和铝,而汽车碰撞是一个高速的动力学过程,在碰撞过程中客车车身各部构件以不同的应变率产生变形。
由于客车模型要在0~60Km/h的侧碰速度下产生不同的变形,因此需要考虑应变率对不同构件力学属性的影响。
为了获取各部构件材料参数和应力-应变曲线,分别采用三种装置对客车进行了碰撞试验:
(1)机械装置或液压伺服装置——准静态条件、应变率为小于0.1/s
(2)液压伺服装置——应变率为0.1~500/s
(3)分离式霍普金森杆装置——应变率为500~1000/s,或更大
2.4连接形式
在HyperMesh里fe-joint选项可以对连接单元进行创造、检查和优化。
由于连接单元默认为建立刚性体之间的连接,因此连接部位要定义为刚性件。
CONSTRAINED-EXTRA-NODES-SET选项用来创建连接单元和刚性件之间的连接。
每个连接单元都存有相关的参数,如连接形式、方位。
整个客车有限元模型共创建了32个连接,如表2.2所示。
这些连接用来定义客车的悬架系。
表2.2客车各构件间的连接形式
2.5悬架系
该客车采用空气悬架,它通过高度控制阀来调节空气弹簧内的体积,使车体保持合适的高度。
图2.4所示为运动学连接在客车前悬架中的应用,其中前轴的四个空气弹簧建模为具有非线性阻尼减震功能的平移连接,两个液压减震器建模时也采用这种形式;控制臂采用球连接;车轮采用卷曲连接。
空气弹簧和减震器的动力学曲线如图2.5所示:
图2.4车轴运动学连接模型
图2.5空气弹簧和减震器的动力学曲线
客车后轴的四个空气弹簧采用具有非线性阻尼减震功能的平移连接;两个液压减震器也采用这种连接形式;四个控制臂采用了八个球连接;车轮采用两个卷曲连接,如图2.6所示。
其空气弹簧和减震器的动力学曲线与图2.5相同。
客车有限元允许改变车体的高度和弹簧的压力。
图2.6客车后轴运动学连接模型
2.6各部件固有振动特性分析
整车模型由几大部件组装而成。
可变形部件之间采用点焊连接,Constrainedrigidbody选项用来定义刚性件与可变形件之间的连接,Rigidbodymerge选项用来定义刚性件之间的连接。
整车模型各组件如图2.7所示。
为了使整车装配合理,需要对各组件进行固有振动特性分析,使用的指令是CONTROL-IMPLICIT-GENERAL和CONTROL-IMPLICIT-EIGENVALUE。
分析时给各组件输入其固有频率,通过分析可以检查各种连接形式是否合理。
进行固有振动特性分析时计算时间短,但占用内存大。
图2.7客车模型组件
2.7整车装配
客车整车有限元模型如图2.8/2.9所示。
图示模型为客车结构模型,去掉了车身内部附件。
若添加车身内部附件,整车模型则有282025个单元、20306个点焊和26种材料,如表2.2所示。
图2.8客车有限元模型装配图
图2.9客车有限元模型
表2.2客车有限元模型参数
2.8加速度计
为了得到各点的加速度数据,以便于与实车碰撞试验进行对比,以及不同碰撞速度下各点加速度的对比,需要安装加速度计。
在客车模型中加速度计的安装位置如图2.10所示。
加速度计为刚性材料(MAT-RIGID)三维单元模型,受客车部件的约束。
一些节点的加速度受多因素的影响,没有参考意义,会带来误导。
在客车地板、侧梁、座椅、保险杠和顶部等重要部位都装有加速度计。
每个加速度计都有协调性的局部系统。
NODOUT文件包含了所有加速度计的位置参数和历史记录数据。
图2.10加速度计安装位置
DATABASE-SECFORC和DATABASE-CROSSSECTIONPLANE-ID指令用来测量剪切力,并以ASCII格式写入SECFORC文件。
2.9创建客车有限元模型乘员生存空间
生存空间在客车翻滚碰撞中保存完好能使乘员的生存率大大提高。
生存空间的定义以ECER66法规为依据。
生存空间的描述如图2.11所示,生存空间应与客车地板(客车翻滚过程中不变形)相连接,确保在客车翻滚过程中其位置不变,通过CONSTRAINED-EXTRA-NODES-SET指令可以完成此操作。
图2.11乘员生存空间
2.10客车模型顶部冲击试验
使用CATIA对车顶冲击架进行三维建模,然后导入HyperMesh划分网格,车顶冲击架如图2.12所示。
其网格尺寸与客车顶部相同,执行MAT-RIGID指令可以完成此操作。
去掉空气调节装置后客车模型的质量为9.712吨,设定车顶冲击架的密度参数使它的质量是客车质量的1.5倍,即为14.568吨。
图2.12车顶冲击架
进行客车顶部冲击试验时,客车地板置于一个刚性底板上,而不是由轮胎支持车体,这样确保了所有载荷都由车身结构承担,不会分配到悬架系上。
底板是固定不动的,即在任何方向上没有自由度。
冲击架的放置如图2.13所示,所有横梁都应在车顶上。
图2.13车顶静压试验
冲击架在自身重力的作用下冲击客车顶部。
冲击后客车各部件之间的接触由AUTOMATIC-SINGLE-SURFACE指令生成。
冲击架与客车顶部的接触由AUTOMATIC-SURFACE-TOSURFACE指令定义,该指令也用来定义客车地板与刚性底板之间的接触。
客车顶部冲击试验结果由生产商提供,如下:
客车顶部达到的载荷为
N
内部构件最大静态偏移量为111.2mm
如图2.14所示,客车翻滚碰撞仿真过程中内部构件的最大动态偏移量达到143mm,弹性变形恢复后,静态的塑性变形量为95~115mm,与真实试验结果相符。
支持底板承受的载荷为:
底板承受的载荷=客车质量+冲击架质量
=9.712+14.568=24.28吨=24280Kg
因此,支撑底板承受的压力=
=
N。
图2.14车顶静压试验结果
图2.15所示为最大偏移量区域受力变化曲线:
图2.15最大偏移量区域受力变化曲线
冲击区域受力变化曲线用来计算应变能,即曲线下方的面积就是碰撞中客车顶部吸收的能量。
如图2.16所示:
图2.16静态冲击试验冲击区域能量变化曲线
3结构分析
3.1基于ADAMSView进行ECER66翻滚碰撞试验
按照ECER66法规要求,翻滚试验台的几何形状如图3.1所示。
试验台应有足够的刚度,并至少可以转动1度的角度。
试验台车轮支撑平面与翻滚沟平面的高差为800
20mm。
翻滚转轴中心线与车轮支撑平面和翻滚沟垂直壁的最大距离为100mm。
车轮侧滑挡块用来防止试验台翻转时车轮发生侧滑,它的位置应尽量靠近翻滚转轴。
试验台的设计应避免车辆发生纵向移动。
翻滚沟碰撞区域应水平、干燥、光滑和足够坚固。
图3.1翻滚试验台
图3.2基于ADAMSView进行ECER66翻滚试验
如图3.2所示,创建一个与客车尺寸相同的重块。
重块的质量为客车去掉空气调节装置后的质量,即为9.76吨。
利用EasyCrashDyna计算客车模型的重心(CoG),由于客车的发动机等大质量部件置于后面,因此CoG靠近后轴。
该重块的重心应与客车模型的重心一致。
按照ECER66法规要求,试验台与碰撞块的高差应为800mm。
设置车轮侧滑挡块,防止试验台翻转过程中车轮发生侧滑。
试验台与地面之间采用卷曲连接。
法规要求试验台的最大翻转角速度为
。
碰撞块与地面之间采用锁销连接,即碰撞块不可动。
试验台车轮支撑平面和碰撞块平面的摩擦系数分别为0.2和0.7。
3.5考虑乘员重量基于LS-DYNA进行客车翻滚碰撞仿真
客车翻滚碰撞仿真在LS-DYNA中的实现如图3.7所示:
图3.7基于LS-DYNA进行客车翻滚碰撞仿真
表3.1无乘员客车模型碰撞前的线速度和角速度
碰撞面为零自由度的刚性面(MAT-RIGID)。
客车模型的线速度和角速度与ADAM仿真相同,如表3.1所示,通过INITIAL-VELOCITY-GENERATION指令进行分配。
客车模型和碰撞面之间的接触由AUTOMATIC-SURFACE-TO-SURFACE指令定义。
碰撞面的静摩擦系数和动摩擦系数都为0.7。
客车各构件之间的接触由CONTACT-AUTOMATIC-SINGLE-SURFACE-ID指令定义。
客车在翻滚碰撞过程中的变形如图3.8所示,可以看出没有发生大变形,不对乘员安全和生存空间的完好性构成威胁。
图3.9所示为仿真前后乘员生存空间的形态,可以看出没有遭受破坏,翻滚后依然完好。
Mises等效应力的分布如图3.10所示。
可以看出大部分碰撞能量被A柱、后纵梁和地板横梁所吸收,A柱和后纵梁在外力的作用下逐渐发生弯曲变形。
图3.8无乘员客车变形过程
图3.9无乘员客车生存空间
3.10无乘员客车Mises等效应力分布
3.11无乘员客车顶部冲击区域应力-应变曲线
客车顶部冲击区域应力-应变曲线如图3.11所示。
当变形量为230mm时应力达到峰值
N。
利用LS-POST将该曲线转化为应变能-应变曲线,如图3.12所示。
当变形量达到250mm时,冲击能量为
Nmm。
可以看出动态翻滚试验下冲击能量和最大变形量都要大于静态客车顶部冲击试验。
但其应变能-应变曲线的斜率要小于静态试验,这是因为静态试验中客车顶部吸收的冲击能量要多于动态翻滚试验。
3.12无乘员客车翻滚碰撞顶部冲击区域应变能-应变曲线
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 地板 大型 公共汽车 翻滚 碰撞 分析 印度 博士 巴巴 萨海布安贝德卡
![提示](https://static.bdocx.com/images/bang_tan.gif)