汽车前保险杠结构及安全性能分析学士学位论文.doc

汽汽车车前前保保险险杠杠结结构构及及安安全全性性能能分分析析 A STUDY ON THE STUCTURE AND THE SAFETY PERFORMANCE OF THE VEHICLE FRONT BUMPER 专专 业 机械设计制造及其自动化 汽车工程 业 机械设计制造及其自动化 汽车工程 摘 要 前保险杠系统作为安全防护装置是现代汽车结构的重要组成部分 它要在 低速碰撞 时 硬 些 防止汽车发生大的变形以致冲击侵入驾驶舱伤害乘 员的安全 同时 当汽车与行人相碰撞时要 软 些 以避免或减少汽车对行 人产生伤害 可见 它的结构对于汽车的安全性能起着至关重要的作用 它能 有效地减轻汽车损坏以及人员伤亡的情况 因此 随着人们安全意识的不断提 高 前保险杠系统成为了国内外汽车被动安全领域的重点研究对象 过去对于乘员的保护已经做了很多研究 所以本论文主要侧重于前保险杠 结构对行人的保护研究 根据我国现行的行人安全保护法规 本论文采用计算 机模拟仿真技术 针对行人下肢的防护 对汽车前保险杠的结构及其安全性能 展开研究 论文把人体下肢损伤机理作为理论依据 利用 CATIA 软件建立两款 汽车前保险杠系统的简化几何模型 应用 HYPERMESH 软件划分网格 引进了本 校损伤生物力学与车辆安全工程中心的人体下肢仿真模型 利用 PAM CRASH 软 件模拟行人下肢 汽车前保险杠碰撞试验并进行相应的求解计算 根据行人保护 法规规定行人下肢损伤的三项评价指标 胫骨加速度 膝关节弯曲角度和膝关 节剪切位移 对两款汽车前保险杠模型的碰撞试验结果进行分析比较 分析结 果表明 车型 1 的前保险杠系统能够有效保护行人的安全 然而车型 2 的前保 险杠系统不能为行人提供有效的防护 本文对其提出了改进的措施 本研究是对我国前保险杠安全性能课题的有益探索 致力于通过整合建模 到模拟试验的整个过程 总结出有效的 有利于行人保护的前保险杠安全性设 计的计算机仿真试验方法 关键词 碰撞安全 行人下肢 前保险杠系统 冲击损伤生物力学 ABSTRACT As a safety devise Front Bumper System is an important part of a modern automobile When a low speed collision occurs the Front Bumper System should be stiff enough to avoid large deformation of the car and protect the passengers from the impact of the collision Meanwhile when a pedestrian is hit by the front part of a car the Front Bumper System should be soft enough in case that it will do harm to the pedestrian Obviously the structure of the Front Bumper System plays a very important role in the safety performance of an automobile Therefore with people s increasing awareness of safety the Front Bumper System has already become a key problem in passive vehicle safety area at home and abroad Large numbers of studies on the protection for passengers have been done previously for which this paper mainly focuses on the one for pedestrians According to the lately pedestrian laws and regulations in China this paper mainly studies on the structure and the safety performance of the automobile by computer simulation for the protection of the lower limbs of the pedestrians This paper takes the injury mechanism of lower limb as theory builds the geometric model by CATIA meshes the model by HYPERMESH introduces the simulation model of lower limbs from the Center for Impact Injury Biomechanics and Automobile Safety in TUST and conducted the simulation test of the lower limb Front Bumper System collision by PAM CRASH On the basis of the evaluation criterions peak acceleration of the tibia bending angle of the knee lateral shear displacement at the knee ruled by the pedestrian laws and regulations comparing the analysis results of the tests it indicates that the Front Bumper System of VEHICLE 1 is able to protect the pedestrians effectively but the one of VEHICLE 2 couldn t provide a valid safeguard for the amelioration of which several methods is presented This paper is helpful to the development of pedestrian protection Through the whole process from model building to simulation test an efficient method of the design with computer simulation in passive vehicle safety is summarised for the benefit of pedestrian protection Key words Impact Pedestrian Leg Front Bumper System Impact Injury Biomechanics 目 录 1 前 言 1 1 1 本课题研究背景和意义 1 1 2 本文主要研究内容及方法 4 1 2 1 研究内容 4 1 2 2 研究方法 4 2 汽车前保险杠及行人损伤安全概述 7 2 1 汽车前保险杠概述 7 2 1 1 结构 7 2 1 2 安全性能 8 2 2 行人下肢损伤生物力学 9 2 2 1 冲击损伤生物力学 9 2 2 2 损伤机理 9 2 3 行人安全法规 9 2 4 本章小结 11 3 汽车前保险杠有限元模型的构建 12 3 1 汽车前保险杠几何模型的构建 12 3 1 1 几何尺寸的确定 12 3 1 2 几何模型的构建 13 3 1 2 1 CATIA 软件简介 13 3 1 2 2 V1 前保险杠系统几何模型的构建 13 3 1 2 3 V2 前保险杠系统几何模型的构建 15 3 2 汽车前保险杠有限元模型的构建 16 3 2 1 HYPERMESH 软件简介 16 3 2 2 有限元网格划分的注意事项 17 3 2 2 1 网格划分 17 3 2 2 2 单元特性 18 3 2 2 3 网格质量的控制 18 3 3 3 有限元网格的划分 19 3 3 3 1 V1 前保险杠系统模型有限元网格的划分 19 3 3 2 2 V2 前保险杠系统模型有限元网格的划分 21 3 4 本章小结 23 4 下肢 前保险杠计算机模拟碰撞试验 24 II 4 1 PAM CRASH 软件简介 24 4 2 仿真试验的前处理 25 4 2 1 基本假设 25 4 2 2 参数量纲 25 4 2 3 材料属性 25 4 2 4 厚度参数 26 4 2 5 边界条件 26 4 2 6 接触定义 27 4 2 7 求解时间 27 4 2 8 碰撞模型 27 4 3 碰撞模拟试验求解 28 4 4 仿真结果的后处理 29 4 5 仿真结果的分析与讨论 32 4 5 1 独立分析 32 4 5 2 比较分析 33 4 6 本章小结 33 5 结论与展望 34 5 1 结论 34 5 2 创新点 34 5 3 展望 35 参考文献 36 致 谢 37 天津科技大学 2013 届本科生毕业论文 1 1 前 言 1 1 本课题研究背景和意义 近现代以来 随着全球经济和科技水平不断提高 汽车工业迅猛发展 在 我国 汽车行业的自主品牌如雨后春笋般崛起 汽车市场不断发展壮大 机动 车保有量逐年猛增 汽车总销量已经超越美国 成为了全世界最大的汽车产销 市场 但随之而来的是各种交通事故的发生 我国道路交通情况复杂 是世界 上典型的以平面混合交通为主的国家 道路交通事故发生率高居世界前列 因 此 在注重汽车驾驶感受和动力性能体现的同时 越来越多的人们开始重视汽 车的安全性问题 汽车前保险杠作为汽车安全的第一道防护线 成为了国内外 汽车安全性研究的重要课题之一 1 早期的前保险杠系统的主要设计目的是防止汽车发动机舱及驾驶室发生过 大的变形从而保护车内乘员的人身安全 对于车内乘员安全的保护性能设计研 究工作已经开展了几十年 然而对于车外行人的安全保护近几年来才逐渐受到 人们的重视 随之汽车安全性能研究工作的不断深入与展开 诸如安全带 安 全气囊等形形色色的汽车被动安全装置都在汽车的安全性结构中得到广范的使 用 当汽车发生碰撞时这些装置能够较好的保护车内乘员的安全 但是道路上 的行人及其他车外人员因为缺乏保护而成为了最容易收到伤害的高危人群 2 某调研数据指出 每年的道路交通死亡事故当中 大约有 1 4 的死亡人口为行 人 而在城市中此比例甚至接近 50 2007 年的该比例竟高达 54 4 在这些 事故中 头部和下肢的损伤又占损失总数的 70 以上 行人头部的损伤是导致 行人死亡的主要原因 而行人下肢的创伤往往是交通事故中导致人员伤残的主 要因素 3 从图 1 1 可见 行人与汽车发生意外碰撞时 行人下肢的损伤机率 最大 而下肢的伤害尤其是膝关节伤害往往会给受害者将来的生活和工作带来 极大的不便 有的甚至会导致行人的终身残疾 前保险杠系统作为汽车车身最 前端的部分 碰撞时直接与行人的下肢发生接触 如图 1 2 所示 人 车碰撞事 故中 发生碰撞的部位主要为前保险杠 因此 前保险杠系统的设计是一辆汽 车能否满足行人安全保护性能要求的关键 4 天津科技大学 2013 届本科生毕业论文 2 图 1 1 行人与汽车前部碰撞的身体部位损伤分布百分比 图 1 2 人 车事故中车辆碰撞部位分布 再者 随着行人安全保护法规的实施 我国的汽车企业对于行人的碰撞保 护产品的设计和开发日益重视 我国各大汽车生产商以及汽车研发机构在设计 汽车前保险杠系统是开始引进行人碰撞保护试验系统 将行人下肢保护产品的 设计融入了开发的流程 2009 年 3 月 中国汽研中心从法国 BIA 公司引进了行 人碰撞保护试验系统 此举为我国相关法规和标准的出台奠定了硬件和技术基 础 随着行人碰撞保护试验系统的引进和相关行人安全保护法规的出台与实施 研究具有行人安全保护性能的汽车提上日程 4 总而言之 行人 汽车碰撞事故 中行人下肢防护技术以及前保险杠安全性结构设计的改进设计研究 对于降低 交通事故发生率 减免事故中的伤亡 减少财产损失等 在当前中国具有着十 分重要的科研意义 天津科技大学 2013 届本科生毕业论文 3 括加速度 应力和变形等运动及损伤情况 目前 对于汽车被动安全性研 究领域的研究 主要应用的软件有 CAL3D MADYMO LS DYNA3D 和 PAM CRASH 等 5 这些软件中 根据建模理论及其分析功能的不同 可以分为采显式有限元 理论建模和采用多刚体系统动力学理论建模两类 其所构建模型的主要特点如 表 1 1 1 所示 表 1 1 两种模型的主要特点 模型类别有限元模型多刚体模型 仿真软件PAM CRASH LS DYNA3D 等CAL3D MADYMO 等 建模理论显示有限元理论多刚体系统动力学理论 分析内容任意单元的位移 速度 加速度 力等动力学响应参数及其应力 应 变 能量分布 刚体单元的位移 速度 加速度 力等动力学响应参数 优点精确地获得模型上不同位置的仿真 结果情况 建模较便捷 计算时间较短 工 作量较小 缺点建模过程较复杂 计算时间较长 工作量较大 描述模型的运动方面 薄弱 较难掌握模型的具体情况的 变 形部件描述的方面薄弱 综合比较上诉软件的优缺点 PAM CRASH 软件有着极大的优势和强大的功 能 本文采用 PAM CRASH 软件进行计算机模拟试验 本文的研究目的是通过计算机模拟仿真 对某两款汽车前保险杠的结构及 安全性能进行分析与比较 得到更有利于保护行人安全的前保险杠结构 并总 结试验方法 主要从以下三个方面展开研究 1 查阅相关文献 对现有文献进行研读 总结国内外关于此课题的研究 现状 了解损伤生物力学的应用以及相关安全法规和标准的规定 2 利用有限元方法进行汽车碰撞数值仿真能够有效的降低运算成本 大 大的缩短设计周期 并且可以实现无样车的模拟碰撞试验 因此 近年来世界 各大汽车生产商正在重点应用有限元法进行设计与分析 本论文的研究试验流 程 依据所应用软件及工作性质先后分为以下六个阶段 图 1 1 依据原型绘制二维图 标注尺寸 据二维图构建几何模型 对几何模型划分有限元网格 天津科技大学 2013 届本科生毕业论文 4 仿真模型前处理 仿真试验求解 求解结果后处理 数据可视化 1 3 基于以上的研究 分析讨论仿真试验结果 最终得出结论 图 1 1 仿真试验流程图 几何模型 CATIA 网格划分 HYPERMESH 前处理 VISUAL HIV 求解 PAM CRASH 后处理 VISUAL VIEWER 二维图 CAD 天津科技大学 2013 届本科生毕业论文 5 图 2 3 行人下肢与汽车前部碰撞是的主要损伤形式 2 3 行人安全法规 目前 世界各组织对这三项指标的极限值有着不同的要求和规定 其中 我国于 2009 年 10 月 30 日发布我国第一部行人安全保护法规 GB T24550 2009 汽车对行人的碰撞保护 于 2010 年 7 月 1 日修改为采用全球技术法规 GTR9 关于机动车碰撞时对行人及弱势道路使用者加强保护和减轻严重伤害的 认证统一规定 2008 年英文版 并正式全国实施 成为我国推荐性的国家标 准 1 11 该行人安全保护法规规定 车辆以 40km h 的速度与行人下肢发生侧面 碰撞时 要求 1 胫骨加速度 at 170g 2 膝关节弯曲角度 k 19 3 膝关节剪切位移 Skt 6mm 表 2 1 4 1 4 图 2 2 下肢碰撞伤害指标 股骨骨折 股骨骨节骨折 韧带断裂 胫骨骨节骨折 胫骨骨折 腓骨骨折 胫骨加速度 膝关节剪切位移 膝关节弯曲 角度 天津科技大学 2013 届本科生毕业论文 6 3 汽车前保险杠有限元模型的构建 3 1 汽车前保险杠几何模型的构建 整车有限元模型的建立及计算需要花费的大量的精力和时间 且需要有充 足的资源支持 因此 在现有的资源和条件下 以两款汽车为建模的原型 并 定义两款汽车分别为 V1 和 V2 如图 3 1 所示 本文仅分别建立 V1 和 V2 的蒙 皮 泡沫及横梁的有限元模型 通过设置相应的材料属性及边界条件等来模拟 碰撞试验 a V1 建模的原型 b V2 建模的原型 图 3 1 两款汽车建模的原型 3 1 1 几何尺寸的确定 本节参考了两款汽车原型的前保险杠部分测量数据 加以修正及简化 并 经过计算和经验补全未知的尺寸数据 绘制出两款车型前保险杠系统各主要部 件的 CAD 二维图 由于汽车前保险杠系统的轮廓相当复杂 圆弧面和圆角众多 而复杂的曲 面往往会增加几何模型建模的难度 甚至降低网格质量 最终影响到试验分析 结果的准确性 因此 为了简化几何建模和网格划分的操作 提高有限元网格 的质量 节省计算时间 优化试验过程 本节对其原型的外形轮廓作以下处理 1 忽略原型中过小的斜度 2 忽略用于尺寸较小且不用于装配的圆孔 3 去除原型中的圆角和斜度较小的倒角 4 尺寸较小的曲线 尽可能采用圆弧代替 如 V1 的车标处 5 对于复杂多曲率曲面拼接的曲面 尽量简化为尽可能少的曲率平滑过 天津科技大学 2013 届本科生毕业论文 7 度的曲面 6 部分弧度较小的曲面用平面来代替 7 对于薄壁件 如 V1 的蒙皮和横梁 在忠于原型的前提下 尽可能使 用均匀的厚度尺寸 8 对于厚度较大的实体 如 V1 和 V2 的泡沫 同样的 在忠于原型的前 提下 尽可能使用相同的厚度尺寸 9 所有部件的模型均视为左右对称 10 将除厚度尺寸外的其他尺寸值圆整为整数 其他细微处的处理将体现于 CAD 工程图中 本文不再一一列举 通过这些 修正和简化 最终分别绘制得到 V1 和 V1 前保险杠系统各主要部件的 CAD 二维 图 3 1 2 几何模型的构建 3 1 2 1 CATIA 软件简介 CATIA 系列软件可以帮助客户实现产品的外型和机构设计 通过环境模拟 对产品性能进行分析和优化等 CATIA 软件已经成为汽车 航空航天 船舶 建筑 电力与电子 消费品和通用机械等产品的设计生产的首要 3D 模拟应用软 件 CATIA 软件能为各类车辆的设计和制造提供了端对端 end to end 的解 决方案 尤其在造型 车身结构及引擎设计等方面具有相当出色的表现 CATIA 的可伸缩性及其并行工程能力能够显著缩短产品的上市时间 凭着人性化和智 能化的功能优势 CATIA 成为了欧洲 北美乃至亚洲各大汽车制造商所用的核 心系统 3 1 2 2 V1 前保险杠系统几何模型的构建 本节应用 CATIA 软件分别构建 V1 前保险杠系统三大组成部分 蒙皮 泡沫 和横梁 的几何模型 V1 蒙皮为左右对称壳体 因此只需创建其左半边壳体的几何模型 其后使 用 镜像 命令即可获得整个 V1 蒙皮的几何模型 其左半边壳体靠近对称面部 分可以大致看成一各半月牙的平面沿曲率较小的弧线 扫略 而成的壳体 其 远离对称面部分可以大致看成一各半月牙的平面沿曲率较大的弧线 扫略 而 成的壳体 因此 本节根据前文确定的 V1 蒙皮几何尺寸 在 xoz 平面上使用 草图 绘制出符合尺寸的半月牙平面 在 xoy 平面上使用 草图 绘制出曲 率较小的弧线 使用 扫略 命令生成截面为半月牙的实体 其远离对称面部 天津科技大学 2013 届本科生毕业论文 8 分的实体生成方法与前者类似 文节不再赘述 左半边实体生成后 使用 凹 槽 命令去除模型中的通孔 完成左半边实体的建模后 使用 镜像 命令得 到整个 V1 蒙皮的实体模型 最后使用 抽壳 命令即得 V1 蒙皮的壳体模型 如图 3 2 所示 图 3 2 V1 蒙皮的几何模型 V1 泡沫的形状较简单 在 草图 中绘制中心平面和引导线后 通过 扫 略 命令获得沿一定弧度延伸的底面为四边形的长柱体 其后使用 凹槽 命 令生成其下表面的三处通槽 最后得到的 V1 泡沫几何模型如图 3 3 所示 图 3 3 V1 泡沫的几何模型 V1 横梁为左右对称的壳体 分为支架和梁两部分 梁部分 首先根据尺寸 数据通过在 草图 中绘制弓形曲线和引导线 在 线框与曲面设计 工作模 块中 扫略 获得 返回 零部件设计 工作模块使用 厚曲面 获得壳体模 型 使用 凹槽 命令生成通孔 支架部分 主要用到了 拉伸 命令生成实 体 使用 凹槽 命令生成通孔 最终得到的 V1 横梁的几何模型如图 3 4 所示 天津科技大学 2013 届本科生毕业论文 9 图 3 4 V1 横梁的几何模型 完成 V1 前保险杠系统的三大部件的几何模型后 在 CATIA 的 装配件设计 模块中 使用 现有组件 将三个部件模型导入其中 主要通过 相合约束 偏移约束 和 角度约束 等命令根据其位置及装配关系赋予各部件间的相 互约束 最终获得 V1 前保险杠系统的几何装配模型 如图 3 5 所示 图 3 5 V1 前保险杠系统的几何装配模型 3 1 2 3 V2 前保险杠系统几何模型的构建 本节应该用 CATIA 软件分别建立 V2 前保险杠系统三大组成部分 蒙皮 泡 沫和横梁 的几何模型 V2 前保险杠系统各部件几何模型的建模步骤及方法与 V1 前保险杠系统中对应的部件相类似 本节不再赘述 建成的 V2 蒙皮 泡沫 和横梁的几何模型分别如图 3 6 7 8 所示 图 3 6 V2 蒙皮的几何模型 天津科技大学 2013 届本科生毕业论文 10 图 3 7 V2 泡沫的几何模型 图 3 8 V2 横梁的几何模型 将 V2 蒙皮 泡沫和横梁的几何模型进行装配后得到 V2 前保险杠系统的几 何装配模型 如图 3 9 所示 图 3 9 V2 前保险杠系统的几何装配模型 3 2 汽车前保险杠有限元模型的构建 3 2 1 HYPERMESH 软件简介 HYPERMESH 软件是美国 Altair 公司研发的产品 其功能强大 集成了设计 与分析所需的各种工具 是世界领先的 CAE 应用软件包 也是一个创新和开发 天津科技大学 2013 届本科生毕业论文 11 的企业级 CAE 平台 在汽车安全领域中 主要应用于碰撞试验的前处理过程 图 3 10 所示为 HYPERMESH 10 0 的工作界面 图 3 10 HYPERMESH 10 0 界面 为了让各行业的用户在对各种模型的划分网格时都能有灵活全面的出色体 验 HYPERMESH 软件提供了三种划分网格的方法 其应用范围列于表 3 1 结合 运用三种网格划分方法能够解决各种模型的网格划分工作 本文只通过自动式 和网格优化两种方法划分网格 6 表 3 1 三种网格划分方法及其应用范围 网格划分方法应用范围 自动式 Auromatic 通过设定一定的网格参数对曲面多而复杂的几何体进行自 动划分网格 交互式 Interactive 主要应用于曲面不多且形状相对简单的几何体 利用这种 方法可以修改单元类型 曲面边界上的网格密度 网格平 滑度等参数 网格优化则是对已经划分好的网格作优化处 理 网格优化 QI optimized 对已经划分好的网格作优化处理 天津科技大学 2013 届本科生毕业论文 12 3 2 2 有限元网格划分的注意事项 3 2 2 1 网格划分 划分有限元网格需要到考虑所要研究的问题的规模 重点和周期 以及能 够为研究所用的计算机硬件等条件 力求计算时间 计算精度和计算成本的最 佳搭配 由于网格密度的大小会直接影响到能量吸收的情况 因此网格划分不 当将导致在受力梯度不大时就得到错误的结果 所以 在碰撞模拟试验前需要 仔细安排模型网格的疏密大小 尤其在变网格密度时更需要谨慎 应该要充分 地考虑到碰撞时模型的受力分布等情况 以避免由人为的单元划分而导致模拟 结构真实性受到不必要的影响 1 本文中有限元模型的整体网格尺寸控制在 10mm 左右 局部细节区域应该控制在 5mm 左右 另外 划分网格时还需注意以下几点 1 尽可能不使用过小的网格单元 以免网格密度过大而缩小时间步长 2 尽量使用四边形单元 减少三角形 四面体 棱柱单元的使用 3 前保险杠模型的每一个方向至少使用三个单元 尽量保证较均匀的网 格大小 4 为保证计算精度 避免集中出现锐角单元和翘曲的壳单元 1 3 2 2 2 单元特性 在 HYPERMESH 软件中 用户可根据分析需要 在其丰富的单元库中选择便 于试验计算与求解的单元类型 其中各主要单元类型的使用范围如表 3 2 所示 不同类型的单元有几十种不同的算法可供选择 本论文采用的单元类型主要包 括壳单元和实体单元 表 3 2 HYPERMESH 软件中的主要单元类型的使用范围 单元类型使用范围 壳单元 SHELL 薄板和曲面结构 体单元 SOLID 三维实体结构 面单元 PLANE 平面结构 梁单元 BEAM 需要承受弯矩的二维结构 杆单元 LINK 只承受轴向力的二维结构 3 2 2 3 网格质量的控制 为了尽量避免单元网格质量不好而导致计算步长的降低而浪费过多的求解 时间以及试验误差的增加 因此 对于一个大型有限元模型 往往需要很多时 天津科技大学 2013 届本科生毕业论文 13 间来检查模型单元的质量 单元质量检查主要包括三个方面 重合节点和重合 单元检查 自由边和自由面检查及单元形状检查 使用 check elems edges 和 faces 等功能进行单元网格质量的检查后 在工作界面上 HYPERMESH 软件会将模型上不合格的单元显示成红色 应用 translate smooth 等功能对其进行修正 如果修正后仍显示红色 则需要对该部件模 型重新划分网格 经检查和修正 本试验模型的所有单元均合格 9 3 3 3 有限元网格的划分 上一章中已经建立了 V1 和 V2 的 CATIA 前保险杠系统装配模型 确定了 V1 和 V2 的 CATIA 前保险杠系统中各部件的相对位置和装配关系 因此本章只需将 上一章所建立 V1 和 V2 的 CATIA 前保险杠系统装配模型 另存为 igs 格式 导入到 HYPERMESH 软件当中 进行网格划分 3 3 3 1 V1 前保险杠系统模型有限元网格的划分 本小节分别对 V1 各部件进行分析和网格划分 对于 V1 蒙皮模型 由于其为薄壁件 而且曲面片较多 轮廓相当复杂 因 此适合采用壳单元划分网格 此时采用壳单元可以提高网格质量及缩短工作时 间 网格划分过程中主要使用了 solids midsurfase quick edit 等 命令抽取中面并对其进行修正 其后使用 automesh 命令对中面划分 2D 网格 最后使用 check elems 等命令检查网格质量并优化网格 获得 V1 蒙皮有限 元模型 如图 3 11 所示 图 3 11 V1 蒙皮有限元网格模型 对于 V1 泡沫模型 由于其厚度较大 但轮廓十分简单 因此适合采用实体 单元划分网格 此时采用实体单元可以更能够提高网格质量及缩短工作时间 在开始划分前 为了避免 V1 蒙皮模型的干扰 本节使用 mark 命令将 V1 蒙 天津科技大学 2013 届本科生毕业论文 14 皮模型暂时隐藏 工作过程中首先通过 edges 命令检查了自由边等影响网格 质量的因素 其后通过 quick edit 等命令修复这些因素 然后使用 solids 命令建立实体 然后使用 automesh 命令对泡沫模型前表面进行 2D 网格的划分 其后使用 quick edit translate 等命令通过改变节点数 和移动节点等方式对 2D 网格进行优化 solid map 中的 general 命令将 上述划分的 2D 网格拉伸至泡沫模型后表面 完成泡沫模型 3D 网格的划分 获 得 V1 泡沫有限元网格模型 如图 3 12 所示 图 3 12 V1 泡沫有限元网格模型 V1 横梁模型的特征与 V1 蒙皮模型的类似 曲面片较多 轮廓相当复杂 为薄壁件 因此同样时候采用壳单元划分网格 同样的 首先使用 mark 命 令将 V1 蒙皮模型和 V1 泡沫模型暂时隐藏 工作过程中主要使用 midsurfase 命令抽取中面 使用 quick edit 等命令并对其进行修正抽 中面后产生的自由边 其后使用 automesh 命令对修正后的中面划分 2D 网格 最后使用 check elems 等命令检查网格质量并优化网格 获得 V1 横梁有限 元网格模型 如图 3 13 所示 图 3 13 V1 横梁有限元网格模型 天津科技大学 2013 届本科生毕业论文 15 最后 将隐藏的 V1 蒙皮模型和 V1 泡沫模型恢复显示 unmark all 检 查各部件的网格大小 质量 完后将 V1 前保险杠系统有限元网格模型 图 3 14 保存 save 待用 图 3 14 V1 前保险杠系统有限元网格模型 3 3 2 2 V2 前保险杠系统模型有限元网格的划分 本小节分别对 V2 各部件进行分析和网格划分 对于 V2 蒙皮模型 与 V1 蒙皮模型类似 由于其为薄壁件 而且曲面片较 多 轮廓相当复杂 因此也适合采用壳单元划分网格 此时采用壳单元可以提 高网格质量及缩短工作时间 工作过程中主要使用了 solids midsurfase 和 quick edit 等命令抽取中面并对其进行修正 其后使用 automesh 命令对中面划分 2D 网格 最后使用 check elems 等命令检查网 格质量并优化网格 获得 V2 蒙皮有限元网格模型 如图 3 15 所示 图 3 15 V2 蒙皮有限元网格模型 对于 V2 泡沫模型 其厚度较大但厚度均匀 上下端面 左右端面均为相互 天津科技大学 2013 届本科生毕业论文 16 平行的矩形 相比 V2 泡沫模型其轮廓更加简单 因此为了用尽可能少的工作时 间获得更佳的网格质量 同样采用实体单元划分网格 在开始划分前 为了避 免 V2 蒙皮模型的干扰 本节使用 mark 命令将 V2 蒙皮模型暂时隐藏 工作 过程中首先使用 edges 命令检查了自由边等影响网格质量的因素 其后通过 quick edit 等命令修复这些因素 然后使用 automesh 命令对泡沫模型 上端面进行 2D 网格的划分 其后使用 quick edit translate 等命令通 过改变节点数和移动节点等方式对 2D 网格进行优化 solid map 中的 general 命令将上述划分的 2D 网格拉伸至泡沫模型下端面 完成泡沫模型 3D 网格的划分 获得 V2 泡沫有限元网格模型 如图 3 16 所示 图 3 16 V2 泡沫有限元网格 相比 V1 横梁模型 V2 横梁模型的的厚度较大 不适合采用壳单元划分网 格 但其轮廓相对平整 可以使用实体单元划分网格 首先使用 mark 命令 将 V1 蒙皮模型和 V1 泡沫模型暂时隐藏 划分前先使用 edges 命令检查模型 的自由边并加以修正 为便于划分及获得质量较好的网格 划分过程中通过 solid edit 命令将模型 切割 成形状简单的 块 其后分别对这些 块 进行划分网格 其主要步骤仍是先用 automesh 命令对 块 的一平整表面 进行 2D 网格的划分 solid map 中的 general 及 ends only 等命令将 上述划分的 2D 网格进行沿特定路径的拉伸 完成该 块 的 3D 网格的划分 最后使用 check elems 等命令检查网格质量并优化网格 获得 V2 横梁有限 元网格模型 如图 3 17 所示 天津科技大学 2013 届本科生毕业论文 17 图 3 17 V2 横梁有限元网格模型 最后 将隐藏的 V2 蒙皮模型和 V2 泡沫模型恢复显示 unmark all 检查 各部件的网格大小 质量 完后将 V2 有限元网格模型 图 3 18 保存 save 待用 图 3 18 V1 前保险杠系统有限元网格模型 3 4 本章小结 本章参考某两款车型 V1 和 V2 前保险杠系统原型的尺寸及形状数据 通过 测量 比例计算等方法确定其各部件的尺寸构造 对原型结构加以修正和简化 后利用 CAD 软件绘制出其二维工程图 参照二维工程图 应用 CATIA 三维建模 软件分别建立了前保险杠系统的主要部件 蒙皮 泡沫和横梁 的简化几何模 型并装配为装配体模型 其后学习并利用 Hypermesh 软件对上述建立的前保险 杠几何模型进行有限元网格划分 最终分别获得 V1 和 V2 两款汽车的前保险杠 系统的有限元模型 完成了人 车碰撞试验中汽车部分的有限元建模 天津科技大学 2013 届本科生毕业论文 18 4 下肢 前保险杠计算机模拟碰撞试验 在目前的行人保护法规中 我国现行的行人保护法规是最贴合我国实际道 路交通现状的 因此 本章按照我国行人保护法规的要求 通过 PAM CRASH 软 件建立行人下肢撞击汽车前保险杠系统的碰撞有限元模型 考察上述两款汽车 前保险杠系统对行人下肢的保护性能 为了表达方便直观 本论文将下肢 V1 前保险杠碰撞仿真试验定义为试验 1 将下肢 V2 前保险杠碰撞仿真试验定义为 试验 2 4 1 PAM CRASH 软件简介 PAM CRASH 是 ESI 集团最新一代的碰撞模拟分析软件 是碰撞损伤评估和 安全性评估的必要工具 工程师们利用它可以减少或摆脱物理样机测试 可以 将更多的时间和精力投入到创建 评估和管理在真实碰撞环境下的模拟样机测 试中 工程师还可以相当方便地改变碰撞环境 从而得到不同碰撞条件下的仿 真结果 PAM CRASH 是基于显示有限元算法的一种计算机三维碰撞冲击仿真模拟系 统 是一个广泛应用于汽车 航空 电子和材料生产等专业领域 有关碰撞 冲击和安全性模拟测试等问题的解决的实用工具 PAM CRASH 能够对大位移 大旋转 大应变 接触碰撞等问题进行非常精确的模拟仿真 能够非常简便的 处理异常复杂的边界约束 其主要显著特点如表 4 1 所示 表 4 1 PAM CRASH 软件的主要显著特点 主要显著特点 三维图形显示属性灵活控制 色彩多样逼真 支持多 CPU 并行计算 DMP 和 SMP 运算效率高 能够简便地处理异常复杂的边界约束 灵活控制计算的时间步长 动态分配内存 无须用户设置 灵活搜寻接触区间 针对大变形材料可采用特有的措施来保证求解的稳定性和精确性 可设定材料的断裂失效条件 简便地定义焊点 铆钉等约束及其断裂条件 可设置阻尼以加快求解弹性接触时的收敛 天津科技大学 2013 届本科生毕业论文 19 针对汽车碰撞而特设指标整形 输出 比较模块 4 2 仿真试验的前处理 4 2 1 基本假设 由于受到时间 资源和经费等多种因素的限制 本论文的人 车碰撞试验中 只对已构建汽车的前保险杠部分模型进行运动力学仿真 本试验主要描述人 车碰撞时汽车的前保险杠系统的结构对行人安全保护性 能的影响 在保证一定精确度的前提下 忽略下肢模型中足掌部与地面之间间 的摩擦系数 由于本文使用计算机软件仿真试验 可以在软件中添加所需要的驱动力 飞行距离等约束条件 因此 本研究不需要设置驱动装置 因此无须考虑下肢 有限元模型与驱动装置在碰撞过程中的接触对试验求解结果的影响 1 4 2 2 参数量纲 在 PAM CRASH 软件中 没有明确的规定所输参数的单位 而是通过各个物 理量单位之间的关系而自动匹配的 因此 为了避免得到错误的计算求解结果 使用的参数量纲必须是相互匹配的 这点尤其重要 本仿真试验中主要使用的 参数量纲组合如表 4 2 所示 表 4 2 常用的量纲组合 基本物理量导出物理量 序号时间长度质量密度重力加速度力应力能量 1s mKgKg m39 81mm ms2N PaN m 2smm t t mm39810mm s2kN MPaN mm 4 2 3 材料属性 有限元分析结果的准确性在很大程度上依赖于材料模型的选择以及材料参 数的确定 花费相当充足的时间了解材料模型并获得准确的材料参数尤其必要 在 PAM CRASH 求解前 需要定义模型的密度 泊松比 弹性模量以及屈服应力 等参数 其中 V1 和 V2 的蒙皮均为 PA 材料 泡沫均使用 EPP 材料 V1 横梁的 材料为 45 钢 V2 横梁为ZAlSi7Mg制造 V1 和 V2 的前保险杠系统各主要部件 的材料参数见表 4 3 和表 4 4 在定义材料时应该注意以下两点 1 并不是材料属性输入窗口中的所有填选项都必须要输入数值 部分不 天津科技大学 2013 届本科生毕业论文 20 影响试验求解的参数可以缺省 2 务必保证所输入材料参数单位的统一 错误的单位制不但会影响到材 料的响应 还会影响到接触刚度等的计算 1 表 4 3 V1 前保险杠系统主要部件的材料参数 参数密度 kg m3 泊松 比 弹性模量 Pa 屈服应力 Pa 切线模量 Pa 失效应 变 蒙皮 PA 22000 3200E6100E6 泡沫 EPP 91 30 330 6E6 横梁 45 78500 3210E91034E62435E60 25 表 4 4 V2 前保险杠系统主要部件的材料参数 参数密度 kg m3 泊松 比 弹性模量 Pa 屈服应力 Pa 切线模量 Pa 失效应 变 蒙皮 PA 22000 3200E6100E6 泡沫 EPP 91 30 330 6E6 横梁 ZAlSi7Mg 27000 372E9345E6690E60 17 4 2 4 厚度参数 V1 的蒙皮和横梁以及 V2 的蒙皮有限元模型均使用壳单元划分网格 在前 处理过程中需要将其厚度参数通过输入定义到壳单元模型中 各模型的厚度参 数见表 4 5 表 4 5 相关模型的厚度参数 模型厚度 mm V1 蒙皮 4 5 V1 横梁 1 47 V2 蒙皮 4 5 4 2 5 边界条件 一项来自国际协调研究机构 IHRA 的调研报告显示 在汽车与行人相碰 天津科技大学 2013 届本科生毕业论文 21 撞的交通事故中 70 的碰撞车速都在 40km h 以下 而且在此车速下的行人的 死亡率也高达 26 2 因此 40km h 为比较理想的并贴近现实情况的碰撞车速 根据试验要求 定义下肢有限元模型的初始速度为 40km h 选取在事故中最易 造成下肢损伤的位置 即车辆中心轴线处 作为碰撞点 由于两款汽车 V1 和 V2 的车身尺寸不同 碰撞点分别处于距离地面 512mm 和 438mm 的位置 4 2 6 接触定义 显示动力分析软件 PAM CRASH 有三种接触面处理算法 单面接触 点面接 触和面面接触 单面接触无需定义接触面和目标面 这种算法十分适用于事先 不知道接触方案的情况 若预先已经指导接触面很小 可以采用点面接触 而 面面接触则适用于一个物体的面穿透另一个物体的面的情况 在软件中 相应 的可用于处理动态接触 碰撞问题的接触面处理算法也有三种 对称罚函数 法 动力约束法和分配参数法 对称罚函数法是软件默认的接触算法 很少激 起网格的灵能模式 不会产生噪音 显示动力分析软件最早采用的算法是动力 约束法 此算法比较复杂 只适用于固连接触问题 而分配参数法则只适用于 滑动界面接触问题的处理 4 本碰撞仿真试验采用的接触包括自动单面接触和 面面接触 模型接触的定义见表 4 6 表 4 6 人体下肢 前保险杠系统碰撞模型的接触定义 接触模型接触类型摩擦系数 1 前保险杠模型全部部件自接触单面接触 0 2 2 下肢模型自身接触自动接触 0 2 3 下肢模型与前保险杠接触面面接触 0 4 4 下肢肌肉与骨骼接触面面接触 0 8 4 2 7 求解时间 由前人的仿真试验可发现 汽车与行人在碰撞 30ms 后 力学指标基本稳定 因此 为了准确地得到整个碰撞系统的力学指标 本论文选取仿真的时间为 30ms 4 2 8 碰撞模型 完成碰撞仿真试验模型的前处理过程后 得到试验 1 和试验 2 的仿真试验 起点模型分别如图 4 1 图 4 2 所示 天津科技大学 2013 届本科生毕业论文 22 图 4 1 试验 1 的仿真起点模型 图 4 2 试验 2 的仿真起点模型 4 3 碰撞模拟试验求解 将上述确定的相关参数输入到 PAM CRASH 软件参数输入窗口中 运行求解 以 V1 横梁的参数输入为例 如图 4 3 所示 将材料编号输入到 IDMAT 将密 度输入到 RHO 将弹性模量输入到 E 将泊松比输入到 NU 将屈服应力 输入到 SIGMAy 等 天津科技大学 2013 届本科生毕业论文 23 图 4 3 V1 横梁的参数输入窗口 4 4 仿真结果的后处理 仿真时间终点 30ms 的模型 分别如图 4 4 图 4 5 所示 通过 PAM CRASH 软件中的 VISUAL VIEWER 模块获得的胫骨加速度 时间曲线 膝关节剪切 位移 时间曲线 膝关节弯曲角度 时间曲线 分别如图 4 5 和图 4 6 所示 图 4 4 试验 1 的仿真终点模型 天津科技大学 2013 届本科生毕业论文 24 图 4 5 试验 2 的仿真终点模型 a 胫骨加速度 时间曲线 b 膝关节弯曲角度 时间曲线 c 膝关节剪切位移 时间曲线 图 4 6 试验 1 的后处理结果 天津科技大学 2013 届本科生毕业论文 25 a 胫骨加速度 时间曲线 b 膝关节弯曲角度 时间曲线 c 膝关节剪切位移 时间曲线 图 4 7 试验 2 的后处理结果 天津科技大学 2013 届本科生毕业论文 26 另外 为了更直观地观察在碰撞过程中下肢骨骼的变形情况 本文将下肢 的皮肤和肌肉隐藏 通过动画截图 获得试验 1 和试验 2 中下肢模型在不同时 刻的变形 分别如图 4 8 和图 4 9 所示 图 4 8 试验 1 中下肢模型在不同时刻的变形 图 4 9 试验 2 中下肢模型在不同时刻的变形 4 5 仿真结果的分析与讨论 4 5 1 独立分析 分别对两款汽车的碰撞试验结果进行以下独立分析 1 对于试验 1 由图 4 6 a 可得胫骨峰值加速度为 145 4g 小于法 规所规定的 170g 从图 4 6 b 可见 碰撞的 0 30ms 间 最大的膝关节弯曲 角度为 15 5 没超出法规所规定的 19 同时由图 4 6 c 得出最大的膝 关节剪切位移为 4 66mm 符合法规中要求小于 6mm 的规定 并从图 4 8 可以直 观的看到 从碰撞开始到 30ms 期间 下肢的变形较小 而且没有发生明显的骨 折和韧带损伤 因此 可以论定 V1 前保险杠系统符合行人保护法规的规定 在 碰撞事故发生时能够为行人提供有效的保护 2 对于试验 2 由图 4 7 a 可得胫骨峰值加速度为 135 2g 小于法 天津科技大学 2013 届本科生毕业论文 27 规所规定的 170g 从图 4 7 b 可见 在碰撞 20ms 左右膝关节弯曲角度开始 超出法规所规定的 19 并继续增大 同时由